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ABSTRACT 


This  report  summarizes  results  of  studies  which  were  conducted  at 
North  Carolina  State  College  under  various  contracts  with  the  Department 

f 

of, the  Army  relative  to  sources  of  dispersion  of  artillery-type  rockets. 
Mathematical  equations  and  formulas  which  are  applicable  to  the  analysis 

of  the  of  various  dispersion-producing  factors  on  both  spin- 

.  *  .*• 

stabilized  and  fin-stabilized  gun-boosted  rockets  are  presented,  and 

*  * 

their  uses ’are  illustrated  by  numerical  examples.  Experimental  techniques 
are  described  and  some  results  are  given.  Implications  of  mathematical 
results  relative  to  rocket  design  are  discussed.  A  mathematical  model 
for  studying  effect  of  launcher  tube  motion  is  given  in  Appendix  B. 
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CHAPTER  1 


INTRODUCTION 

The  purpose  of  this  report  is  to  summarize  the  work  of  the  Rocket  Research 
Group  at  North  Carolina  State  College  for  the  total  period  of  duration  of  its 
investigations,  and  to  present  the  results  obtained  in  a  unified  and  related 
manner  so  that  they  may  be  of  use  in  the  design  and  development  of  new 
artillery-type  rockets.  This  work  was  done  for  Army  Ordnance  under  Con¬ 
tract  Nos.  DA -01-021  0RD-3190,  DA -01-021  0RD-U592,  DA-36-03U-509  OED-25, 
and  DA-01-009  0RD-1022,  and  it  spans  the  period  from  June,  1952  to  August, 

1962.  This  research  group  had  as  its  basic  objectives  the  following: 

(1)  To  investigate  sources  of  dispersion  of  spin-stabilized  and  fin- 
stabilized  rockets  by  theoretical,  computational  and  experimental 
methods . 

(2)  To  assist  in  the  design  and  development  of  experimental  methods 
of  measuring  certain  parameters  that  are  of  significance  in  the 
study  of  rocket  accuracy. 

(3)  To  act  in  a  consulting  capacity  for  other  groups  who  might  be 
working  on  the  design  and  development  of  artillery-type  rockets. 

In  initiating  the  theoretical  studies  of  sources  of  dispersion  of  spin- 
stabilized  rockets,  this  group  first  made  an  exhaustive  study  of  the  existing 
literature  on  the  theoretical  treatment  of  the  motion  of  spin-stabilized 
rockets  in  order  to  determine  the  most  suitable  existing  mathematical  model 
which  might  be  used  as  a  springboard  for  its  activities.  The  theoretical 
treatments  given  by  Davis  [DFB]*,  Follin  [F],  Galbraith  [Ga],  Harrington  [HI], 
and  Herz  [H],  and  Rankin  [Ra]  were  compared  and  found  to  be  essentially 
equivalent  under  uniform  assumptions  relative  to  the  various  physical  quan¬ 
tities  associated  with  the  motion.  This  group  consequently  adopted  the 
notation  and  approach  used  by  Harrington  as  a  takeoff  point  for  further 
study  and  development  of  the  theory  of  rocket  motion. 


Bracketed  expressions  refer  to  the  list  of  references  at  the  end  of  the  report. 


CHAPTER  2 


FACTORS  WHICH  CONTRIBUTE  TO  ROCKET  INACCURACY 

Before  going  into  the  theoretical  basis  for  the  study  of  rocket  motion, 
we  shall  examine  first  certain  factors  which  may  lead  to  rocket  inaccuracy. 
In  later  chapters  mathematical  expressions,  called  characteristic  functions, 
will  be  introduced  by  means  of  which  the  effect  on  angular  deviation  for 
a  unit  amount  of  any  one  of  these  disturbing  factors  may  be  computed. 

All  of  the  factors  to  be  discussed  here  contribute  in  varying  degrees 
to  the  inaccuracy  of  both  spin-stabilized  and  fin-stabilized  artillery-type 
rockets  except  that  of  fin  misalignment,  which  clearly  applies  only  to  the 
latter  type  of  rocket. 

In  the  course  of  the  following  discussion  it  will  be  necessary  to 
frequently  make  use  of  the  term  ge ometric  axis  of  the  rocket.  Although 
the  rocket  would  have  to  be  perfectly  formed  for  a  geometric  axis  to  exist, 
it  suffices  for  practical  purposes  to  define  such  an  axis  as  the  line  of 
centers  of  two  circular  bands  or  bourrelets,  one  placed  around  the  rocket 
near  the  rear  of  the  rocket  and  the  other  near  the  forward  part  of  the 
rocket  near  the  point  where  the  body  of  the  rocket  begins  to  taper  off 
toward  the  nose. 

2.1.  Initial  Cross -Spin, 

At  the  instant  the  rear  end  of  the  rocket  emerges  from  the  launcher, 
the  rocket  is  describing  a  transverse  rotation  about  an  axis  passing  through 
the  center  of  gravity  and  perpendicular  to  the  axis  of  the  rocket.  Thus  the 
angle  ($)  which  measures  the  direction  in  which  the  rocket  is  pointing  with 
reference  to  some  fixed  direction  (such  as  the  bore-line  of  the  launcher) 
is  undergoing  a  time  rate  of  change  at  launch  which  is  denoted  by  (the 
subscript  zero  being  used  to  indicate  here  a  value  at  launch).  This  cross¬ 
spin  (or  transverse  angular  velocity)  has  long  been  recognized  as  one  of  the 
significant  contributors  to  deviations  of  the  rocket  from  the  desired 
direction  of  flight. 

The  cross-spin  at  launch  may  be  attributed  to  various  causes,  one  of 
which  is  certainly  the  tip-off  effect  induced  by  gravity  because  of  the 
fact  that  for  a  short  interval  during  the  launching  phase  the  center 
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of  gravity  of  the  rocket  is  already  outside'  the  launcher  while  the  rear  end 
of  the  rocket  is  still  engaged.  Other  causes,  are,  attributable  to  behavior 
of  the  rocket  within  the  launcher  itself,  where  precessional  motion  of  the 
center  of  the  rear  band  of  the  rocket,  static  and  dynamic  imbalance,  and 
misaligned  thrust  (see  later  sections)  can  cause  a  build-up  through  the 
tip-off  phase  of  a  significant  amount’  of  cross-spin  at  launch. 

It  is  highly  probable  that  launcher  reaction  also  contributes  signifi¬ 
cantly  to  cross-spin  at  launch,  especially -in  the  case  of  gun -boosted 
rockets.  Some  preliminary  experimental :  work  in  this  area  indicates  that 
this  is  true,  but  definitive  results  still  have  not  been  obtained. 

2.2.  Initial  Yaw,  A  . 

- - - — o 

The  yaw  of  a  rocket  is  defined  as  the  angle  which  the  geometric  axis  of 
the  rocket  makes  with  the  tangent  to  the  curve  described  by  the  center  of 
gravity  of  the  rocket  during  flight.  This  curve  is  usually  referred  to  as 
the  trajectory  of  the  rocket.  It  is  clear  that  the  larger  this  yaw  angle 
becomes,  the  more  tendency  there  is  for  the  rocket  to  be  driven  off  the 
desired  flight  path  by  the  thrust  imparted  by' the  rocket  jets.  Thus  it  is 
that  any  yaw  which  the  rocket  attains  during  the  tip-off  period,  and  hence 
has  acquired  at  the  instant  the  rocket  becomes  disengaged  from  the  launcher, 
will  be  a  factor  in  determining  the  direction  of  flight  of  the  rocket  through¬ 
out  the  burning  period.  It  is  true  that  under  .desirable  conditions  the 
varying  yaw  outside  the  launched  may  damp  out  rather  rapidly,  but  its  value 
at  launch  still  has  an  effect,  that  endures  to  the  end  of  burning . 

The  same  factors  that  contribute  to  initial  cross -spin,  as  discussed 
in  Section  2.1  above,  can  cause  yaw  at  launch. 

2.3.  Cross-Wind*  w  .  ■  ■■  - 

- r  — C 

The  component  of  the  wind  velocity  which  is  perpendicular  to  the  tra- 

■  1:  .  *  "  ‘ 

jectory  of  the  rocket  is  a  significant  'factor  in.  causing  rocket,  inaccuracy, 
especially  if  the  winds  are  quite  gusty.  ■•It  will  be  seen  later  that  a  con¬ 
stant  cross-wind  has  an  effect  equivalent,  to-'an  initial  yaw.  It  will  also 
be  shown  that  the  major  contribution  of  cross-wind,  to  deviation  of  the  rocket 


from  the  desired  direction  of  motion  occurs  in  the  interval  immediately 
following  the  launching  of  the  rocket  pa^icbi^ly  if  the  launch  velocity 
is  not  large.  Another  significant  result  that '  shows' up  in  the  theory iis 
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that  wind  effect  may  be  diminished  by  decreasing  the  variation  in  the  velocity 
during  the  flight  of  the  rocket. 

2.1;.  Dynamic  Unbalance , 

Dynamic  unbalance  arises  from  the  fact  that  the  rocket  jets  are  attempt¬ 
ing  to  rotate  the  rocket  about  one  axis  (the  geometric  axis  for  a  perfectly 
symmetrical  rocket  with  perfectly  aligned  jets),  while,  due  to  the  unsymme- 
trical  distribution  of  the  mass  of  the  rocket,  the  inclination  of  the  rocket 
is  to  rotate  about  another  axis  (the  longitudinal  axis  of  inertia).  The 

angle,  (3  ,  between  these  two  axes  is  taken  as  the  measure  of  dynamic  unbalance, 
c 

The  dynamic  unbalance  of  an  unloaded  round  may  be  measured  on  a  machine 
designed  for  this  purpose,  such  as  the  Tinius-Olsen  balancing  machine,  and 
the  amount  of  unbalance  may  range  in  magnitude  up  to  several  mils  (or 
thousandths  of  a  radian).  This  magnitude  depends  to  a  great  extent  on  the l 
tolerances  maintained  in  manufacturing  and  assembling  the  metal  parts  of  i 
the  rocket. 

Clearly  as  t  he  rate  of  spin  of  the  rocket  is  increased  the  effect  of 

a  given  amount  of  dynamic  unbalance  on  the  direction  of  motion  of  the  rocket 

.  \ 

will  be  increased.  A  high  rate  of  spin  may  also  cause  a  break-up  of  the 
propellant  in  the  rocket  motor  and  thus  cause  an  increase  in  the  amount  of 
dynamic  unbalance  of  the  rocket  as  the  burning  progresses. 

2.5.  Static  Unbalance,  r  . 

- '  — c 

Ideally,  the  mass  distribution  of  the  rocket  should  be  such  that  the 

center  of  gravity  (henceforth  referred  to  as  the  c.g.)  is  on  the  geometric 

axis  of  the  rocket.  If  this  fails  to  be  true,  the  rocket  is  said  to  have 

static  unbalance,  which  is  measured  by  the  distance,  r  ,  from  the  c.g.  to 

c 

the  geometric  axis. 

Static  unbalance  is  brought  about  by  the  same  factors  that  cause  dynamic 
unbalance,  as  listed  in  Section  2.U  above.  The  chief  contribution  of  static 
unbalance  to  rocket  inaccuracy  is  made  through  its  effect  on  the  behavior 
of  the  rocket  during  the  launching  (or  tip-off)  phase,  and  thus  shows  up  in 
the  form  of  additions  to  and  Aq,  the  quantities  discussed  above  in 
Sections  2.1  and  2.2. 

2.6.  Linear  Thrust  Misalignment.  L^. 

If  the  resultant  line  o£  thrust  due  to  the  rocket  jets  fails  to  pass 
through  the  c.g.  of  the  rocket,  the  rocket  is  said  to  have  linear  thrust 
misalignment,  which  is  measured  by  the  distance,  Lq,  from  the  c.g.  to  this 
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line  of  thrust.  When  the  rocket  has  such  a  misaligned  thrust,  it  is  clear 
that  there  is  a  tendency  of  the  thrust  to  rotate  the  rocket  about  a  line 
through  c.g.  and  perpendicular  to  the  axis  of  the  rocket,  and  this  torque 
can  cause  significant  deviations  of  the  rocket  from  the  desired  direction 
of  flight. 

Linear  thrust  misalignment  can  be  caused  by  lack  of  uniformity  in 
the  manufacture  of  the  jet  nozzles,  by  improper  insertion  of  these  nozzles, 
or  by  tilting  of  the  nozzle  plate  in  assembling  the  rocket.  This  type 
of  disturbance  could  also  be  caused  by  erosion  of  the  jet  nozzles  during 
burning,  by  the  improper  attachment  of  a  diverter  plate,  or  by  the  plugging 
of  nozzles  by  pieces  of  propellant  during  burning. 

In  cases  where  the  rocket  spins,  the  effects  of  misaligned  thrust  tend 
to  cancel  out,  and  hence  this  factor  may  then  become  relatively  negligible 
in  cases  where  reasonable  tolerances  are  maintained  in  the  manufacturing 
and  assembling  of  the  rocket.  For  the  purpose  of  counteracting  this  effect, 
it  is  thus  desirable  to  impart  some  spin  to  fin-stabilized  rockets. 

2.7.  Angular  Thrust  Misalignment ,  a  . 

The  angle,  a  ,  which  the  resultant  line  of  thrust  of  the  rocket  jets 
c 

makes  with  the  geometric  axis  of  the  rocket  is  called  the  angular  thrust 
misalignment  of  the  rocket.  Such  misalignment  can  be  attributed  to  the 
same  causes  as  those  listed  above  in  Section  2.6.  as  .causes,  of  linear,. thrust 
misalignment. 

Practically  speaking,  the  effect  of  angular  thrust  misalignment  on 
the  motion  of  spin-stabilized  rockets  is  negligible,  but  it  may  be  signifi¬ 
cant  in  the  case  of  fin-stabilized  rockets. 

2.8.  Fin  Misalignment,  p^. 

If  a  perfectly  made  fin-stabilized  rocket  is  suspended  in  an  air 
stream  (such  as  in  a  wind  tunnel)  which  initially  is  parallel  to  the  geometric 
axis  of  the  rocket  and  in  such  a  way  that  the  rocket  is  free  to  rotate  about 
the  c.g.,  then  the  air  stream  will  cause  no  change  in  the  direction  of 
the  geometric  axis.  However,  if  there  are  bent  fins  or  other  imperfections 
in  the  fin  assembly,  the  rocket,  under  the  above  conditions,,  may  assume  a 
position  in  which  its  axis  makes  an  angle,  p  ,  relative  to  the  direction, 
of  the  air  stream.  This  angle  is  taken  as  the  measure  of  fin  misalignment . 
Thus  a  rocket  with  misaligned  fins  has  a  natural  tendency  to  assume  a  yaw 
relative  to  the  resultant  direction  of  air  flow  during  its  flight,  and  this 
can  be  a  significant  source  of  inaccuracy  for  fin-stabilized  rockets. 


CHAPTER  3 


MATHEMATICAL  BASIS  FOR  STUDY  OF  THE  MOTION  OF 
A  SPIN-STABILIZED  ROCKET  DURING  BURNING  OUTSIDE  THE  LAUNCHER 

3.1.  Reference  Systems 

■  In  setting  up  the  differential  equations  which  describe  the  motion  of  a 
rocket,-  several  coordinate  systems  are  used.  For  complete  details  as  to  all 
of  these  coordinate  systems,  reference  is  made  to  lH-1]  or  [CH],  For  the 
purposes  of  the  present  report,  it  is  sufficient  to  say  that  the  location  ' 
of  the  center  of  gravity  of  the  rocket  in  space,  and  the  rotational  motion 
of  the  rocket  about  the  c.g.  are  completely  described  by  use  of  the  three 
rectangular  coordinates  X^,  Yq,  Zq  of  the  center  of  gravity,  and  by  the 
three  complex  angles,  #,  0,  A,  denoting  respectively  orientation,  angular 
deviation,  and  yaw  of  the  rocket.  All  of  these  quantities  will  be  described 
more  fully  in  what  follows. 

The  coordinate  system  to  which  Xq,  Yq,  Zq  are  referred  is  a  right-handed 
rectangular  coordinate  system  0q-  with  the  X^-axis  along  the  bore¬ 

line  of  the  launcher  and  Xq  measured  positively  in  the  direction  of  travel 
of  the-  rocket,  the  Y^-axis  running  vertically  with  Yq  measured  positively 
upward,  and  the  ZQ-axis  running  horizontally  with  Zq  measured  positively 
to  the  right  (as  viewed  from  the  rear  of  the  launcher).  The  origin  0q  of 
this  system  is  fixed  on  the  boreline  in  such  a  position  that  the  X^-coordinate 

of  a  point  (denoted  by  X  )  at  the  muzzle  of  the  launcher  is  given  by 
2  ox. 

X  =.v  /(2G)  (See  Fig„3.1)where  v  is  the  launch  velocity  (velocity  at 
ox  °  o 

the  end  of  tip-off)  of  the  rocket,  and  G  is  the  acceleration  of  the  rocket 
outside  the  launcher.  It  is  assumed  here,  as  in  much  of  the  later  work  of 
this  chapter,  that  G  is  constant  during  burning.  In  actual  practice,  the 
average  value  of  G  over  this  interval  will  serve  as.  an  adequate  approximation 
for  use  in  this  mathematical  model.  It  should  be  noted  that  for  rockets 

fired  from  an  open  launcher,  X  can  be  interpreted  as  the  length  of  the 

ox 

launcher,  but  for  boosted  rockets,  X  varies  with  the  amount  of  boost  and 

.  ox 

is  thus  often  called' the  effective  launcher  length. 

Since- Xq  measures  distance  down  range  along  the  boreline,  one  sees  that 

for  the  interval  immediately  following  launch,  where  there  is  little  bending 

of  the  trajectory,  the  XQ-coordinate  of  the  c.g.  during  burning  is  given 


COORDINATE  SYSTEM 
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by 


x 

o 


=  X 

01 


4  ^-v2 


+  V  (t-t  ), 
0  0 


where  (t-tQ)  represents  the  time  elapsed  since  the  c.g.  left  the  muzzle  of 
the  launcher.  The  other  two  space  coordinates,  Y  ,  Zq  of  the  center  of 
gravity  (denoted  by  c.g.)  are  dealt  with  most  conveniently  as  components 
of  a  complex  number  R,  where 

R  =  Y  +  iZ 
o  o 

This  quantity  is  called  the  linear  deviation  of  the  rocket,  and  it  clearly 
determines  how  far  and  in  what  direction  the  c.g.  deviates  from  the  bore¬ 
line  of  the  launcher.  Later  in  this  report,  in  order  to  take  care  of 
situations  where  there  is  appreciable  change  in  the  direction  of  the  tra¬ 
jectory  from  launch  to  burnout,  the  meaning  of  R  will  be  extended  to  denote 
the  linear  deviation  of  the  rocket  from  an  ideal  trajectory  that  the  center 
of  gravity  would  traverse  under  the  actions  of  jet  thrust,  aerodynamic  drag, 
and  gravitational  force,  with  all  other  forces  neglected.  In  general,  we 
shall  thus  write 

R  =  Ry  +  iRz, 

so  that  Ry  denotes  the  upward  or  downward  (negative)  displacement  (in  ft. ), 

and  R^  the  horizontal  displacement  of  the  c.g.  from  the  ideal  trajectory. 

This  means  that  the  reference  line  for  R  in  this  case  is  the  tangent  to  the 

trajectory  rather  than  the  boreline  OX.  If  the  trajectory  is  essentially 

o  o 

a  straight  line,  clearly  Ry  ~  Yq  and  Rz  -  Zq.  The  magnitude  of  R  is  clearly 
given  by 

|R  1  =  (Ry2  +  Rz2)1/2. 


Thus,  for  example,  if  one  has  given  that 


R  =  -  O.li  +  i  0.3 

for  a  certain  point  on  the  trajectory.,  this  means  that  the  c.g.  of  the 
rocket  is  0.1).  ft.  downward  (measured  perpendicular  to  the  tangent)  and  0.3 
ft.  to  the  right  of  the  desired  path  of  flight  (the  ideal  trajectory).  The 


The  symbol  =  is  used  to  denote  "approximately  equal  to.  " 
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magnitude  of  this  linear  deviation  is 

!r|  =  [(-o.ii)2  +  (0.3)2  1//2  =  0.5  ft., 

so  that  at  the  point  in  question  the  rocket  is  actually  a  distance  of  0.5  ft. 
from  the  desired  trajectory. 

In  order  to  study  the  rotational  motion  of  the  rocket  about  its  c.g., 
as  well  as  to  deal  with  the- direction  of  motion  of  the  rocket  during  burning, 
a  second  right-handed  rectangular  coordinate  system  0-XYZ  is  used,  with  0 
at  the  c.g.  of  the  rocket,  and  with  OX,  OY,  OZ  fixed  in  direction  parallel 
respectively  to  O^X^,  O^Y^,  and  0qZq  as  given  above.  Also  the  complex  angles 
9,  A,  and# are  introduced.  The  complex  angle  0  is  defined  as  the  angle  bet¬ 
ween  the  forward  tangent  to  the  actual  trajectory  of  the  rocket  and  a  re¬ 
ference  line  which  at  first  will  be  taken  as  the  boreline  OX  of  the  rocket 
launcher.  This  angle  is  measured  positively  from  the  positive  reference 
line  to  the  forward  tangent  to  the  trajectory.  This  angle  will  be  referred 
to  as  the  angular  deviation  of  the  rocket,  and  clearly  defines  the  direction 
of  motion  of  the  c.g.  of  the  rocket.  Later,  in  order  to  take  care  of  cases 
where  there  is  appreciable  change  in  direction  of  the  trajectory  during 
burning,  the  reference  line  for  measuring  0  will  be  taken  as  the  forward 
tangent  to  the  ideal  trajectory. 

The  angle  ©  is  used  as  a  measure  of  rocket  accuracy,  since  a  large 
magnitude  for  ©  at  burnout  (denoted  by  indicates  a  large  deviation  of 
the  rocket  from  the  desired  trajectory,  and  this  deviation  will  in  turn  be 
reflected  in  the  deviation  of  the  impact  point  from  the  desired  target 
point.  Furthermore,  lack  of  reproducibility  of  ©^  from  round  to  round  in 
a  given  firing  of  successive  rounds  will  be  reflected  in  the  dispersion  of 
the  impact  points  for  this  group  of  rounds. 

It  is  convenient  to  represent  0  as  a  complex  angle,  where,  as  in 
dealing  with  linear  deviation,  R,  above,  the  term  complex  refers  to  the 
familiar  complex  numbers  of  the  form  a  +  i  b  which  are  dealt  with  in  algebra. 
One  can  then  embody  in  one  expression  for  the  angle  ©  both  the  size  of  the 
angle,  or  magnitude  (denoted  by  J ©|  ),  and  the  location,  or  orientation,  of 
the  plane  in  which  the  angle  is  measured.  This  plane  rotates  as  the  rocket 
moves  down  range.  Thus  we  shall  write,  in  general, 

©  =  ^  +  i©2. 
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where  the  components  ®y  and  0^  are  in  radians,  0^  being  measured  positively 
upward  in  the  vertical  plane  through  the  reference  line,  and  0  measured 

L4 

positively  to  the  right  in  the  plane  (through  the  reference  line)  perpendi¬ 
cular  to  that  of  0y. 

The  magnitude  of  the  complex  angle  9  is  then  given  by 

M  =  (Qy2  +  ®Z2)1/2, 

and  the  orientation  of  the  plane  (called  argument  of  9)  in  which  this 
magnitude  is  measured  is  given  by 


Arg  3  =  arc 


tan  ®z/®yo 


This  latter  angle  is  measured  clockwise  from  the  vertical  as  one  faces  down 
range. 

For  example,  if 

0  =  -  0.0070  +  i  0.0050 


at  a  certain  point  on  the  rocket  trajectory,  and  if  the  boreline  of  the 
rocket  launcher  serves  as  references  line,  then  at  the  point  in  question 
the  forward  tangent  to  the  trajectory  is  pointing  0.00^0  rad. (approximately 
5  mils)  to  the  right  and  0.0070  rad.  (or  7  mils)  downward  from  the  boreline 
of  the  launcher  (or  from  the  desired  direction  of  motion).  Otherwise  put, 
we  compute  the  magnitude  of  3, 


0 


(-  0.0070)2  +  (0.0050)2 


1/2 


=  0.0086  rad. 


and  the  orientation  of  0, 


Arg  0  =  arc  tan  (O.OCgO/-  0.0070)  =  H4*.l*° 

(the  choice  of  angle  is  clear  from  the  description  given  by  the  components 
above,  namely,  wto  the  right  and  downward")  and  get  the  result  that  the 
tangent  to  the  trajectory  at  the  point  in  question  lies  in  a  plane  passing 
through  the  reference  line  and  rotated  114*. ij0  from  the  vertical  in  a  clock¬ 
wise  sense.  Furthermore,  the  magnitude  of  0,  or  the  actual  angle  which  the 
tangent  to  the  trajectory  makes  with  the  reference  line,  is  ,0.0086  rad. (or 
8.6  mils).  We  thus  say  that  the  angular  deviation  (or  direction  of  motion) 
of  the  rocket  at  this  point  is  8.6  mils,  this  deviation  being  in  a  direction 
indicated  by  Arg  0  =  li*l**U°j  as  already  described* 
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The  complex  angle  <f>,  called  the  orientation  of  the  rocket,  is  the  angle 
determined  by  the  geometric  axis  of  the  rocket  and  the  reference  line  used 
in  measuring  0  .  This  angle  is  measured  positively  from  the  positive  reference 
line  to  the  forward  directed  rocket  axis.  (The  term  geometric  axis,  also 
called  the  bourrelet  axis,  is  used  to  mean  the  axis  of  the  cylindrical  body 
of  the  rocket,  assuming  that  it  is  a  perfect  circular  cylinder).  This 
angle  gives  the  direction  in  which  the  rocket  is  pointing,  and  is  repre¬ 
sented  in  complex  form  by 


<£  =  4> 


i 


‘Y  :  *  V 

so  that  the  magnitude  of  #  is  given  by 

I  #  !  =  ($Y2  + 

and  the  location  of  the  plane  in  which  #  is  measured  by 


2)1/2, 


Arg  <£  =  arc  tan  #z/§y. 

As  noted  in  Chapter  II,  the  time  rate  of  change  of  $  at  launch,  denoted  by 
4>q,  is  called  the  initial  cross-spin  and  is  one  of  the  significant  factors 
contributing  to  rocket  inaccuracy. 

A  third  complex  angle  that  is  of  significance  in  describing  rocket 
motion  is  the  angle  A,  called  the  yaw,  which  is  determined  by  the  geometric 
axis  of  the  rocket  and  the  tangent  to  the  trajectory.  In  component  form, 
it  is  represented  by 

and  its  magnitude  and  direction  are  found  similarly  to  those  of  0  and  # 
above.  It  is  related  to  0  and  #  by  the  equation 

A  =  #  -  ®. 

3«2.  Characteristic  Functions 

Fe  use  the  term  characteristic  function  to  refer  to  that  mathematical 
expression  (or  'formula)  which  determines  at  any  point  during  burning  the 

•  ». 

effect  of  a  unit  amount  of  a  disturbing  factor  (such  as  initial  cross-spin  It ) 
on  one  of  the  significant  quantities  related  to  rocket  accuracy  (such  as  the 
angular  deviation,  0,  discussed  above).  There  is  thus  a  separate  characteristic 
function  for  ®  corresponding  to  each  of  the  factors  causing  rocket  inaccuracy, 


which  were  discussed  in  Chapter  2  of  this  report.  Likewise,  there  is  a  set 
of  such  characteristic  f auctions  for  yaw.  A,  and  for  linear  displacement,  R. 
In  order  to  distinguish  between  the  characteristic  functions  due  to 


the  various  disturbing  factors,  we  use  a  notation  in  which  the  subscript 

indicates  which  factor  applies.  Thus,  for  example,  a  subscript  q(as  in  ©  ) 

Q 

always  refers  to  a  characteristic  function  related  to  initial  cross-spin. 


The  characteristic  function  giving  the  effect  on  ©  of  a  unit  amount  of 
cross-spin  (namely,  one  radian  per  second)  is  denoted  by  0  /<§o.  Thus, 
in  the  example  of  0  given  on  Page  12,  if  the  value  of  ©  given  there  repre¬ 


sents  the  value  of  the  characteristic  function  giving  the  effect  on  angular 
deviation,  9,  of  initial  cross-spin,  $  ,  we  would  write 


8  /§  =  -  0.0070  +  i  O.OO^O 
q'  o 

and,  unless  otherwise  indicated,  the  units  would  be  rad/(rad/sec).  Thus 
multiplying  this  by  1000  would  convert  the  units  (approximately)  to 
mils/(rad/sec).  Again,  as  wre  saw  on  Page  12,  the  magnitude  of  this  complex 
value  is  i  0  /#  |  =  0.0086  rad/ (rad/sec),  or  8.6  mils/(rad/sec).  This 

means  that  an  initial  cross-spin  of  0.5  rad/sec  would  produce  U-3  mils  of 


angular  deviation  at  that  point  of  the  trajectory  for  which  the  characteristic 
function  was  computed.  One  notes  that  the  notation  is  quite  natural, 

since  the  total  amount  of  ®  due  to  a  given  amount  of  #  ,  divided  by  that 

•  q  .  o 

$  gives  the  amount  of  9  due  to  a  unit  of  #  ,  which  is  what  we  mean  by  the 
o  q  - o 

characteristic  function  value. 


It  is  desirable  to  deal  with  characteristic  functions  in  connection 
with  accuracy  computations,  since  values  of  these  functions  depend  only  on 
parameters  such  as  physical  measurements  of  the  round,  launch  velocity 
(angular  and  linear),  acceleration  produced  by  the  thrust,  etc.  Thus,  in 
dealing  with  accuracy  computations  for  a  group  of  rounds  of  a  given  rocket 
type,  with  consistent  launching  parameters,  the  value  of  the  characteristic 
function  depends  essentially  on  the  location  down  range  (during  burning)  for 
which  it  is  computed,  and  hence  does  not  have  an  appreciable  round-by-round 
variation.  On  the  other  hand,  the  value  of  cross-spin  at  launch,  varies 
from  round  to  round  (see  [C])  even  under  the  most  carefully  controlled 
firing  conditions.  Thus  a  complete  analysis  of  rocket  behavior  falls  into 
two  quite  separate  categories,  first, the  determination  of  values  of  the 
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characteristic  functions  for  prescribed  design  values  and  firing  conditions 

(launch  velocity,  acceleration,  etc . ) ,  and  secondly,  the  determination  of 

« 

the  distribution  of  by  experimental  measurements.  The  first  phase  of 
this  analysis  is  of  chief  concern  in  this  report. 

Characteristic  functions  which  appear  later  in  this  chapter  will  be 
designated  as  follows; 


Aq/^orad. / (rad. /sec . ) 
0q/*o  rad. /(rad. /sec.)  > 
R  /#o  ft. /(rad. /sec. )  ^ 


Effects  on  A,  Q,  R  of  initial 
cross-spin,  #Q(rad./sec.). 


0p/Pc (rad. /rad. ) 

R  /p  (ft. /rad.) 

P  c 


Effects  of  dynamic  unbalance, 

S  (rad.). 

'  c 


®lAc  (rad.  /ft. ) 
RLAc(ft./ft.) 


Qg/ (rad.  /rad.  ) 
R^A  Q(ft./rad. )  | 


0/^  rad. /(ft. /sec.  ) 
R w/vrc  ft. /(ft. /sec. ) 


Effects  of  linear  thrust  mis¬ 
alignment,  ( f  t . ) . 


Effects  of  initial  yaw,  A^Grad. ). 


Effects  of  constant  cross-wind, 
w  (ft. /sec.). 

C 


0r/rc(rad./ft.) 


Effect  of  static  unbalance  r  (ft.). 

c 


0aAc  (rad.  /rad. ) 


Effect  of  angular  thrust  misalignment 
a  (rad.). 

v 
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Formulas  for  these  quantities  are  derived  from  differential  equations 
which  will  be  given  in  Section  3 . 1+. 

3.3.  Notation 

In  the  differential  equations  and  formulas  for  characteristic  functions 
which  appear  later  in  this  chapter,  the  following  notation  is  used: 

2 

A  =  axial  moment  of  inertia  -  (slugs  -  ft  ). 

2 

B  =  transverse  moment  of  inertia  (slugs  -  ft  ). 
q  =  A/(2B). 

m  =  mass  of  the  rocket  (slugs). 

k  =  B/m  =  transverse  radius  of  gyration  (ft). 

2 

G  =  acceleration  of  the  rocket  (ft /sec  ). 
v  =  velocity  (ft/sec)  of  rocket  at  time  t  (sec). 
vo=  velocity  (ft /sec)  at  t  =  tQ(sec)(at  launch), 
co  =  axial  spin  rate  (rad/sec), 
n  =  co/v(rad/ft)  =  "spin  ratio", 
p  =  m/qn(ft). 

s  =  v2/2G(ft.),  s.  »  v  2/2G(ft). 

0.  Or) 

G  =  acceleration  (ft/sec  )  due  to  jet  thrust. 

-L  2 

n^G-^  =  spin  angular  acceleration  (rad/sec  )  due  to  jet  thrust. 

r  =  s/p  =  v2/2Gp,  rQ  =  SQ/p  =  VQ2/2Gp. 

K^=  aerodynamic  overturning  moment  coefficient  (see  reference  [KM]). 

Kjj=  aerodynamic  lift  coefficient. 

Kp=  aerodynamic  axial  drag  coefficient. 

Kg=  aerodynamic  damping  moment  coefficient. 

K  =  aerodynamic  spin  decelerating  coefficient, 
p  =  density  of  the  atmosphere  (lb/ft^). 
d  =  diameter  of  the  rocket (ft). 

J.  =  pd^K./m,  i  =  M,  D,  N,  H,  A. 

1  222 

S  =  qnk/i  =  aerodynamic  (gyroscopic)  stability  factor. 

o  =  a-i/sV/z. 

%  =  PJN/d>  °H  =  pdVk2»  °M  =  pVk2* 

°B  ’  Vd<  °A  =  dVk2' 
c " 1  * 2  Vo- 
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=  rr(l  +  cr),  h2  =  tt(1  -  cr). 
kx  =  tr/q  -  k2  =  rr/q  -  ty 

wo  =  V2hlro/Tr  ’  -o  =  ^h2ro/Tr* 

w  =  v'2h^r/rr  ,  w  =  V^h^r/rr  . 
w 

C(w)  =  J  cos(nx2/2)dx  =  Fresnel  integral  (see  references  [JE]  or  [D] 

°  for  tabulation) . 
w 

S(w)  =  /  sin(nx  /2)dx  =  Fresnel  integral. 

Jo 

E(w)  =  C(w)  +  i  S(w). 

rc(x)  =  rr(x)  -  i  r j (x)  =  the  complex  conjugate  of  the  rocket  function 

rc(x)  (See  references  [RNG]  amd  [RC]  for 
tabulation).. 

D  =  1  r~(hlr0)- 

.  2 

g  =  acceleration  (ft/sec  )  due  to  gravity. 


e  =  angle  of  elevation  of  tangent  to  trajectory. 

rj  =  I  arit  =  spin  angle  (after  launch). 

_  o 

3»U.  Differential  Equations  Of  Motion  Of  A  Spin-St abili zed  Rocket 

differential  equations  of  motion  of  a  spin-stabilized  rocket  during  the 
burning  period  outside  the  launcher  are  derived  by  Harrington  in  [H-l]  and 
[H-2]  by  making  use  of  fundamental  principles  of  dynamics.  The  results, 
with  certain  insignificant  terms  omitted,  are  given  here  without  details 
of  the  derivations. 

Basic  equations  defining  the  motion  of  spin-stabilized  rockets  and 
whose  solutions  yield  the  characteristic  functions  defined  in  Section  3*2 
the  following : 


•  •  • 

#  -  (2iqco  -  GtfV/p)#  -  (Cl-v/pHvA  +  w  ) 

tl  M  C  1  - 

*=  -(G1LcA2)eiT)  +  (l-2q)pc((D2  -  ifi>)eiT1,  .  (3-U.l) 


are 
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v#  -  (VA)  -  (CNv/p)(vA  +  wc) 

=  -  g  cos  e  +  ®lace^  +  rc(°°^  -  iaOe^,  (3 * U- 2 ) 

®  =  §  -  4, 

v=Gx  -  CDv2  -  g  sin  e,  (3-U-3) 

oo  =  n^  -  C^vto,  (3*4.4) 

R  =  v®,  (3.4.5) 


where  the  dot  indicates  differentiation  with  respect  to  time,  ' 

••22 

t(e.g.,  $  =  d  #/dt  ).  In  this  system  of  equations  the  quantities  A,  $,© 
have  as  reference  line  the  boreline  OZ  of  the  0-ZYZ  coordinate  system  as 
discussed  in  Section  3.1,  and  the  complex  quantity  R  is  referred  to  the 
OIZ  plane  of  that  system. 

The  differential  equations  of  the  ideal  trajectory  referenced  in 
Section  3.1  are  given  by 

2 

v  =  Gx  -  CDv  -  g  sin  e  , 

v£  =  -g  cos  e.  (3.4.6) 

This  pair  of  equations  is  obtained  by  assuming  that  the  only  forces  acting 
on  the  rocket  are  due  to  the  jet  thrust,  the  aerodynamic  drag,  and  gravity. 
It  is  also  assumed  that  this  trajectory  lies  in  the  vertical  plane  through 
the  boreline  of  the  launcher  and  that  the  rocket  axis  (geometric  axis)  re¬ 
mains  tangent  to  the  trajectory.  In  this  situation  it  is  clear  that  the 
complex  ysfw.  A,  remains  zero,  and  that  the  complex  orientation,  3?,  has 
the  value  41?  =  e-  e  ,  where  e  is  the  angle  of  elevation  of  the  boreline 
(OZ).  •  • 

It  is  now  assumed  that  equations  (3.4.6)  have  been  solved  to  obtain 
v  and  e  as  functions  of  t.  With  v  and  e  thus  determined  as  known  functions 
of  t,  and  with  the  value  co  resulting  from  equation  (3*4*4),  one  could  solve 
equations  (3*4.1)  and  (3.4.2)  for  A  and  arid  from  these  results  obtain 
®  -  9  -  A.  However,  it  is  advantageous  first  to  obtain  new  values  of  ® 
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and  §  by  subtracting  from  each  the  angle  e-  bq  representing’  the  change  in 
direction  of  the  ideal  trajectory.  The  new  #  and  ©  thus  obtained  (for 
convenience,  we  shall  not  change  notation)  are  now  referred  to  the  tangent 
to  the  ideal  trajectory  instead  of  the  boreline.  The  angular  deviation 
thus  becomes  the  deviation  from  the  ideal  trajectory,  as  pointed  out  in 
Section  3*1.  It  should  be  noted  that  the  new  linear  deviation  R  (given 
by  equation  (3.i|.5>))  is  now  measured  in  a  coordinate  plane  perpendicular 
to  the  tangent  to  the  ideal  trajectory,  and  having  its  origin  on  that 
tangent  line;  so  that  R  now  represents  the  complex  distance  from  the  ideal 
trajectory. 

In  addition  to  the  changes  made  in  the  variables  $  and  ©,  we  find  it 
convenient  to  change  the  independent  variable  in  equations  (3.U.1),  (3*U»2), 
and  (3 - U- 5)  to  the  new  variable  r  =  s/p  representing  normalized  distance 
along  the  trajectory.  Thus  we  make  the  change 

4>  =  d#/dt  =  (d#/dr)  (dr/dt)  =  (v/p)#' , 

and  similarly  for  the  change  in  #  ,  etc.  After  making  these  changes  of 
variables,  we  eliminate  4?  from  equation  (3.I1.I)  to  get  a  new  equation  in 
the  quantity  vA  ,  by  use  of  We  also  substitute  the  value  =  ©  +  A 

from  equation  (3.I4..U)  into  equation  (3.1*. 2). 

After  all  of  these  changes,  the  resulting  equations  of  interest  for 
later  work  thus  become 

(vA)«  +  (C^  +  CR  -  2irr)  (v  $  1  -  /S)  (vA  +  wc) 

=  -  j^p2G1Lc/(k2v)j  e1T1“inpPc  (l-2^)  (ve11^)  '  (3-U.7) 

+  (g  cos  e  )  £  (p/v)  (CH  +  CD  -  2iir  -  pOj/v2)  j  , 

v©P  =  v’A  +  C^(vA  +  w  ) 

+  (G1p/v)acein  -  inrc(ve1,n)»  ,  O.ij.,8) 

R’=P0,  (3.U,9) 

where  the  primes  indicate  differentiation  with  respect  to  r. 
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3*5.  Formulas  For  Characteristic  Functions 

Me  now  solve  equations  (3.1*. 7)  and.  (3 . ii. 8 )  under  the  assumption  that 
G,  the  acceleration  of  the  rocket  during  burning,  is  essentially  constant, 
so  that  we  may  use  the  relation 

v=  v^Gpr  .  (3.5.1) 

It  is  also  assumed  that  n^  is  constant  and  equals  to  .  Since  the 

effect  of  gravity  is  reproducible  from  round  to  round,  the  terms 

containing  g  do  not  appreciably  affect  the  dispersion  of  a  group  of 

rounds,  and  hence  these  terms  are  neglected. 

Because  of  the  superposition  property  of  linear  differential  equations, 

one  may  solve  equation  (3. It. 7)  considering  separately  the  effects  of  initial 

conditions  A  and  ,  and  then  the  effects  of  p  ,  w  ,  and  L  respectively. 

The  resulting  value  of  vA  for  each  separate  effect  (and  thence  A  )  may  then 

be  substituted  in  equation  (3.I4..8)  and  the  resulting  equation  solved  for  0 

(w  being  omitted  except  in  solving  for  wind  effect).  Finally,  this  value 
c 

of  0  is  substituted  in  (3.I4..8),  and  the  formula  for  R  is  thus  found  by 
evaluating  the  integral 

R  =  p  ®dr.  (3.5.2) 

Such  work  is  carried  out  in  detail  in  reports  issued  by  this  project  group 
and  referenced  as  [H-l],  tH-2],  [B-2],  [BT-1*],  [BT-5],  [F-l],  and  hence  the 
complete  details  of  the  derivations  will  not  be  repeated  here.  Formulas  for 
characteristic  functions  which  are  the  end  results  given  in  the  referenced 
reports  are  listed  below.  Note  that  a  normalized  form  of  notation  is  intro¬ 
duced  for  convenience  in  relating  the  functions  to  each  other. 


Effect  of  Initial  Cross-Spin  _on  Yaw 


_  .  ih_  (r-r  )  ilu  (r-r  ) 

2^  =  y5G/P  Aq/#o  =  {±/2voyfr)le  0  -e  0  ]. 


Effect  of  Initial  Cross-Spin  on  Angular  Deviation 


s_  = 


(*V»/p  )0q/$o 


=  [  C/(2iro-)j  e  ±^1  0  |e(w)  -  E(wo)j 

e^2*0  Je(w)  -  E(wo)] 

.  (V 


(3.5.3) 
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Tables  of  Fresnel  integrals  are  given  in  references  [JE]  and  [D] . 

One  may  also  write  0  in  terms  of  the  rocket  function  rc(x)  as  follows: 

q 

f  l±u  (r-r  ) 

®q  =  (iC/2Tro)  <vh^  Ire  (h^)  -e  1  0  rcCh-^r)] 


_  _  ih  (r-r  )_ 

-  Vh o  [rc(h0r  )  -e  rc(h0r)] 


iiu  (r-r  )  ih~(r-r  )  I 

+  (l/vfl  le  1  0  -e  2  0  ]>. 


(3.5.3a) 


Tables  of  rocket  functions  are  given  in  references  [ENG]  and  [RC] 
Effect  Of  Initial  Cross -Spin  On  Linear  Deviation 

Rq  =  (  J&fih  RqAQ 

=  r0q  +  jiC/Qiricxh-^]  ^>vb"  e  ^  °^E(w)  -  E(wq)]| 

-  [iC/(UTt7h2)j  e  E(w)  -  E(wq)]|.  (3 

In  terms  of  rocket  functions,  this  is 

-  -  [■—  ih-.  (r-r  )__  1 

Rq  =  r®q  -  (C/l^oV^)  j^rc  (tL^)  <-e  rc(h^r)  J 

r  _  ih_(r-r  ) _  1 

+  (C/l;Tro\/h^)  rc(h2rQ)-e  °rc(h2r)j.  (3 


(3.54) 


(3.54a) 


®6  = 


Effect  Of  Initial  Yaw  A_q  On 


Deviation 


=  [  -iC ]/r^/ (2Tfo)  j  |h2  [  E(w)  -E(wq)] 


-h^  yfetthp  e  '  2  0  £  E(w)  -  E(wj)'j 

r  ih,  (r-r  )  ih  (r-r  ) 
*  (i/VV)[h2e  1  0  -hjS.  2  0 


r 


(3,5,5) 


Note  that  the  0u  used  here  has  the  value  -  1  at  r  =  r  . 

o  0 
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Or,  in  terms  of  rocket  functions, 


j  [-  ih  (r-r  ) _  -i 

®6  =  (C  Vr^/2trCT)  ^  h2  vE^  |rc(h^ro)  -e  0  rc(h^r)J 


r _  ih?(r-r  ) _ 

-  |^rc(h2ro)  -e  0  rc(h2r) 


(i/vE) 


h2e 


ihu  (r-r  )  ih  (r-r  )n 
^  “-l^e  2  ”1 


*  2CnV 


(3-5-5a) 


An  approximation  for  this  expression  is  given  by 


ig  =  -ih1vE0  8  -  c  rc(h1ro)  +  2G^tq.  (3-5-5b) 


Effect  Of  Initial  Yaw  A_  _on  Linear  Deviation 
Rg  =  (l/p)Rg/Ao  ±  -il^  v(T  Rq  +  (CD  -  l)(r-rQ).  (3-5.6) 

Effect  Of  Dynamic  Unbalance  On  Angular  Deviation 

®o  “  ®q/P  -  ±K  Vf  ®  +  CD.  (3-5-7) 

j3  p'  rc  3L  o  q 

Effect  Of  Dypamic  Unbalance  _0n  Linear  Deviation 

R  =  (l/p)R  /p  =  ikL  R  +  CD  (r-r  ).  (3-5-8) 

P  P  c  x  o  q  o 

Effect  Of  Constant  Cross-Wind  _0n  Angular  Deviation 

8  w  =  V^Gp  %/*c  =  1/  v?  +  (1/  v^)  (3-5-5) 

Effect  Of  Constant  Cross-Wind  w  .On  Linear  Deviation 


\  =  v^gTT  R/Wc  =  2(v?  -  vT)  +  (l/vT)^.  (3-5-10) 
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Effect  Of  Linear  Thrust  Misalignment  On  Angular . Deviation , 
®L  =  -(k2/p)  Sj/L,  =  1  ^o(k2r0^/^2  '^2}  ®q  +  CD/(2rQk1k2) 


(3.5.11) 


Effect  Of  Linear  Thrust  Misalignment  1^  On  Linear  Deviation 
=  -(k2 /p2)]^/!^  =  [  i  rc(k2rQ)/(2  v£P  +  CD (r-rQ) 7(2^1^^) 


(3.5.12) 


¥e  have  listed  here  only  one  of  the  characteristic  functions  for  yaw.  A, 

but  others  may  be  found  in  [BT-l],  [BT-l*],  [BT-5]. 

Setae  further  remarks  are  in  order  relative  to  the  manner  in  which  the 

characteristic  functions  listed  here  were  obtained.  To  obtain  formula 

(3.5*3),  one  first  drops  all  the  forcing  terms  (w  ,  L  ,  j3  ,  g,  a  ,  r  ) 

c  c  c  c  c 

appearing  in  equation  (3.I4..9)  and  also  assumes  that  =  C^  =  0.  The 
resulting  differential  equation  is  then  solved  under  the  boundary  conditions 
r  =  rQ,  vA  =  0,  (vi)'  =  p#o>  This  result  for  vA  is  then  substituted  in 
equation  (3.!;. 8),  also  using  (3.5.1),  and  the  resulting  equation  is  inte¬ 
grated  using  the  condition  r=rQ,  0  =  0.  The  resulting  value  of  0,  which 
is  really  ©^,  is  now  multiplied  byihe  constant  C  (appearing  in  the  list 
of  notation)  to  give  formula  (3-5.3).  The  effect  of  putting  C^=  C^=  0 
is  to  neglect  the  aerodynamic  lift  and  damping  moment,  but  it  is  shown 
in  [CH]  and  [B-2]  that  multiplication  of  the  end  result  by  the  constant  C 
then  gives  an  excellent  approximation  for  the  solution  resulting  from 
keeping  and  in  the  differential  equation.  Once  0^  (or  any  0 -function) 
has  been  determined,  it  is  a  routine  matter  to  substitute  the  r  esult  in 
(3.5.2)  to  find  R^(or  any  R-function). 

Formula  (3-5-5)  results  from  a  procedure  similar  to  that  for  finding 
0^,  except  that  the  solution  of  the  modified  vA  differential  equation  is 

carried  out  with  boundary  conditions  r  =  r  ,  vA  =  v  A  ,  (vA)  ’  =  0.  After 

0  00 

the  resulting  Qg/AQ  is  found,  in  order  to  incorporate  the  effects  of 
aerodynamic  lift  and  damping  moment  ,  one  must  multiply  the  result  by  C 
and  add  (C  -  l).  This  is  because  of  the  fact  that  in  the  formula  of 


this  chapter  the  initial  value  (at  r  =  r  )  of  ®g/ Aq  is  maintained  as  -1, 
so  that  the  actual  change  in  angular  deviation  during  burning  as  a  result  of 
a  unit  of  initial  yaw  is  obtained  by  adding  one  to  the  expression  given.'  here 
for  <%/Ao. 

Formula  (3-5-9)  results  from  neglecting  terms  in  L  ,  p  ,  g  in  equation 

c  c 

(3.h-7),  again  putting  =  0^  and  integrating  for  vA  using  b oundary 

values  r  =  r  ,  vA  =  0,  (vA)'  =  0.  Use  of  (3.I4..8)  then  gives  an  expression 
for  ©j/wc  which,  after  having  been  multiplied  by  C,  is  in  terms  of  Aq 
as  shown. 

Formulas  (3.b*7)  and  (3. 3.11)  are  excellent  approximations  for  the 
quantities  in  question  for  points  on  the  trajectory  for  which  (r  -  rQ)>0.5. 

The  details  of  their  derivation  and  the  validity  of  the  formulas  are  com¬ 
pletely  discussed  in  IBT-U],  [BT-3]  and  [B-2]. 

In  Chapter  U  the  use  of  these  formulas  for  making  computations  of  values 
of  characteristic  functions  will  be  illustrated. 
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3.6.  Graphs  Of  Characteristic  Functions 

By  assigning  values  to  the  parameters  occurring  in  the  coefficients  of 
the  system  of  differential  equations  (3.U.7)}  (3»U«8),  and  (3-U*9),  one  may 
obtain  solutions  by  use  of  analog  or  digital  computers.  In  this  section 
analog  computer  graphs  of  some  of  the  characteristic  functions  referenced 
in  Section  3*2  will  be  shown  and  discussed. 

The  parameters  assigned  here  for  the  purpose ' of  computing  are  based 
on  measurements  taken  from  rockets  of  types  that  were  fired  on  experimental 
programs  carried  out  in  connection  with  the  rocket  research  work  of  this 
project  group.  Most  of  the  graphs  shown  here  are  connected  with  gun-boosted 
rockets.  For  the  sake  of  comparison,  a  few  graphs  relating  to  unboosted 
rockets  are  given. 

Graphs  exhibited  here  represent  complex  plots  of  the  quantities  in 

question,  with  the  real  axis  R  taken  vertical  and  the  imaginary  axis  I 

taken  horizontally.  Furthermore,  points  indicated  along  the  graphs  by 

small  circles  are  shown  at  intervals  of  0.2  for  values  r  -  r  ,  the 

o 

normalized  distance  along  the  trajectory  from  the  launch  position  tq.  These 
points  are  labelled  at  integer  values  of  r  -  rQ. 

For  instance.  Figure  3-2  shows  two  graphs  of  the  yawing  motion  due 
to  initial  cross-spin  in  the  case  of  an  unboosted  rocket  for  two  different 
launcher  lengths,  one  of  which  is  four  time  the  length  of  the  other.  Typical 
lengths  here  would  be  3  ft.  and  12  ft.  One  notes  that  there  is  little 
difference:  in  the  magnitudes  of  yaw  (represented  by  the  distances  of  corres¬ 
ponding  points  such  as  r  -  rQ  =  0.2,  0.^,  0.8,  etc.  shown  on  the  graph  from 
the  origin  of  coordinates),  with  the  shorter  launcher  showing  larger  magni¬ 


tudes  during  the  phase  up  to  r  -  r  “  1.  To  illustrate  a  reading  from  this 
graph,  note  that  a  radial  line  OP  is  drawn  from  the  origin  to  the  point 


where  r  -  tq  =  0.6,  and  the  length  of  this  line  measured  to  the  indicated 
scale  of  either  axis  represents  the  magnitude  of  V2G/p  A  /<£  ,  or  the 
value  of  | A  /§ 

_ _q  o  ? 

case  y2G/p  ]  A  /4.  ]  «  0.52.  Thus  if  G  =  21*00  ft. /sec.  and  p  =  120  ft. 

for  this  rocket,  one  would  have  iy*0i  -  o  . 082  rad. / (rad.  /sec )  of  initial 

cross-spin.  Thus  at  p(r  -  rQ)  =  120(0.6)  ==  72  ft.  along  the  trajectory 


i  multiplied  by  the  normalizing  factor  V2G/p.  In  this 


the  yaw  due  to  one  unit  of  is  82  mils.  Also,  the  plane  in  which  this 
yaw  takes  place  is  located  by  the  angle  {~ll5°)  which  the  radial  line  makes 


FIGURE  3.2 
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with  the  "vertical  axis  (measured,  clockwise  from  the  R(+)  axis).  One  could 
also  arrive  at  these  same  results  by  reading  the  real  and  imaginary  com¬ 
ponents  of  the  normalized  characteristic  function  to  get 

V^7p"  ■  A  /#  =  -  0.16  +  i  0.50  . 

q  0 

This  is  the  form  in  which  results  are  given  by  the  characteristic  function 
formulas  in  Section  3-U-  One  would  now  determine  the  magnitude  to  be 

VWF\  Aq/#0  I  =  [(-0.16)2  +  (0.50)2]1//2=  0.52, 

as  before. 

Figure  3*3  shows  a  yaw  graph  for  a  gun-boosted  rocket.  The  general 
behavior  is  much  the  same  as  for  the  unboosted  rocket,  but  magnitudes  are 
considerably  less.  One  should  keep  in  mind,  however,  that  the  value  of 
V2G/p  would  still  need  to  be  reckoned  with  in  measuring  the  actual  yaw. 

One  also  notes  that  the  rate  of  damping  of  the  yaw  is  a  little  slower  in 
Figure  3.3  and  the  precession  rate  is  faster  than  for  the  unboosted  case. 

Figure  3.U  shows  a  family  of  graphs  based  on  varying  launcher  length 
(or  launch  velocity  in  case  length  is  fixed)  and  showing  the  effect  of 
initial  cross-spin  4>q  on  angular  deviation  for  an  uriboosted  rocket.  One 
may  interpret  these  graphs  as  showing  that  as  launcher  length  is  increased 
the  effect  of  a  unit  (rad. /sec. )  on  magnitude  of  angular  deviation  at 
a  given  value  of  (r  -  rQ)(or  distance  along  the  trajectory)  becomes 
smaller.  Otherwise,  one  may  conclude  that  for  a  fixed  launcher  length 
a  higher  launch  velocity  gives  a  relatively  smaller  magnitude  of  angular 
deviation  due  to  unit  cross-spin  at  launch. 

Figure  3*5  shows  a  family  of  graphs  for  the  normalized  linear  devia¬ 
tion  corresponding  to  the  same  set  of  parameters  r^  used  in  Figure  3-U- 
These  two  figures,  3.U  and  3.5,  illustrate  respectively  results  of  com¬ 
putation  one  could  make  by  use  of  formulas  (3-5*3)  and  (3. 5-U)  with  C  =  1. 

Figures  3-6,  3-7>  and  3-8  display  graphs  of  normalised  angular  devia¬ 
tion  due  to  unit  for  cases  of  gun-boosted  rockets  with  three  different 
launch  velocities,  vq (recall  that  v  =  \/2GpF).  These  graphs  are 
practically  identical  in  form,  with  points  corresponding  to  a  given  value 
of  r  -  rQ  almost  in  phase,  but  magnitudes  decrease  with  increasing  launch 


!03  jft  WITH  r0=  7.477 
S  =  1.67  C  =  03469  C  =  07477 
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FIGURE  3.6 
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FIGURE  3.7 


FIGURE  3£ 
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velocity.  From  the  theoretical  formulas  developed  in  connection  with  Section 
3.2,  one  concludes  that  with  other  parameters  fixed  the  value  of  ©  varies 

O  /o  ^ 

approximately  as  (l  +  2C^r  ) /r .  Checking  this  result  against  the  three 
graphs  in  question,  using  r  - r  =  1,  the  ratios  of  successive  values  of 
i  ®q/^0 1  for  r0  =  7-U77,  15,  25  are  2.03,  1.56,  while  for  (l  +  2GNroV*’0'3/2 
the  corresponding  ratios  are  2.12,  1.60.  One  is  thus  led  to  conclude  that 
higher  launch  velocities  lead  to  higher  accuracy,  but  it  is  not  yet  known 
what  the  effect  on  4?o  of  higher  velocities  might  be.  Other  sources  of 
inaccuracy,  such  as  dynamic  unbalance,  might  be  enhanced  by  increased 


velocity.' 

In  connection  with  Figures  3-6,  3-7>  and  3.8,  one  also  notes  that  as 

r  -  rQ  increases  each  of  the  graphs  of  ©^/5q  eventually  winds  around  and 

spirals  inward  toward  a  limiting  position.  Thus,  in  Figure  3*6,  for  instance, 

if  the  rocket  bums  out  at  r  -  r  =  10  (a  distance  of  lOp  =  10(70)  =  700  ft. 

o  • 

fran  launch),  then  the  magnitude  of  Q^/#q  would  be  pretty  well  approximated 
by  the  magnitude  of  the  limiting  value.  Estimates  of  these  limiting  values 
due  to  various  disturbing  factors  are  given  in  the  next  section. 

Figure  3.9  illustrates  one  of  the  significant  results  arising  from  the 


theoretical  and  computational  work  of  this  project  group.  It  shows  two  graphs 
of  103  one  of  which  (the  graph  labeled  0^)  results  as  the  true  graph  for 
the  case  of  a  boosted  rocket,  and  the  other  of  which  (the  dotted  graph 
labeled  results  from  using  the  formulas  applying  to  unboosted  rockets 

to  compute  points  on  the  graph.  In  the  latter  case  the  effects  of  aero¬ 


dynamic  lift  and  damping  are  neglected.  One  notes,  however,  that  the 
essential  difference?  in  the  graphs  is  in  magnitudes.  Note  that  radial 
lines  through  the  origin  almost  pass  through  corresponding  (r  -  r  )  points 
of  the  graphs,  as  illustrated  for  (r  -  r  )  =  1,  and  (r  -  tq)  =  8.  Further¬ 
more,  the  constant  factor  (l  +  2(3^1*  )  (in  this  case  2. 0I4.)  multiplied  by  the 
value  for  dotted  graph  (the  case  of  the  simpler  unboosted  rocket  theory) 
gives  an  excellent  approximation  to  the  value  of  the  solid  graph  (resulting 
from  the  more  complicated  boosted  rocket  theory).  In  £b-2]  similar  re¬ 
sults  were  found  to  hold  for  other  effects  (dynamic  unbalance,  initial  yaw, 
etc.  ),  so  that  a  very  simple  transition  can  be  made  from  theoretical 


formulas  for  unboosted  to  the  formulas  for  boosted  rockets  given  in  Section 
3.5  above. 
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Figures  3.10,  3.11,  3.12,  3.13  show  respectively  the  graphs  for  ®g/A0> 

Qp/(Bc,  v£Gp“©w/wc,  and  -106 (k2/p)  S^/l^  (effects  of  initial  yaw,  Aq,  dynamic 

unbalance,  |3  ,  cross-wind,  w  ,  and  linear  thrust  misalignment,  L  ),  for  the 

case  of  rQ  =  15,  thus  corresponding  to  the  rocket  of  Figure  3*7.  Let  113 

look  at  results  given  by  these  graphs,  assuming  that  we  are  dealing  with 

a  rocket  for  which  the  following  parameter  values  hold:  G  =  300  ft. /sec.  , 

p  =  100  ft.,  k2  =  0.5  ft.2,  r^  =  25(e.g.  burnout  distance  from  launch  is 

(r,  -r  )  =  10,  so  that  s  =  lOp  =  1000  ft.), 
bo 

From  Figure  3*7  we  read  at  (r  -  r  )  =  10  the  value 


103  v5g/p  0q/$o  =  -  2.25  -  i  1.65,  (3-6.1.) 

so  that 

|Q/$o|  =  l.l(lO_3)rad./(rad./sec.  )  =  1.1  mils/ (rad. /sec.)  . 

Thus  =  1  rad. /sec.  produces  1.1  mils  of  angular  deviation  at  burnout. 
From  Figure  3-10  we  read  at  (r  -  rQ)  =  10  the  value 

S^/Ao  =  -  1.028  +  i  0.055,  (3-6.2) 

and  since  the  initial  value  was  -1,  this  gives  a  change  in 

magnitude  0.062  mil/mil,  which  is  the  effect  of  one  mil  of  initial  yaw. 

From  Figure  3*H  we  read  at  (r  -  rQ)  =  10  the  value 


,/pc  =  0.85  -  i  1.25, 

so  that 

|ypc  I  =  1.5  mils/mil. 

From  Figure  3.12,  using  (r  -  rQ)  =  10  we  read 
/  2Gp  Q  =  "0*066  +  i  O.Oll;, 


(3.6.3) 


(3.6.U) 


jQ^/wJ  =  0.00027  rad. /(ft. /sec «.) 
=  0.27  mil/ (ft. /sec.  ), 


so  that 


R  (+) 


.02 


a 

g  WITH  r=l5.0 

S  =1.67  ,  CN=  .03469,  CH=  .0747  7 


0 


I(+) 
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FIGURE  3.12 


FIGURE  3.13 


ill 


From  Figure  3.13  we  read  at  (r  -  r  )  =10  the  value 

-106(k2/p)®LAc  -  U  -  i  1.9,  (3.6.5) 

so  that 

|8lAcI  =  l|.7(l0"k)rad./ft. 

=  O..I4.7  mil/ft. 

From  these  results  it  is  clear  that  for  the  hypothetical  rocket  in 
question  the  factors  having  significant  effect  on  angular  deviation  at  the 
end  of  burning  are  initial  cross-spin,  #  ,  dynamic  unbalance,  Pc,  and 
cross-wind,  w  ,  with  the  latter  the  dominant  factor.  This  is  based  on 
the  assumption  that  initial  cross -spin  would  not  exceed  one  radian  per 
second  and  dynamic  unbalance  would  not  exceed  one  mil.  Then  a  10  mi.  /hr. 
wind  could  cause  about  2*. 5  mils  of  angular  deviation  at  burnout,  as 
against  about  one  mil  for  each  of  the  other  effects. 

Figure  3. Hi  is  included  to  show  the  effect  on  linear  deviation  of 
initial  cross-spin.  Again  at  (r  -  r  )  =  10  it  is  seen  that 

\/2G/p3  Rq/$o  =  -  0.031  +  i  0.0055, 

so  that  for  the  parameters  used  above, 

|  Rq/$0  I  =1*3  ft. /(rad  . /sec.)., 

A  rocket  having  the  qualities  that  would  lead  to  the  results  given 
here  would  thus  be  considered  extremely  accurate. 


v/2GP  Ra 

P2  i 


WITH  C  =  I5.0 


S  =  1.67  ,  C  =  .03469  ,  C  =  .07477 


3-7.  Asymptotic  Estimates  Of  Angular  Deviation 

It  is  noted  from  the  graphs  of  the  characteristic  functions  given  in 
Section  3-6  that  as  r  increases  (or  as  burning  time  increases)  the  curves 
representing  angular  deviation  approach  a  limit  point.  Thus  if  a  rocket 
burns  long  enough,  this  limit  point  will  furnish  a  fair  approximation  to 
the  desired  value  of  the  characteristic  function  at  burnout.  In  this 
section,  under  proper  restrictions,  quite  simple  expressions  giving  good 
approximations  to  these  limit  points  are  listed.  The  details  of  the 
derivation  of  these  formulas  may  be  found  in  reference  [K-$  ]  and  will  not 
be  repeated  here. 

Additional  notation  used  here  involves  the  following  aerodynamic 
parameters  with  typical  magnitudes: 

2  — li  —2 

Overturning  moment :  c^  =  jkA~  h(io“4)(ft.-"). 

Normal  force:  cN  =  JN/d~  3(10”^)  (ft.”1)  . 

Assumptions  made  in  arriving  at  the  results  listed  later  are  that 
during  the  burning  period  outside  the  launcher  we  have 

-  0.08(rad.  1),  |G/v^|  -  Cyt'ft.”1),  v  -  3>00(ft./sec. ) . 

(3.7-1) 

Under  these  conditions  formulas  which  follow  give  approximate  values  for 
the  limit  points  of  the  indicated  characteristic  functions  for  angular 
deviation. 


Effect  Of  Initial  Cross-Spin 

Where  one  is  interested  in  magnitude  only,  the  estimate  is 

®q/$0~  “(?|j  +  G/T02)/cMv0^rad*  per  rad*/sec')*  (3-7.2) 

Effect  Of  Initial  Taw  A 
— - - ■  -  - - o 

p  -  "V 

*  0g/AQ »  -2iqn(cN  +  G/vq  )/c^(radi/rad.  or  mils/mil).  .  ;  .1  ■  (3 •7.3) 

Note:  The  used  here  has  the  value  0  at  r  -  r  ,  and  thus  takes 

t  rrr*  ■  Q  '  0  . 

account  of  the  change  in  ®  due  to  Aq# 


■hk 


Effect  Of  Dynamic  Unbalance 


«  imo  ®^/#o  +  ®g/^o(rad./rad.  or  mils/mil.) 


-  inPn(l-2q)(cN  +  G/v^)/cm. 


(3-7.ii) 


Effect  Of  Constant  Cross-Wind  w 


0  /w  =  (l/v  ) 
f  c  o 


-  (l-v^/v) j  (rad.  per  ft. /sec.).  (3- 7 -5) 


This  is  the  same  formula  as  that  given  in  (3.5*9),  since  this  is  an 

exact  relationship.  It  becomes  an  estimate  if  (3*7*3)  is  used  as  the  value 

of  SL/A  . 
o'  o 


Effect  Of  Linear  Thrust  Misalignment  1^ 

0/L  «  -  (G/k2co  2)(ifl>  0  /$  +  0../A  )  (rad. /ft.).  (3.7.6) 

c  o  o  q  o  00 

Effect  Of  Angular  Thrust  Misalignment  a_ 

*  C 

_  For  n>V lOc^  (e.g.  n>0.1) 


0a/ac  ~  ~i(G/v0ffi0)  (l  ~  ®g/A  o  -  /-va>)  (rad. /rad.  or  mils/mil). 

(3*7*7) 


Effect  Of  Static  Unbalance  r 

- —r, - C 

@r/rc»in(l  -  ®g/0  -e1T])(rad./ft.).  (3,7*8) 

The  latter  two  characteristic  functions  were  not  listed  in  Section  3*5 
because  of  the  fact  that  for  spin-stabilized  rockets  they  are  negligible 
during  the  burning  period.  This  is  indicated  by  the  estimates  given  below 
using  the  formulas  just  listed. 

As  an  example  of  the  use  of  these  formulas  as  estimates  for  values  at 
the  end  of  burning,  consider  a  boosted  spin-stabilized  rocket  with 
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G  =  3!?0  ft. /sec.2,  vq  =  1000  ft. /sec.,  n  =  1.3  rad. /ft.,  and  so<LO°.  We 
note  that  with 

G/vo2*  3(l0~k),  k2«0.3r  •  - . "  '  . . 

the  estimates  of  formulas  (3.7.2)  and  (3- 7-3)  give 

I®  /§  1*1.6  mils  per  rad. /sec.  of  |#  |, 
q  o  o' 

)@_|  wO.l  mil  per  mil  of  |  A  |  . 
o  o 


For  a  comparable  shell,  assuming  no  significant  change  in  (cn_cd^/cm-5 
the  corresponding  estimates  would  be  half  as  large  as  for  the  rocket,  since 
the  term  in  G/v  arising  from  rocket  thrust  would  be  missing. 

Looking  at  the  other  sources  of  dispersion  for  this  rocket,  we  note 
from  (3-7-5)  that  if  at  the  end  of  burning  the  rocket  has  burned  long 
enough  for  the  velocity  to  be  twice  the  launch  velocity  vq,  then  the 
term  is  (l  -  vq/v)  becomes  the  major  contributor  to  ® ^/w^.  Since  the 
term  in  ®6A0  is  almost  pure  imaginary,  as  indicated  by  (3-7-3),  we  would 
then  have 

0 w/wc  =0.001(-  0.5  +  i  0.1), 

or  approximately  0.£  mil  per  ft. /sec.  of  wind. 

For  the  effect  of  dynamic  unbalance,  formula  (3-7-h)  yields  a  value 
on  the  order  of  2  mils  per  mil  of  |  pj  . 

In  the  cases  of  linear  and  angular  thrust  misalignment,  one  notes 
from  (3-7-6)  and  (3-7-7)  that  for  the  rocket  of  this  example 


and 


|  0^/L^  I  «8  mils/ft. 


|©a/oJ  »  0.0002  mil/mil. 


Since  one  would  usually  have  |  L  |  <  0.00  3  ft.  and  ja  |  <  5  mils  for  a  rocket 

.  .  •  •  o  .  c  -  . .  . . 

manufactured  with  reasonable  tolerances,  it  is  clear  that  the  effects  of 
thrust  misalignment  are  negligible  for  this  rocket.  / " 

From  formula  (3.6.8)  the  parameters  assigned  to  this  rocket  yield  the 
value  ’  •  '  ’  \ 


I  ®r/rJ  a*  1.28  rad. /ft. , 


U6 

so  that  if  the  amount  of  static  unbalance  is  limited  to  |  r  |<  0.0005  ft. 

..'C 

the  effect  on  linear  deviation  here  is  less  than  about  0.5  mil. 

It  is  clear  from  this  discussion  that  the  estimates  made  by  use  of 
the  very  simple  formulas  listed  in  this  section  show  that  the  significant 

causes  of  angular  deviation  for  the  rocket  treated  here  are  initial  cross- 

« 

spin,  $  ,  cross-wind,  w  ,  and  dynamic  unbalance,  {3  .  If  one  should  use 
o  c  c 

the  more  complicated  formulas  given  in  Section  3.5,  much  more  accurate 
quantitative  results  would  be  obtained,  but  qualitatively  one  would  arrive 
at  the  same  conclusions  relative  to  significant  sources  of  dispersion  for 
the  subject  rocket  as  were  reached  here. 

3,8.  Application  Of  Theoretical  Results  To  Design  Of  Spin -Stabilized  Rockets 

The  use  of  theoretical  results  listed  in  this  chapter  as  guides  in  de¬ 
signing  an  accurate  rocket  would  usually  begin  at  the  stage  where  the  rocket 
configuration  was  already  prescribed  in  the  sense  that  its  dimensions,  loca¬ 
tion  of  center  of  gravity,  moments  of  inertia,  nose  ogive,  etc.,  would  have 
already  been  determined.  Considerations  such  as  purposes  for  which  the 
rocket  was  to  be  used,  total  weight  desired,  limitations  on  length-td-diameter 
ratio,  range  desired  etc.,  would  have  been  used  in  determining  design 

characteristics  up  to  this  point.  It  would  still  remain  to  determine  such 

/ 

parameters  as  launch  velocity,  vq,  launch  spin  rate,  coq,  and  for  the  burning 
period  outside  the  launcher,  acceleration,  G,  spin-to-velocity  ratio,  n,  and 
burning  time. 

The  following  disturbing  factors,  which  were  discussed  in  Chapter'. 2, 
are  the  prominent  causes  of  inaccuracy  of  spin -stabilized  rockets  in  general: 

1.  Initial  cross-spin, 

2.  Cross-wind,  w  . 

3  c 

3.  Dynamic  unbalance,  j3  . 

c 

U.  Initial  yaw,  Aq. 

5.  Linear  thrust  misalignment,  L  . 

c 

On  the  basis  of  estimates  made  by  use  of  formulas  in  Section  3*7,  or 
by  use  of  the  more  accurate  formulas  given  in  Section  3*5,  one  can  conclude 
that  for  typical  gun-boosted  rockets  the  effects  of  initial  yaw,  and 
linear  thrust  misalignment,  Lq,  are  practically  negligible.  This  las  ; 
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already  been  noted  also  in  connection  with  graphical  results  presented  in 
Section  3*6* 

For-  the  sake  of  more  specific  discussion,  let  us  assume  that  in  the 

remainder  of  this  section  we  have  under  consideration  a  fairly  typical  gun- 

boosted  rocket  with  launch  velocity  fixed  at  v  =  1000  ft. /sec.  and  with 

aerodynamic  constants  C^~]4.(10-4),  ~  3  (ICT4),  and  other  parameters  2 q ~0 . OI4, 

n  ~1.2.  In  order  to  make  use  of  the  estimates  of  Section  3*7,  we  should 

keep  G/vo  -  CM  =  U (10  4),  which  means  that  G  -  I4.OO  ft. /sec.  .  The  value 

G  =  I4OO  ft. /sec. ^  is  still  considerably  higher  than  one  would  want  it  to  be 

in  order  to  attain  best  results  from  the  standpoint  of  reducing  that  angular 

deviation  of  the  rocket  at  burnout.  Equations  (3*7.2)  and  (3*7*3)  show  that 

a  significant  contribution  to  the  values  of  0  /§q  and  is  made  by  the 

quantity  (C^  +  G/v  ^),  and  since  and  remain  nearly  constant  except 

near  sonic  velocity,  decreasing  G  would  decrease  the  effect  of  initial 

cross-spin  on  both  of  these  characteristic  functions  (0/$  and  0k/A  ). 
r  q  o  o  o 

As  a  result,  other  characteristic  functions  which  depend  directly  on  these 
would  be  correspondingly  decreased.  However,  consideration  of  desired  range 
characteristics  might  well  dictate  a  lower  limit  to  which  one  could  go  in 
assigning  G. 

Note  that  in  particular  the  value  G  =  1*00,  along  with  other  parameters 
listed  here,  gives  estimates  of  |0^/#J  ~1*75  mils  per  rad. /sec.  and 
I Q6/Aq [~0.07  mil/mil.  Since  such  an  excessive  initial  yaw  as  10  mils  would 
lead  to  only  0.7  mil  of  angular  deviation,  one  concludes  that  initial  yaw 
is  not  a  significant  factor  in  determining  accuracy  in  this  case. 

The  advantage  of  a  low  value  of  G  is  further  borne  out  in  the  con¬ 
sideration  of  the  effect  of  constant  cross-wind.  From  formulas  (3*7*5) 
the  wind  effect  on  angular  deviation  is  given  by 


V"o'v-  [v,  - (1  -  t<A>]» 

O 

wherein  ®gAQ  may  be  estimated  by  use  of  formula  (3*7*3)*  This  expression 
shows  that  the  nearer  the  velocity  v  at  any  time  during  burning  is  to  the 
initial  velocity  v  ,  the  smaller  the  second  term  involving  (l  -  vq/v)  will 
be.  This  means,  in  other  words,  that  the  nearer  to  zero  the  total  accelera¬ 
tion  G  is,  the  less  the  effect- of  wind  on  angular  deviation  will  be,  and 
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hence  wind  sensitivity  would  be  minimized  by  use  of  a  sustainer  type  of  gun- 
boosted  rocket. 

Considering  the  expression  for  the  effect  of  dynamic  unbalance  given 
in  equation  (3.7»1*),  we  note  again  that  a  low  acceleration  is  desirable. 
Furthermore,  since  the  effect  of  dynamic  unbalance  clearly  increases 
directly  as  the  spin-to-velocity  ratio  n,  and  since  a  rather  high  rate  of 
spin  is  essential  to  maintain  stability  in  flight,  another  recourse  would 
be  to-  increase  launch  velocity  while  maintaining  the  minimum  spin  rate 
required  for  stability.  Remaining  means  of  reducing  inaccuracy  due  to 
dynamic  unbalance  are,  of  course,  the  assignment  of  tolerances  in  manu¬ 
facturing  metal  parts  and  the  choice  of  a  propellant  that  will  not  break 
up  under  the  conditions  of  high  spin  rate. 

Summing  up  the  results  of  the  above  discussion,  we  again  see  that  the 
chief  sources  of  inaccuracy  for  gun-boosted  spinner  rockets  are  initial 

cross-spin,  ,  cross-wind,  w  ,  and  dynamic  unbalance,  f3  .  If  we  combine 
o  c  c 

the  consequent  effects  on  ®^,  the  angular  deviation  at  the  end  of  burning, 
we  get 

®b  '  ®o  +  1  “oPc)(®q/%)  *  (wc/vo)(l  '  Vv)>  (3-8’l) 


where  ®  is  the  direction  of  motion  of  the  c.g.  as  the  rear  end  of  the  rocket 
o  .  & 


clears  the  launcher, 
by  equation  (3. 7. 2). 


and  ®  /#q  is  the  unit  effect  whose  magnitude  is  estimated 
High  launch  velocity  and  the  maintaining  of  a  small 


change  in  velocity  during  burning  tend  to  minimize  the  effects  of  the 

referenced  disturbing  factors.  The  cross-spin  at  launch,  and  the 

dynamic  imbalance,  {3  ,  appearing  in  the  coefficient  of  equation  (3-8.1) 

c 

are  statistical  quantities  which  vary  from  round  to  round  in  firing  a  series 
of  rockets.  For  pc  this  variation  can  be  controlled  by  manufacturing t oleranc e s 
There  is  little  experimental  Evidence  as  t  o  the  behavior  of  for  boosted 
rockets.  Its  magnitude  clearly  depends  on  such  things  as  precessional  mo¬ 


tion  of  the  rear  of  the  round  during  the  tip-off  period  and  launcher  motion. 

From  the  above  remarks,  it  would  appear  that  a  pure  sustainer  type  of 
rocket  (thrust  exactly  cancelling  drag)  would  be  the  proper  choice;  however, 
for  purposes  of  attaining  maximum  range,  it  appears  that  for  certain '  quad¬ 
rant  elevations  a  rocket  with  acceleration  somewhat  higher  than  that  of  a 
pure  sustainer  might  furnish  an. optimum  combination  of  range  and  accuracy.. 
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In  this  connection,  for  high  angle  fire,  it  is  desirable  to  maintain  burning 
almost  to  summit,  for  it  is  desirable  in  all  cases  to  have  as  large  a  velo¬ 
city  at  the  summit  as  possible  under  the  conditions  being  used.  This  follows 
from  the  fact  that  the  magnitude  of  the  yaw  of  repose  is  given  essentially 
by 

2 

Ar  =  2qng  cos  e/C^v  , 

where  e  is  the  angle  of  elevation  of  the  trajectory  (s  =  0  at  summit). 

Thus  small  velocities  vg  at  summit  would  lead  to  larger  relative  varia¬ 
tions  in  A  due  to  small  variations  in  v  than  would  large  values  of  v  . 

2?  S  S 

The  consequent  drift  on  the  downward  path  would  be  less  systematic  for 

small  v  than  for  large  v  and  hence  lead  to  higher  dispersion  of  impact 
s  s 

points. 


CHAPTER  k 


COMPUTATIONS  ILLUSTRATING  ACCURACY  ANALYSIS  FOR  A 
GUN-BOOSTED  SPIN-STABILIZED  ROCKET 

In  the  accuracy  analysis  illustrated  in  this  chapter,  it  is  assumed 
that  the  rocket  has  already  been  designed  and  that  the  following  set  of " 
typical  values  has  been  assigned  to  significant  parameters.  Qaethen  wishes 
to  compute  effects  of  disturbing  factors  by  use  of  the  formulas  listed  in 
Section  3*5-  We  illustrate  such  computations  here, 
li.l.  Rocket  Data 


At  Launch 


v  =  300  ft. /sec. 
o  ' 


At  Burnout 
Vj_i=2000  ft. /sec 


co  =  100  rev. /sec. 

0  ? 
A  =  llj.0  lb. -in. 

B  =  3000  lb. -in. 2 

14;  lb. 

18  in.  from  nose 


0^=390'  rev. /sec 
A  =110  lb. -in. 2 
B  =2300  lb. -in. 
3h  lb. 

17  in. from  nose 


Linear  Velocity: 

Angular  Velocity: 

Axial  moment  of  inertia: 

Transverse  moment  of  inertia: 

Weight : 

Center  of  gravity: 

Other  parameters  are 

Burning  Time:  2.5  sec. 

Diameter  of  rocket:  d  =  U-5  in. 

Aerodynamic  constants  (for  vq  =  500  ft. /sec.):  K^=  2,K^=  l,Kg  =  6 
Other  physical  constants  to  be  used  are 

Air  density:  p  =  0.002335  slug/ft. ^ 

2 

Gravitational  constant:  g  =  32.17  ft. /sec. 

Computation  Of  Basic  Parameters 

We  confute  the  acceleration  of  the  rocket  (assumed  constant)  by 


U.2. 


_  b  o  _  1500  ft. /sec,  _  jh,  /  2 

*  -r — r—  — — 0  - L  600  ft. /sec. 

i^-t  2.5  sec.  ' 


G 


By  referring  to  the  list  of  notation  in  Section  3»3>  we  compute 
the  following  basic  quantities  for  use  in  the  characteristic  function 
formulas,  noting  that  here  launch  values  are  used  for  moments  of  inertia 
weight,  and  for  computing  n.  One  could  use  average  values  here.' 

In  the  following  computations,  free  use  is  made  of  numerical  - 
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results  appearing  in  [R]. 


n  -  »A  -  XS&Jzr'A  .  1.26  rad./it., 

500  ft. /see. 

q  =  A/2B  =  -^0  — p-  =  0.0233, 

2(3000  lb. -in.  ) 

k2  =  B/M  = - 3Qp0  lb. -in_. 2 -  =  0.1*73  ft.2, 

(ill*  lb. )  (liiJU  in.  /ft.  ) 


p  =  tr/qn  = 


TT 


(0.0233) (1.26  rad. /ft.) 


=  107  ft.. 


=  v.  /2Gp  = 


(2000  ft. /sec./ 


2(600  ft. /sec.2) (107  ft.) 
\2 


=  31.2, 


r  = 


(500  ft. /sec.)* 


0  2Gp  2(600  ft. /sec.2) (107  ft.) 


=  1.95, 


pd3/m  -  1°-°Q2?35.  ,  o.soodo)4*, 

(l|li/32.17)  Slug 


JM=  (pd3/**)!^  =  0. 900  (10)' ~k  (2)  =  1.80(10)"^, 

JN  =  (pd3/*)!^  =  0.900(10)^(1)  =  0.900(10)"^, 

Jfi  =’  (pd3/ui)Kjj  =  0.900(10)^(6)  =  540(10)“^. 

i*.3  Computation  Of  Quantities  Occurring  In  Formulas. 

From  the  above  basic  values  the  following  quantities  are  now 

computed.  ?  p  p 

S  =  o2n2k2 /J  =  (°»°233)  (O.U73  ft.*) (1.26  rad./ft.)  _  2  26 

M  ~  1.80  do)"4  '  ’  ’ 


£2 


C  =  PJ  M  -  (107  ffc.)I0.9QQ(lO;4JI  =  0.0257, 

N  n  (U.5A2)ft. 


=  pdJ  /k2  =  (107  ft. )[  (l*.5A2)ft.  ][540(10)~U]  =  0.0I4.58, 


H 


O.U73  ft. 


a  =  (1-1/S)1/2  =  (1-1/2. 26)1/2  =  (0.557522)2  =  0.71*667!*, 


hx  =  tt(1+ct)  =  ir(1.7U67U)  =  2rr(0. 873337)  =  5A873U, 

h2  =  tt  (1  -a)  =  tr(0. 253326)  =  2tt(0. 126633)  =  0.79581*7, 

hxrb  =  (5.1*8731*)  (31.2)  =  171.205, 

h1rQ  =  (5.1*8731*)  (1.95)  =  10.7003, 

h^-r^  =  2n(0. 873337)  (29.25)  =  2tr(25.51i5l), 

h2rb  =  (0-79581i7)(31.2)  =  2l*.630l*, 


h2ro  =  (0. 79581*7)  (1.95)  =  1.&190, 

h0(r -r  )  =  2rr(0. 126663)  (29.25)  =  2*r(3. 701*89). 

2  b  o 

We  new  use  these  values  to  evaluate  the  exponential  functions  and 
rocket  functions  which  appear  in  formulas  (3«5.3&),  (3«5.1*a),  and 
(3.5.5a).  Note  that  we  now  take  r  =  rfe  in  order  to  compute  values 
at  burnout.  For  the  exponential  functions  we  reduce  results  to 
rectangular  form  and  get 

ih_.  (r  -r  )  i2ir(25.51*5l)  i2fr(25)  i2w(0.5i*5l) 

.  1  b  0 _  =  „ 


i2Tr(0. 51*51)  13-U25 


=  cos  3.4*25  +  i  sin  3-14-25, 


=  -  0.960108  -  i  0.279629, 
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i2Tfj 

wherein  use  has  been  made  of  the  property  that  e  =1  for  j  an 
integer,  and  the  Euler  formula 

ix  .  . 

e  =  cos  x  +  x  sin  x. 

In  like  manner, 

ih0(r,  -r  )  i2m(3. 701*89)  i2Ti(0. 7OI4.89) 

„  2  b  o 

e  =  e  =  e 

-  e^*1*290  =  -  0.279611  -  i  O.96OII3. 

To  evaluate  the  rocket  functions  rc (x)  appearing  in  the  characteristic 
function  formulas  referenced  above,  one  may  use  tables  of  these  functions 
given  in  [RC]  or  [ENG],  and  if  values  fall  beyond  the  ranges  of  these 
tables,  the  following  series  forms  are  available  for  use  in  making  the 
computations: 

rc(x)  =  l/Vx  [  (1  -  0.?5/x2  +  6.5625/x^  +  ...) 

-  i  (0.$/x-  1.875/x3  +  30/x5  +  ...).  (U.3.1) 

From  tables  in  [RC}, 

rc(h1rQ)  =  rc (10. 7003)  =  O.30383I  -  i  O.OI3867, 

rc(h2ro)  =  rc  (1.55190)  =  0.701*530  -  i  0.156599. 

By  use  of  the  series  form  given  in  equation  (l*.3.l),  one  finds 

rc(hirb)  =  rc  (171. 200)  =  0.0761*21;  -  i  0.000223, 

rc(h2rb)  =  rc(2l*.83°l*)  =  0.2001*1*1  -  i  0.001*017. 

Remaining  quantities  in  formulas  (3. 5* 3a),  (3-5.1*a),  and  (3. 5* 5a) 

are 

C  =  1  +  2  C^r  =  1  +  2(0,0257) (1.95)  -  1.10, 

iG/2ircr  = - il',-° -  =  i  0.231*1*67, 

2m(0. 71*6671*) 

=V  5.1*873U  *  2.31*251, 

=  0.79581*7  =  0.892103, 

1  /yftT  =  vi/31. 2  =V  0.0320513  =  0.179029. 


■  Substituting  the  above  values  in  (3- 5- 3a)  now  gives 

®q  =  (±0- 23I4.J4.67)  |(2.3U25l)IO. 303831  -  i  0.013867 

-  (-0.960108  -  i  0. 279629)  (0. 076U214.  -  i  0. 000223)] 

-  (0. 892103)10. 70i;530  -  i  0.156399 

-  (-0.279611  -  i  0. 960113)  (0.200441  -  i  0.004017)] 

+  (0.179029)1-0.960108  -  i  0.279629) 

-  (-0.279611  -  i  (0. 960113)  ]| 

=  (i  0.234467) (0.079975  +  i  0.107924), 

and  thus 

0  =V%i7p  ®  /*  =  -  0.025304  +  i  0.018752  (4.3.2) 

Q  Q  O 

The  unnormalized  form  of  (4*3.2)  thus  becomes,  on  using  the  values 

of  p  and  G  given  above, 

•  _ 

9  /$  =  Vp/2G  (-0.0  25304  +  i  0.018752) 

(j  o 

,  \ 

=  (O.2986O8) (-0.025304  +  i  0.018752), 

so  that  we  finally  get  for 

Angular  Deviation  Due  To  Initial  Cross -Spin  : 

Qq/$^  =  -  0.007556  +  i  0.005599  rad. /(rad. /sec. ) 

=  -  7-556  +  i  5.599  mils/ (rad. /sec. ),  (4.3*3) 

and  hence  for  the  magnitude  we  get 

]©q/i0]  =  K-o.007556)2  +  (0.005599)2]1/2 

=  0. 009404  rad«/(rad. /sec.)  (4-3-4) 

=  9.404  mils/ (rad.  /sec . ) . 
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In  a  similar  manner,  substitution  into  formulas  (3. 5 -4a)  and  (3 *5. 5a) 
of  quantities  computed  above  leads  immediately  to  the  following  results: 

Linear  Deviation  Due  To  Initial  Cross -Spin 

Hq  =vWp3  (R  /iQ)  =  -0.7079  +  i  0.5892,  (4.3.5) 

Rq/i>o  =  -22.62  +  i  18.83  ft. /(rad. /sec.  )  ,  (4.3.6) 

|  R  /$Q  |  =  29.43  ft. /(rad. /sec.  )  .  (4-3-7) 

Angular  Deviation  Due  To  Initial  Yaw  _Aq: 

®6  =  ®gA o  =  "  °‘81*96  +  1  °*2437  mil/mil,  (4-3.8) 

I  ©6/  A  ^  |  =  0.8839  mil/mil.  (4.3.9) 

The  remaining  characteristic  functions  whose  formulas  are  listed  in 
Section  3-5  are  expressed  in  terms  of  the  normalized  functions  whose  values 
are  given  here  in  (4-3.2),  (4- 3 - 5) ,  and  (4-3-8),  and  hence  their  com¬ 
putation  is  quite  simple.  Additional  quantities  entering  these  formulas 
are 

=  a/1755"  -  1.396,  y/F  ^y/JUT  =  5.585, 

D  =  1  -  Vh1ro  ScO^r  )  =  1  -  Vl0.70  (0.3038  -  i  O.OI387) 

=  0.006126  +  i  O.04536, 
k±  =  m/q  -  h±  =  (3.1416/0.0223)  -  5-487 
=  129.3, 

k2  =  tf/q  -  h2  =  134.832  -  5.487  =  129.3, 

V^2  =  y/FITTT  -  11.58, 
rc(k0r  )  =  rc (261.4) 

2  o  .  .  . 

=  -  0.0624 99  -  x  0.0001198. 
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Making  use  of  R^  from  (b«  3- 5)  then  gives 

Rg  =  (l/p)(Rg/Ao)  =  -ih^vT  +  (CD-l)  (r-rQ) 

=  -  i(5.U87)(l. 396) (-0.7079  +  i  0.5892) 

+  [1.10(0.006126  +  i  0.0^36)  -  1]  (29.25) 

=  -  2ii.5U  +  i  6.883. 

Thus  the  unnormalized  form  gives  the  following 
Linear  Deviation  Due  To  Initial  Yaw  A^  : 

R6/Aq  =  (l07)(-2ii.5U  +  i  6.883) 

=  -  2626  +  i  736.1;  ft. /radian 


(li.3.10) 


so  that 


=  -  2.626  +  i  0.736  ft. /mil. 


|R6/Aj  =  2.727  ft. /mil. 


(U.3.U) 


Similar  substitution  into  formulas  (3-5-7)  -  (3-5-12)  respectively 
then  gives  the  remaining  results  as  follows. 


Angular  Deviation  Due  To  Dynamic  Dnbalance 

®p/Pc  -  -3.379  -  i  li-519  mils/mil , 

iQp/Pcl  =  5-6I;3  mils/mil. 

Linear  Deviation  Due  To  Dynamic  Unbalance 

Ra/Pn  =  -  11360  -  i  13520  ft. /rad. 

p  c 

=  -  H.36  -  i  13.52  ft. /mil, 
|Rp/pc  |  =  17.66  ft/mil. 


(U.3.13) 

(U-3.1U) 


(J4..3.15) 

(U.3.16) 


Angular  Deviation  Due  To  Constant  Cross-Wind 

@w/wc  =  -  0.001399  +  i  20.000U89  rad/ (ft. /sec.) 

=  -  1.399  +  i  0.1+89  mils/ (ft. /sec. ), 

|e Jw  I  =  0.0015  rad. /(ft  ./sec. ) 

=  1.5  mils (ft. /sec. ) . 

Linear  Deviation  Due  To  Constant  Cross-Wind  w  : 

R Jvc  =  -  2.71+8  +  i  1.1+73  ft. /(ft. /sec.) , 

|R yWc  |  =  3-118  ft. /r(ft  ./sec. ). 

Angular  Deviation  Due  To  Linear  Thrust  Misalignment  1^: 

9/L  =  0.1171+  +  i  0.1521  rad.  /ft .  (mil/lO.^ft . ) , 
J_i  c 

|®TA*  I  =  0.1921  mil/lO“3ft. 

Jj  c 

Linear  Deviation  Due  To  Linear  Thrust  Misalignment  1^: 
RlAc  =  1+5.38  +  i  39.27  ft. /ft., 

|rl/lcJ  =  60.078  ft. /ft. 


(I+.3.17) 

0+.3.18) 

(1+.3.19) 

(1+.3.20) 

(1+.3.21) 


(1+.3.22) 

(1+.3.23) 


CHAPTER  5 


FIN-STABILIZED  ROCKETS  WITH  SLOW  SPIN 


5.1.  Differential  Equations  Of  Motion 

In  this  chapter  a  mathematical  basis  for  study  of  the  motion  during 
burning  of  a  fin-stabilized  rocket  with  slow  spin  will  be  introduced. 
Differential  equations  and  the  resulting  formulas  for  the  characteristic 
functions  expressing  unit  effects  of  the  various  disturbing  factors  des¬ 
cribed  in  Chapter  2  will  be  given. 

The  quantities  used  to  describe  the  motion  of  a  fin-stabilized  rocket 
and  the  coordinate  system  to  which  they  are  referred  are  the  same  as  those 
described  in  Chapter  3  for  the  case  of  spin -stabilized  rockets. 

The  equations  of  motion  to  be  considered  constitute  a  three-dimensional, 
small-yaw  representation  of  the  motion  of  fin-stabilized  rockets.  The  deriva¬ 
tion  of  the  equations  closely  parallels  the  corresponding  derivation  for 
spin-stabilized  rockets  in  [H-l]  and  [CH].  Closely  related  material  con¬ 
cerning  the  equations  of  motion  and  their  derivation  is  to  be  found  in  [MKR] 
and  [DFB]. 

The  following  four  equations  describe  the  motion  of  the  rocket  and  re¬ 
late  respectively  to 

(a)  the  velocity  of  the  rocket  in  its  trajectory, 

(b)  the  spin-rate  about  the  rocket  axis, 

(c)  the  angular  rotation  of  the  rocket  about  an  instantaneous 
transverse  axis  through  the  center  of  gravity, 

(d)  the  motion  of  the  center  of  gravity  at  right  angles  to  OX. 


G  = 


"  V 


g  sins; 


(5.1.1) 


<B  =  n^G^  -  c^  vta;  (5.1.2) 

$  -  (2iqco  -  Cgv)#  +(cMv  +  icT<»)(vA  ,+  w  )=;  (5.1.3) 
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■v#  -  vi  -  vA  -(c^v-ic^)  (vA  +  .w  )=  -g  cos  e+  Fc/m;  (5-l.U) 

in  which  much  of  the  same  notation  as  for  spinner  rockets  appears,  namely 
v  =  rocket  velocity  (ft. /sec.), 

2 

G  =  acceleration  of  the  rocket  (ft. /sec.  ), 
co  =  axial  spin  rate  (rad. /sec.), 

=  complex  orientation, 

A  =  complex  yaw, 

0  =  4?  -A  =  complex  angular  deviation, 

2 

g  =  gravitational  constant  (ft. /sec.  ), 

G^=  acceleration  (ft. /sec.  )  due  to  rocket  thrust  outside  the 
launcher, 

2 

n^G^  =  axial  angular  acceleration  (rad. /sec.  )  such  as  might  he 
provided  by  canted  nozzles, 

e  =  angle  of  elevation  of  tangent  to  trajectory, 

2  2 
B  =  mk  =  transverse  moment  of  inertia  (slugs -ft.  ), 

2q  =  A/B  =  ratio  of  axial  and  transverse  moments  of  inertia  (~  0.02) 

w  =  cross-wind  velocity  (ft. /sec.), 

M  =  resultant  of  cross-torques  due  to  misalignment  (and  perhaps 
c 

unbalance), 

Fc  =  resultant  of  cross-forces  due  to  misalignment  (and  perhaps 
unbalance). 

Aerodynamic  parameters  with  representative  magnitudes  (for  rockets 
with  diameters  of  the  order  of  1*  or  inches)  are  the  following: 

Drag:  c^  =  JD/d~  5(10'*)  (ft.-1). 

Spin-deceleration:  c^  =  dJ^/k  ~  (lO'^Hft.  ). 

Damping  moment:  cn  =  dJ  /k  ~  2 (lO-^) (ft.  ). 

■  n  ii 

Stabilizing  moment:  c^  =  J^/k  ~  (l0-j))  (ft.  ). 

Magnus  torque:  c^,  =  d<J^/k^~  (10~^) (ft.~“). 

Normal  force:  c^  =  J^/d ~  6 (10  ^)(ft.'^). 

Magnus  force:  Cj»  =  Jj,  ~(10'^)  (ft.-1). 
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in  which  d  denotes  the  projectile  diameter  and  k  the  axial  radius  of  gyration. 
It  should  be  noted  that  in  (5.1.3)>  the  term  c^v  A  is  preceded  by  a  +  since 
it  enters  in  a  stabilizing  torque  rather  than  in  a  de -stabilizing  (over¬ 
turning)  torque  as  in  the  spinner  case. 

In  the  development  below,  we  shall  have  use  for  the  following 
synib  ols : 

n  =  co/v  =  ratio  of  spin-rate  to  velocity  (rad. /ft.). 


\  =  2tr/\TcJJ“  =  wave-length  of  yaw  (ft.). 


ft 

r)  =  J  c»  dt  =  spin-angle  (after  launch) . 


The  exponential  e  will  appear  in  the  representations  of  those  forces  and 
torques  which  rotate  with  the  rocket.  Thus 

/ 


2^  in 


-(O^A' ")  e 


(due  to  thrust  misalignment) 


-i  mcA  =  < 


+  c#cv2eiT1’ 


(5.1.5) 


F  /m  =  G-,  a  e1T*. 
c  1  c 


(due  to  fin  misalignment) 

(due  to  angular  thrust  misalignment)  ^  6) 


{ 


1 

{ 

\ 


I 


j 

4 


In  these  formulations,  L  ,  \i  ,  and  a  are  the  complex  parameters,  described  | 

in  Chapter  2,  which  represent  "measures"  of  the  respective  misalignments.  i 

They  incorporate  both  a  magnitude  and  an  initial  orientation .  (i.e. ,  at 
launch,  t  =  t  ).  |l  (corresponds  to  a  distance  (ft.);  |  p  |  and  |a  |  ■ 

O  C  C  C  i 

correspond  to  angles.  Thus  p  represents  that  angle  of  yaw  (measured 

c  j 

relative  to  the. rocket  axis)  at  which  the  cross-torque  due  to  cross-velocity 
reduces  to  zero.  These  parameters  will  be  considered  constant. 

A  torque  due  to  dynamic  unbalance  and  a  cross-force  due  to  static  I 

unbalance  could  be  included  in  (5.1.5)  and  (5.1.6).  However,  these  effects  1 

(assuming  that  reasonable  tolerances  are  maintained)  should  not  be  signifi-  i 

cant  for  fin-stabilized  rockets  with  slow  or  moderate  spin  rates,  where  the . 

spin  is  employed  primarily  to  "average  out”  possible  misalignment  effects. 

The  formulations  up  to  this  point  have  indicated  that  consideration  , 
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is  to  be  given  to  fin -stabilized  rockets  with  slow  spin.  It  should  be  pointed 
out  that  the  results  apply  equally  well  to  finner  rockets  with  no  spin.  To 
relate  the  equations  of  motion  explicitly  to  the  case  of  no  spin,  one  merely 
sets  co  =  0,  n^  =  0,  &  =  0,  n  =  0  in  equations  (0.1.1)  -  (0.1.10,  and  removes 
the  exponentials  e11^  in  (0.1.0)  and  (0.1.6). 

Just  as  was  done  in  the  case  of  spin-stabilized  rockets,  we  now  intro¬ 
duce  a  new  #  and  a  new  ©  obtained  by  subtracting  the  angle  (  s-  eq)  from 
both  the  old  $  and®  ,  where  we  recall  that  (  e  —  e  )  represents  the  change  in 
direction  relative  to  the  launch  direction,  e  ,  of  the  ideal  trajectory 
discussed  in  Section  3-h-  The  new  #  and  ®  thus  represent  the  orientation  and 
angular  deviation  relative  to  the  tangent  to  the  ideal  trajectory.  In  so 
doing,  we  shall  assume  that,  during  the  burning  period  of  the  rocket,  the 
curvature  of  the  trajectory  remains  sufficiently  small  so  that  we  can  ignore 
the  slow  rotation  of  our  new  moving  axis  system  OXIZ  with  OX  tangent  to  the 
trajectory. 

Equations  (0.1.3)  and  (0.1. Ij.)  are  now  rewritten  in  terms  of  the  new 
^•(notation  is  kept  the  same  for  convenience)  and  a  new  dimensionless  inde¬ 
pendent  variable  r  =  s/x  ,  where  s  is  arc  length  along  the  trajectory  (in 
ft.)  and  X  is  the  wave-length  of  yaw  (in  ft.).  Furthermore,  it  can  be 
readily  shown  that  and  c^.  have  a  negligible  effect  on  solutions  of 
these  equations  for  the  cases  of  interest  here  and  hence  they  are  neglected. 
The  resulting  equations  are  then 

+  (v'/v  -  2iqnX  +  XCg)3M  +  (Iir^/v)  (v.A  +  wc) 

•  =  -  iMcX2/Bv2  -  £  g(cos  e)/v2  j  j2iqnX  -  Xcg+v'/V”  Xg(sin  e)/v2j  , 

(0.1.7), 

-  (vA)'/V  -  (XCjj/v)(vA  +  wc)  =  XF^/mv2,  ~  (0.1.8) 

wherein  the  primes  represent  differentiation  with  respect  to  r. 

The  gravity  effects  reflected  in  the  bracketed  terms  at  the  end  of 
equation  (0.1.7)  will  not  be  further  considered  in  this  report.  To  the 
extent  that,  for  a  given  type  of  rocket,  v  ,  to  and  eQ  can  be  reproduced 
from  round  to  round,  the  gravity  effect  upon  the  deviations  of  the  rocket 
will  be  reproducible  and  hence  is  riot  a  significant  source  of  dispersion. 
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If  one  thus  neglects  the  gravity  terms  in  equations  (5.1.7)  and  (5*1.8)  and 
then  eliminates  from  the  resulting  pair  of  equations,  the  following  equa¬ 
tion  results: 

(vA)"  +  ^(c^+Cjj^iqn)  (vA)  '  +  (l|Tr^-2iqnl?cK)  (w  +  vA) 

=  (v'/v  +  2iqnl  -X Cjj)  ( 1  F^/mv)  -iMc X^/Bv-(X/v)  ^  (Fc/m).  (5.1.9) 


Furthermore,  by  making  use  of  the  relation 

<5  =  0  +  A 

in  equation  (5.1.8),  one  may  write  the  equivalent  equation 


&v*= 


(1  +  21  c^r)  A/2r  +  1  c^wc  /v  +  1  Fc/mv 


(5.1.10) 


Once  0  is  determined,  the  linear  deviation  may  be  found  from  the 
equation 

(5.1.H) 


R  /  ©dr 


Jp 

xo 


Equations  (5.1.9),  (5.1.10)  and  (5.1.11)  serve  as  basic  equations  in 
determining  the  characteristic  functions  given  in  the  next  section. 


5.2.  Characteristic  Function  Formulas 

In  solving  the  linear  differential  equations  (5.1.9)  and  (5.1.10),  one 
may  make  use  of  the  superposition  principle  to  consider  separately  the  effects 
of  initial  launch  conditions  (initial  yaw,  Aq,  initial  cross-spin,  initial 
angular  deviation,  ®Q)  and  of  the  misalignments  given  in  equations  (5.1.5) 
and  (5.1.6).  The  unit  effects  of  each  of  these  on  angular  and  linear  devia¬ 
tion  are  listed  in  this  section. 


Effect  Of  Initial  Cross-Spin  $  On  Angular  Deviation 


8  =a/2G/T  0  /* 

Cj  Q  O 


=  (Cy^tr)  e  ^  0  ^(w^)  -  E(w1Q)  j 


,  uiur  r 

+  V2t®i2  e  °  ®(w2)  -  S(w20) 


*  (i/Vr)  ' 


-  in^Or  -  rQ)  -iaigCr  -  r  ) 
e  -e 


] 


(5.2.1) 
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where  the  following  additional  notation  has  been  used  here  or  occurs  later 
in  the  report: 

r  =  s/k ,  where  s  is  arc  length  along  the  trajectory. 

2 

r  =  v  /2G 1  =  launch  value  of  the  dimensionless  variable  r , 
o  o 


“l = 

=  2TT  +  qn  X, 

=  2w  ■ 

-  qn  1 

> 

—  n\  “■  rrL^j 

II 

OJ 

% 

nl  + 

™2’ 

W1 = 

:  V2 my /it  , 

o 

II 

V2my 

/ft  , 

o'  ’ 

W2  = 

:  V2m2r/rr  , 

W20 

y/2.  my 

7^> 

o 

u> 

II 

■  1  +  2\  c^r 

o’  D3 

=  1  - 

vfiyT  rc(m1ro). 

E(w)  =  C(w)  -  i  S(w),  conjugate  of  E(w)  as  used  in  Chapter  3. 


In  terms  of  rocket  functions,  formula  (3-2.1)  may  be  written  as  follows 


iq  -  v5S7T  a/*, 


=  (iCyliir)  ^ 


r—  inLjCr  -  r  ) _  1 

a/ rcGuyO  -  e  °  rc(in^r) 

r-  -iuip(r  -r  )  1 

-  rc  (m^)  -  e  °  rc  (n^r) 

♦  CWr)[  --^(r-ro>' 


-  e 


Effect  Of  Initial  Cross-Spin  _on  Linear  Deviation 


R  =  V2G/X3  H  /# 
q  q'  o 


-iitL.  r  r 

=  r  0^  +  (iC^/Smi^)  V2toil^  e  y  Efa.^)  -  E(w10) 


+  (iCy^TOt^)  e  2  0j^S(w2)  -  E(w2Q)j  , 

or  in  terms  of  rocket  functions 

im,(r-rn) 

re< 


Rq  =  r®  q  “  re(m1ro)  -  e1™1  °  rcCay) 

+  (C^/8rfm2)  £  rcO^r^  -  0  x 


(3.2.1a) 


(3.2.2) 


(3.2.2a) 


6k 


Effect  Of  Initial  Yaw  _Aq  on  Angular  Deviation 


w  rJ. 

®  =@/A’ 

o  o ( 


t 

■**  ITL-|  y  £—>>.11 


■*Vo 


=  (iG3vTAtr)  |m2v?rm^'  e  10  J^E^)  -  E(w1Q)  j 


im  r  f 

V 2rfm0  e  0  E(w2)  -  £(w2Q) 


(i/V?)  [ 


m2e 


inn  (r  -  r  )  -inn  (r  -r  ) 
1  o  .  ^  o 


+  m^e 


] 


or  in  terms  of  r  ocket  functions 


iml(r  -  r0)~, 


r\t  r _  -uiin  U  "  x-  /  ^  -1 

^  =  -(C^vf^A^)  ^m2  "’i  [rc^miro^  *  e  rc  (m^r)J 

_  r  -inigCr  -  rQ) 

+  n^Vn^  rc(m2rQ)  -  e  rcfri^r) 

in^Cr  -  rQ)  -im^r  -  r\ 


(l//r) 


m2e 


+  ime 


m^e 


+  C  . 

3 


Note  that  the  ® c  used  here  has  the  value  0  at  r  s  r  . 

o  o 

An  approximation  for  this  expression  is  given  by 


=  -  innVr*  ®  +  CD. 

5  1  o  q  3  3 


Effect  of  Initial  Yaw  A_q  _on  Linear  Deviation 


_*6  -  Rg  +  C^tr  -ro). 

Effect  Of  Constant  Cross-Wind  _on  Angular  Deviation 

®»  ■  ®A>  =  ®6/tc  *  1/t  -  1/v0- 


(5.2.3) 


(5.2.3a) 


(5.2A) 


(5.2.5) 


(5-2.6) 


1 
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Effect  Of  Constant  Cross-Wind  w  On  Linear  Deviation 


R  =  R  Aw  =  R_/v  +  2  (\/r~  -  -  (  r  -  r  )/v  . 

w  w  c  o'  o  o  °  0  (5  2  7) 

Effect  Of  Linear  Thrust  Misalignment .1^  On  Angular  Deviation 

®L  =  k2®LALc  -  (-i/2TrH2  vT)Sq  -  C3D3/2Tr2H1^rej.  (5.2.8) 
For  boosted  rockets  this  is  almost  equivalent  to 


SL  -  (-1/21*^)®^ 


which  in  turn  implies  that 


®lAc  -  -(iGX/TTk2H2vo)QqAo 


(5.2.8a) 


Effect  Of  Linear  Thrust  Misalignment  On  Linear  Deviation 

^  =  RlALc  =  -(i/2nH2vT)Rq  -  C3D3(r  -  r^An^H,^.  (5-2.9) 


Effect  Of  Fin  Misalignment  p.  'On  Angular  Deviation 


0  =  0  /p 

=  (2Ao)fr /RnVr~  +  12*  VT)8 

£  £  0  OQ 


+  (C3D3/H1H2)(il  +  inX(l-q)/rrH1H2ro). 


(5.2.10) 


Effect  Of  Fin  Misalignment  jx,  to.  Linear  Deviation 

R  =  R  Ap. 
ji  c 


=  (2  A2)  (iA2V^  +  i2wvF  )Rq 


°3D3(r-r.)/lilH2  ][k  +  in l(l -q ) / rrEL^Hgr  Qj  .  (5.2.11) 
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Effect  Of  Angular  Thrust  Misalignment  On  Angular  Deviation 

0  =  ®  /a 
a  a!  c 

=  |(i/2\/fp  1  -  Vt®~P  rc  (TfH2rQ)  j 


+  nX  (l  -  2q)  rc  (trH2ro)/2 
-  iC  D  nX  (1  -  2q)/rr2H1H2ro. 


(5.2.12) 


For  boosted  rockets  this  is  equivalent  to 


0at  2nX(l  -  2q)/H^]©( 


=  (liGnX/H2vo)0q/#o. 


(5.2.12a) 


Effect  Of  Angular  Thrust  Misalignment  a  On  Linear  Deviation 


R  =  SAa 
a  a'  c 


=  |(i/a/rp  1  -V  TJH2ro  rc  ^TfH2ro^  ] 

+  2TfnX(l  -  2q)  rc  (mH^r^/V  rrH^  Rq 
-  iUC  D^nX(l  -  2q)(r  -  t  ) 


(5.2.13) 


5.3.  Remarks  On  Derivation  Of  Characteristic  Functions 

Although  full  details  of  the  derivation  of  the  characteristic  functions 
listed  in  Section  5.2  will  not  be  given  here,  some  remarks  will  be  made  rela¬ 
tive  to  the  procedure  for  obtaining  these  formulas. 

The  results  are  based  on  the  following  two  fundamental  assumptions: 

(l)  The  acceleration  G  of  the  rocket  is  constant  during  the  burning 
period,  so  that  one  may  use  the  relation 

v  =V  S\Gr  .  (5.3.1) 
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(2)  The  spin  angular  velocity  co  is  proportional  to  the  linear  velocity, 
so  that  co  =  nv(n  constant)  and  hence 

n  "t  p  "b 

r)  =  |  codt  =  /  nvdt  =  i jf  nldr  =  n\(r  -  r  ).  (5-3*2) 

to  ^  o  ^o 

To  find  the  effect  of  initial  cross-spin,  ^  ,  given  vin  (5-2.1),  one 
puts  cjg  =  c^  =  w^  =  g  =  Mc  =  Fc  =  0  in  equations  (5-1-9)  and  (5-1-10)  and 
first  solves  equation  (5-1-9)  Tor  vA  using -the  conditions  that 

at  r  =  r  ,  vA  =  0,  (vA) 1  =  \  4>q, 

the  last  condition  being  a  consequence  of  equation  (5-.1-8).  The  expression 
for  A  thus  obtained  is  then  substituted  in  (5-1-10)  and  the  resulting  equa¬ 
tion  solved  for  A  under  the  condition  that 

at  r  =  r  ,  0  =  0. 
o 

The  result  is  finally  multiplied  by  (as  defined  on  p.63)  to  account  for 
CN  and  CH  in  accordnace  with  the  result  given  in  [CH],  p-70,  and  [B-2] 
modified  for  the  fin  rocket  case.  The  argument  for  doing  this  carries 
through  in  exactly  the  same  manner  as  for  spin-stabilized  rockets.  The 
resulting  expression  is  then  formula  (5- 2.1),  which  gives  a  very  good 
approximation  for  0^/4^. 

Formula  (5-2.2)  results  from  substitution  from  (5.2.1)  into  (5.1.11) 
and  applying  integration  by  parts.  This  same  procedure,  of  course,  applies 
for  all  cases  where  linear  deviation  R  is  to  be  found  after  having  found 
the  angular  deviation  0. 

The  procedure  for  finding  the  effect  of  initial  yaw  A  q  as  given 
in  (5-2.3)  is  identical  with  that  for  finding  (5.2.1)  except  that  in 
determining  vA  the  boundary  conditions  that 

at  r  =  r  ,  vA  =  v  A  ,  (vA) 1  =  0 
o*  0  0 

are  used^  with 


v  = 
o 
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The  effect  on  angular  deviation  of  cross-wind  w  in  equations  (5.2.6) 

c 

results  immediately  from  writing  equation  (5.1.9)  as  a  second  order  equation 
in  (vA  +  w  )  and  proceeding  as  in  finding  the  effect  of  A  ,  except  that  the 
boundary  conditions  that 

at  r  =  r  ,  (v  A  +  w  )  =  w  ,  (v  A  +  w  ) 1  =  0 
o*  c  c*  c 


are  used  in  this  case. 

To  find  the  effect  of  linear  thrust  misalignment  L^,  we  take 
c„  =  c  =  w  =  F  =  0  in  equations  (5.1*9)  and  (5.1.10)  and  put 


-iMc/B  =  -(GLcA  )e 


inX  (r  -  r  ) 


The  equation  thus  resulting  from  (5.1.9)  is  then  solved  for  vA  using  the 
boundary  conditions  that 

at  r  =  r  ,  v  A  =  (v  ^  '  =  0. 
o’ 


The  resulting  expression  for  A  (denoted  by  AT  )  may  then  be  expressed  in 
terms  of  A^  (the  yaw  due  to  initial  cross -spin  $q)  just  as  was  done  in 
[BT-5] (pp.6-7).  This  expression  for  A^  is  then  substituted  for  A  in 
equation  (5.1.10)  and  the  same  approximation  made  as  used  in  [CH] (pp. 69-70) 
to  take  account  of  having  omitted  c^  and  c^  above.  The  significant  inte¬ 
gral  is  then  evaluated  and  negligible  resulting  terms  are  dropped  to  give 


formula  (5.2.8)  for  The  fact  that  tends  to  zero  rather  rapidly 

as  rQ  increases  accounts  for  the  approximation  given  in  (5.2.8a),  which 
is  obtained  from  (5.2.8)  by  taking  C  =  0. 

The  formulas  for  ©  At  and  0  /a  result  from  a  procedure  similar  to 
that  for  finding  0^,  except  that  the  appropriate  or  J1  is  chosen  from 

(5.1.5)  or  (5.1.6). 

#v 

It  should  be  noted  that  once  the  value  of  ®  has  been  computed  by  : 

N  N  Q 

use  of  (5.2.1),  the  remaining  ®'s  ( (^, etc.)  can  be  computed  rather 


easily.  Tables  of  Fresnel  integrals  for  computing  are  listed  as 
IJE]  and  [D]  in  the  references  at  the  end  of  this  report. 
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5.ij..  Graphs  Of  Characteristic  Functions  For  Fin-Stabili zed  Rockets  With 

Slow  Spin 

In  Section  5.2  formulas  were  given  representing  particular  solutions  of 
equations  (5-1.9)  and.  (5-1.10)  corresponding  to  particular  boundary  conditions 
and  particular  forcing  terms  representing  various  factors  which  are  sources 
of  dispersion.  If  these  same  equations  are  solved  under  the  same  conditions, 
with  certain  basic  parameters  assigned,  by  use  of  an  analog  computer,  one 
obtains  graphs  of  the  type  exhibited  in  this  section.  Then  if  one  is  making 
a  study  of  a  rocket  whose  parameters  come  fairly  close  to  matching  those 
for  which  the  graphs  were  computed,  a  quick  graphical  estimate  of  angular 
deviations  due  to  the  various  disturbing  factors  can  be  made.  Furthermore, 
a  succession  of  such  graphs  obtained  by  varying  only  one  of  the  parameters 
involved  will  give  a  picture  of  how  this  particular  parameter  affects  the 
angular  deviation  due  to  a  certain  factor.  For  example,  the  effect  of  vary¬ 
ing  launch  velocity  is  reflected  in  some  of  the  sequences  of  graphs  here. 
However,  we  simply  give  enough  graphs  to  indicate  some  orders  of  magnitude 
and  to  show  how  a  more  extensive  program  of  computing  might  be  used  to  draw 
conclusions  about  accuracy. 

We  shall  first  give  descriptions  of  the  various  graphs  shown  in  this 
section  and  make  some  general  remarks  about  them.  Later  in  the  section  we 
shall  discuss  some  quantitative  results  which  can  be  obtained  from  these 
graphs.  Use  of  these  graphs  in  discussing  rocket  design  will  also  be  made 
in  a  later  section. 


Figure  5-1  shows  a  graph  of  the  yaw  A^  due  to  initial  cross-spin 
One  notes  that  the  plane  of  yaw  rotates  slowly  in  the  direction  of 


the  spin,  and  that  the  yaw  oscillations  gradually  damp  out. 

Figure  5-2  shows  graphs  of  angular  deviation  due  to  initial 
cross-spin  for  two  cases  of  constant  spin  (see  [P]  for  a  discussion  of 
this)  and  one  case  where  spin  rate  f)  is  proportional  to  linear  velocity 
v.  Note  that  an  increase  in  the  constant  spin  rate  from  25  rad. /sec.  to 
75  rad. /sec.  causes  little  change  in  the  magnitude  of  angular  deviation, 

but  gives  a  greater  deviation  to  the  right  of  the  desired  direction  of 

•  .  -  ....  „ . .  ‘ 

motion,  if  is  directed  upward.  These  graphs  are  typical  of  unboosted 
fin-stabilized  rockets. 
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Figures  5-3  -  5-7  show  graphs  of  angular  deviation  0  due  to  initial 
cross-spin  for  various  effective  launcher  lengths  rQ  and  for  two  different 
spin-to-velocity  ratios  n  (n  =  fyv) .  These  graphs  are  typical  for  boosted 
fin-stabilized  rockets  with  slow  spin  where  spin  rate  is  proportional  to 
linear  velocity. 

Figure  5.8  shows  angular  deviation  due  to  initial  yaw  Aq  for  two 

cases  of  constant  spin,  fi  =  25  rad. /sec.  and  r)  =  75  rad. /sec.,  and  for 

one  case  where  spin  is  proportional  to  velocity.  Launch  velocities  used 

here  were  in  the  unboosted  range  (r  =  O.l)  and  one  notes  that  angular 

o 

deviations  of  the  order  of  0.3  mil/mil  are  induced.  As  in  Figure  5*2  most 
of  the  deviation  is  upward  from  the  desired  direction  of  motion,  if  A  is 
upward. 

Figure  5-9  shows  the  angular  deviation  0g  due  to  initial  yaw  which 

is  typical  of  a  gun-boosted  fin-stabilized  rocket  (launch  velocity  v  *>*500 

2  ° 
ft. /sec.,  G~670  ft. /sec.)  with  spin  proportional  to  velocity.  Note  that 

in  comparison  with  Figure  5-8  magnitudes  are  much  smaller.  In  fact  the 

maximum  angular  deviation  here  is  of  the  order  of  0.06  mil  per  mil  of  initial 

yaw  and  hence  is  a  negligible  effect. 

Figure  5-10  shows  graphs  of  angular  deviation  6  due  t o  fin  misalign - 

M- 

ment  for  unboosted  spin-stabilized  rockets  having  spin-to-velocity  ratios 
ranging  through  the  values  n  =  0.02,  0.05,  0.1,  0.2.  As  n  increases  magni¬ 
tudes  of  angular  deviation  clearly  decrease,  with  magnitudes  of  limiting 
values  roughly  inversely  proportional  to  spin-ratio. 

Figures  5«H  -  5- lii  show  angular  deviations  9.  due  to  fin.  misalign- 

& 

ment  which  are  typical  of  gun -boosted  fin-stabilized  rockets  with  constant 
spin-to-velocity  ratio  n.  Graphs  in  this  group  include  cases  of  three 
different  effective  launcher  lengths  tq  and  two  different  spin-to-velocity 
ratios. 

Figure  5-l5  shows  a  graph  of  angular  deviation  0^  due  to  linear 

thrust  misalignment  for  a  rocket  with  spin-to-velocity  ratio  n  =  0.2 

2 

and  effective  launcher  length  r  =  1.5.  One  notes. that  with  k  ~  0.5, 

0  cf 

A~  180,  the  maximum  value  of  i  ®[/^c  I  ~  6(l0”'?)mil/ft. ,  and  hence  this 
effect  is  negligible  in  this  ease. 

In  turning  now  to  a  further  discussion  of  implications  that  can  be 
drawn  from  graphs  presented  here,  we  first  point  out  that  there  is  some 
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flexibility  in  the  use  of  these  graphs  because  of  the  formula 

v  =  V2GX  r  .  (5-4-1) 

0  0 

Thus  a  graph  corresponding  to  a  fixed  value  of  r^  can  indicate  results  for 
various  launch  velocities  if  one  simply  assigns  successive  values  to  G,  the 
acceleration  outside  the  launcher  (during  burning). 

For  example,  in  Figure  5*6  the  effect  of  1  rad. /sec.  of  initial 
cross-spin,  3>q,  on  angular  deviation,  0^,  is  shown  for  the  case  where 
rQ  =  1.5  and  rj  =  0.2v(n  =  0.2).  With  X~l80  ft.  (wave-length  of  yaw), 
equation  (5«4-l)  gives 

vq  =V3Eog  .  (5.4.2) 

2 

Thus  an  acceleration  G  =  667  ft. /sec.  corresponds  approximately  to  a  launch 
velocity  vq  =  600  ft. /sec.  A  point  on  the  graph  of  Figure  5-6  at  which 
(r  -  rQ)  =  4  gives  a  reading 

• 

yT 2G/X  0q/#Q  =  0.0095  +  i  0.0005,  (5.4.3) 

which  then  leads  to  the  unit  effect 

0q/$o  =  v480/1334  (0.0095  +  i  0.0005) 

=  0.00350  +  i  0.000184  rad/(rad/sec).  (5.4.4) 


Thus,  at  this  point,  which  is  at  a  distance  s  = X  (r  -  r  )  =  (l80)(4)= 

720  ft.  from  launch,  and  if  burnout  has  not  yet  been  reached,  the  angular 
deviation  due  to  =  1  rad. /sec.  of  initial  cross-spin  has  a  magnitude 

I®  /#  1«  0.00350  rad.  »3*5  mils.  (5.4*5) 

q  0 


In  this  case  an  initial  cross-spin  of  0.75  rad. /sec.,  which  might  be  a 

reasonable  value  of  #  ,  would  lead  to  a  deviation  of  about  2.6  mils  at  720 

o 

ft. 


From  the  behavior  of  the  graph  in  question,  it  is  also  clear  that  the 
magnitude  of  ®^/$o  will  not  vary  by  more  than  10  °/0  over  the  remainder 
of  the  burning  period. 


Figure  5«7  shows  the  same  type  of  behavior  for  the  case  / where  rQ  =  5> 
which  means  that  if  vq  is  still  kept  at  the  value  600  ft./sec*  "the  thrust 
has  been  lowered  to  give  a  G  value  arising  from  (using  equation  (5.U.1) 

600  =  v£g(i8o)(5)  , 

2 

so  that  here  G  =  200  ft./sec*  For  the  same  amount  of  rocket  propellant  as 
before,  ohe  would  then  have  a  longer  burning  rocket,  and  hence  at  burnout 
would  be  close  to  the  limiting  value  of  the  graph  in  Figure  5»7>  which  clearly 
gives  a  magnitude  reading  of 

j  0  /$  1=0.00214-  (5.U.6) 

q  o 

which  when  unnormalized  by  using  G  =  200,  \  =  180  gives 

=  0.0016  rad. /(rad. /sec.  )  =1.6  mils/ (rad. /sec*  )>  (5*U-7) 

and  hence  the  angular  ^Leviation  in  this  case  is  about  one-half  of  that  in 

the  case  corresponding  to  Figure  5.6.  This  would  then  indicate  that  for  a 

fixed  launch  velocity,  a  lower  acceleration  outside  the  launcher  leads  to  less 

angular  deviation  due  to  initial  cross-spin. 

One  may  interpret  the  results  obtained  above  from  Figures  5-6  and  5-7 

in  a  different  way.  Suppose  in  the  case  of  Figure  5*7  we  keep  the  accelera-  . 

2 

tion  at  G  =  667  ft./sec.  instead  of  fixing  the  velocity  v  .  Then  from 
(5.U.1)  it  follows  that  for  r  =5  the  value  v  =  1100  ft./sec.  results. 

The  reading  from  Figure  5-7  at  (r  -  r  )  =  I4  is  then 

V  2G/\  0q/$o  =0.0025  +  i  0.000I4, 

and  hence  the  unnormalized  magnitude  is 

j0 /#  1=0.00092  rad. /(rad. /sec.  ) 

q  o 


=  1  mil/ (rad. /sec. )  (5-U-8) 

Comparing  this  result  with  that  of  (5.1|..i4),  one  sees  that  raising  the  launch 
velocity  from  v  =  600  ft./sec.  to  vq  =  1100  ft./sec.,  keeping  G  fixed,  lowers 

the  angular  deviation  for  this  set  of  parameters  from  3.5  mils  to  1  mil  for 

•-  .  _  _ 

the  effect  of  1  rad. /sec.  of  #  . 

'  0 
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Figures  5.3  and  5.6  (as  well  as  Figures  5«h  aitjd  5.7)  show  the  results 
of  increasing  spin  relative  to  angular  deviation  due  to  initial  cross-spin. 
These  graphs  indicate  that  stepping  up  the  spin-t  o-veiovity  ratio  from 
n  ®  0.1  to  n  *  0.2  has  little  effect  an  magnitudes  of  angular  deviation  ®^. 
From  Figures  5*12,  5*13,  and  5.1i|  W©  note  the  effect  on  angular  deviation 
®  due  to  fin  misalignment  of  increasing  launch  velocity  with  other  para- 

H* 

meters,  including  G,  held  fixed.  At  (r  -  r  )  e  I;,  for  instance,  we  note  the 

°2 

following  readings  based  on  G  B  66?  ft. /sec.  and  X  ■  180  ft. 

For  rQ  **  0.5  (Fig.  5.12),  so  that  vq  «*  3^6  ft. /sec.. 


I®  /u  I  *  0.0260  mil/mil. 

r  C 


For  rQ  °  1.5  (Fig.  5.13),  so  that  vq  **  600  ft. /sec., 

I®  M  I  -  0.0130  mil/mil. 

|i  .  C 

For  rQ  ■  5.0  (Fig.  5.1U),  so  that  vq  =  1100  ft. /sec., 

1®^C|  -0.00725  mil/mil. 

Thus  the  angular  deviation  due  to  fin  misalignment  decreases  with  increased 

launch  velocity,  and  for  the  parameters  used  here  it  is  clearly  negligible 

unless  a  rather  large  misalignment  develops. 

Comparison  of  Figures  5.H  and  5*lU  indicates  that  doubling  the 

spin-t o-veloeity  ratio  (in  the  given  range  of  parameters)  practically 

reduces  the  magnitude  of  angular  deviation  due  to  fin -misalignment  by 

half.  Figure  5*1®  shows  a  similar  result  for  (r  -  r  )>  2  in  the  three 

o 

cases  of  n  0  0.©5,  0.1,  0.2. 


.  FIGURE  5.2  GRAPHS  OF  ©tt  SHOWING  THREE  DIFFERENT  CASES  OF  SPIN 


77 


53  GRAPHS  OF  a  SHOWING  THREE  CASES  OF  SPIN 


'suhhhiSSB 

iSiiiP  mi — 


■■■■■■■■■nan 

wmaam 


i— — wbm—bmbmI 


FIG.  ftio  §rGRAPH  FOR  FIN  MALALIGNMENT  > 


FIGURE  5.11 


R(+) 


FIGURE  5.13 


r.  =  5.oo 
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$.$.  Asymptotic  Estimates  Of  The  Angular  Deviation 

Just  as  was  the  case  for  spin-stabilized,  rockets,  one  may  write  very- 
simple  formulas  for  estimates  of  the  values  of  the  characteristic  functions 
listed  in  Section  $.2  for  large  values  of  r,  and  these  short  formulas  may 
be  used  for  qualitative  accuracy  estimates  in  the  case  of  fin-stabilized 
rockets  with  slow  spin.  Fe  njerely  list  the  estimation  formulas  here,  but 
details  of  the  Mathematical  derivation  of  these  formulas  are  given  in 
Appendix  A  of  this  report. 

Restrictive  Conditions  for  using  these  formulas  for  reasonable  esti¬ 
mates  are 

v  -  $00  ft. /sec., 

|  v/v2\  -  10"3  ft.-1, 

<  2  < 

10  CM  -  n  -  0.2$. 

For  the  formulas  listed  here,  n  is  assumed  constant. 

Effect  Of  Initial  Cross -Spin 

Where  magnitude  is  of  chief  concern, 

y*0  -(oK  +  Vto2)/omV  (5*5.1) 

Effect  Of  Initial  Taw  A 
—  '  -  ' -  — o 

Where  one  is  interested  primarily  in  magnitude , 

V4o  *  ‘  2iqn(°H  *  VTo2)/°H  '  -  2i«“oV*o-  ■  (5-5-2) 

Effect  Of  Constant  Cross-Wind  w 

- — . - * - -  — c 

%/wc  =  -(l/v0)(l  -  vo/v  -  ®g/A0).  .  ($.$.3) 

This  formula  is  exact  in  the  form  given  and  becomes  an  estimation  if  the 
value  of  ®gA0  from  ($.$.2)  is  used. 
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Effect  Of  Fin  Misalignment  p 


2 

®>/pc*  I°M/(n  -cM)]  (icoo® /#Q  +  (^/Aq)  (rad.  /rad.  ’  or  mils /mil) . 


Effect  Of  Linear  Thrust  Misalignment  L 


2  2 


®  l/Lc  «  "(G/k  +  ®fi/0(rad./f‘t*)- 


o  o  q'  o  o'  o 


(5.9.9) 


Effect  Of  Angular  Thrust  Misalignment  a 


0  /a  «(iG/v  co  )[(l-  0C/A  -(v  co  /vaOe1^ ]  (rad. /rad.  or  miis/mil). 

</  o  '  °  °  1  S'  o  o  o'  (5.5.6) 


These  formulas  are  useful  in  making  a  quick  evaluation  of  the  suscepti¬ 
bility  of  a  given  rocket  to  the  disturbing  factors  discussed  in  Chapter  2, 
as  well  as  to  determine  which  of  these  factors  are  significant  in  causing 
inaccuracy.  Furthermore,  they  enable  cue  to  make  qualitative  estimates  of 
the  relative  contribution  to  angular  deviation  made  by  these  significant 
factors. 

5.6.  Use  Of  Theory  In  Rocket  Design 

Except  for  symbols  denoting  certain  parameters  relating  solely  to  fin- 

stabilized  rockets,  (e.g.,  the  fin-misaiignment  parameter,  p  ),  the  notation 

c 

to  be  used  in  this  .section  will  be  the  same  as  that  of  section  3«3-  Among 
the  dispersion-producing  factors  which  can  be  significant  for  finners  are 
several  (4?  ,  A,®  ,  L  ,  a,w)  which  have  already  been  discussed  in  re- 
lation  to  spin-stabilized  rockets.  The  chief  additional  factor  is  that 
of  fin-misalignment  noted  above.  The  relative  significance  of  the  various 
factors  depends  both  upon  the  magnitudes  of  the  parameters  which  'measure1 
the  respective  factors  and  upon  a  variety  of  conditions  pertaining  to  the 
nature  of  the  rocket,  such  as  the  method  of  launching,  the  launch  velocity, 
the  magnitude  and  duration  of  the  acceleration  due  to  rocket'  thrust  after 
launch,  and  the  amount  of  slow-spin,  if  any,  which  is  imparted  to  the 
round.  *  '  .. 

Before  one  can  consider  specifically  the  possible  significance  of  a 
dispersion-producing  factor,  it  should  be  noted  that  certain  information 
concerning  the  rocket  under  design  must  already  have  been  prescribed,  at  - 
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least  as  order-of -magnitude  specifications,  within  reasonably  narrow  margins. 
Such  information  would  include  the  over-all  dimensions  and  weight  of  the 
rocket,  the  probable  location  of  the  c.g.  and  the  magnitudes  of  the  moments 
of  inertia,  and  the  estimated  amount  of  rocket-impulse  needed  for  the  desired 
pay-load  and  range  characteristics.  All  of  these  considerations  relate  to 
the  type  of  weapon  desired,  its  over-all  purpose,  and  the  manner  in  which  it 
is  expected  to  be  used.  It  should  be  noted  that  these  same  considerations 
relate  to  the  method  of  launching.  The  weight,  portability,  and  complexity 
of  the  launcher  are  important  factors  in  the  design  of  the  weapon,  and,  in 
turn,  the  nature  of  the  launcher  cannot  be  ignored  in  the  discussion  of 
sources  of  inaccuracy.  It  can  be  significant  in  determining  feasible  launch 
velocities  and  spin-rates  (for  slow-spin),  and  also  the  interactions  of  the 
rocket  and  launching  process  determine  the  launch  parameters,  $  ,  A  ,  and  ©q. 
The  round-to-round  variations  in  these  parameters  in  either  magnitude  or 
direction,  together  with  variations  in  launch  velocity,  v  ,  can  be  among  the 
chief  contributors  to  dispersion. 

It  has  already  been  pointed  out  that  the  purpose  of  imparting  slow  spin 
to  a  fin-stabilized  rocket  is  to  prevent  the  deviation,  arising  from  some 
source  such  as  thrust -malalignment  or  fin-misalignment,  from  building  up 
steadily  in  one  particular  direction.  The  slow  spin  provides  a  mechanism 
for  averaging  out  at  least  partially  the  deviations  due  to  such  a  malalign¬ 
ment.  Figure  5*10  shows  how  increases  in  spin-rate  reduce  significantly  the 
magnitude  of  angular  deviation  due  to  fin-malalignment.  It  should  be  pointed 
out  that  the  spin-rate  should  not  be  such  that,  over  any  extended  portion  of 
the  trajectory,  the  distance  required  for  one  spin-revolution  ns  approxi¬ 
mately  equal  to  the  "wave-length  of  yaw"  of  the  rocket.  (See  Appendix  A). 
Otherwise,  there  would  be  possible  resonance  effected  due  to  forcing  functions 
rotating  with  a  frequency  near  that  of  the  "natural  yawing  frequency"  of  the 
rocket. 

Slow-spin  does  not  serve  to  reduce  the  effects  of  launch  conditions, 

A  ,  3?  .  or  of  wind.  For  a  reduction  of  these  effects,  the  considerations'-  : 

O  O'  .  ■  .  --  ■  •  —  •  -  ■  ■  .  .  . 

of  the  role  of  r^,  the  effective  launcher  length,  are  much  the  same  as  for 
spinner  rockets.  The  graphs  shown  in  this  chapter  indicate  closely  the 
effect  of  increasing  tq)  and  increases  in  rQ  can  be  accomplished  by 
changing  vq  or:  G  or  both.  Invreases  in  rQ  also  serve  to  reduce  the  effects 
of  misalignments. 
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9-7*  Accuracy  Analysis  for  a  Gun-Boosted  Fin-Stabilized  Rocket 

In  this  section  we  illustrate  an  accuracy  analysis  for  a  gun-boosted, 
fin-stabilized  rocket  in  which  the  rocket  configuration  is  determined,  but 
certain  significant  parameters  (G,  X,  etc.)  are  allowed  to  take  various 
values.  Computations  made  in  this  analysis  are. based  on  the  formulas  for 
asymptotic  estimates  given  in  Section  3-3  above. 

Based  on  the  given  dimensions,  shapes  and  weights  of  the  component 
parts  of  the  rocket,  calculations  of  the  axial  and  transverse  moments  of 
inertia  are  made  for  the  two  cases  of  the  rocket  with  propellant  and  without 
propellant.  The  following  results  are  thus  obtained  in  terms  of  the  nota¬ 
tion  given  at  the  beginning  of  this  chapter: 

2 

With  propellant :  A  =  126  lb.  -  in.  , 

B  =  3710  lb.  -  in.2, 

k2=  B/M  =  0.38  ft.2, 

2q=  A/B  =  0.03U. 

Without  propellant:  A  =  102  lb.  -  in.2, 

B  =  3330  lb.  -  in.2, 
k2=  0.63  ft.2, 

2q=  0.029. 

It  is  clear  from  the  estimation  formulas  given  in  Section  3-3  that 
the  aerodynamic  parameters  c^  and  c^  are  of  particular  significance  in 
determining  values  of  angular  deviation  due  to  various  factors.  We  note 
that 

CN  =  Vd  =  ’ 

where  the  magnitude  of  the  normal  force  c orresponding  to  a  small  yaw  | A |  is 

Pd2V2jA|  . 

Also 

~  =  Pd3^/11^  » 

where  the  magnitude  of  the  restoring  (stabilizing)  moment  is 

pd3^2  |a!  . 
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The  size  and  nature  of  the  fins  will  affect  both  the  normal  force  and 
the  restoring  moment.  One  simple  way  of  relating  the  two  is  to  consider  a 
distance  measured  from  the  center  of  mass  to  the  point  at  which  it 
would  be  necessary  to  apply  the  normal  force  so  as  to  produce  the  restoring 
moment.  This  position  is  not . identical  with  the  center  of  pressure  corres¬ 
ponding  to  the  normal  force  but  is  probably  fairly  close  to  it.  In  terms 
of  one  can  readily  write 

and 

=  (k2/4>c 


N 


M 


The  restoring  moment  leads  directly  to  the  wave  length  of  yaw,  X;  that , 
is,  the  trajectory  distance  in  which  one  yaw-oscillation  is  completed.  Re¬ 
call  that 

X  =  2n//c^=  2Tr/7pd3KM/B  . 


In  the  results  tabulated  below,  wave-lengths  of  yaw  ranging  from  100  ft.  to 
1|00  ft.  have  been  used  and  Table  3-1  shows  corresponding  values  of  c^  and 
K^.  Reliable  values  of  and  for  t  he  rocket  will  probably  be  obtainable 
only  by  aerodynamic  testing.  There  is  no  simple  way  of  estimating  the  re¬ 
storing  moment  unless  the  rocket  is  quite  similar  to  another  rocket  for 
which  the  wave-length  of  yaw  is  known.  See,  for  example,  the  discussion 
in  Chapter  2  of  Exterior  Ballistics  of  Rockets  by  Davis,  Follin,  and 
Blitzer  [DEB]. 

In  order  to  keep  c^  and  c^(or  and  K^)  inter -related,  we  have  adopted 
the  simple  expedient  of  estimating  what^j  might  be.  With  the  center  of  mass 
some  15  or  16  inches  from  the  base  of  the  rocket,  one  can  expect  that,  with 
fins  of  adequate  size  (to  correspond  to  wave-lengths  of  yaw  in  the  range 
mentioned  above),  the  distance  might  be  of  the  order  of  10  or  12  inches. 
For  convenience,  the  value  =  1  foot  has  been  used  in  the  computations 

below.  If  were  reduced  to  0.8  ft.,  then  the  significant  ratio  c^/c^  • 
would  be  increased  by  25  percent.  As  the  fin  area  is  increased  to  reduce 
the  wave-length  of  yaw,  can  be  expected  to  increase  somewhat.  No 
attempt  has  been  made  to  estimate  the  variation  of  with  X,  since  the 
basic  JL.,  as  used,  is  at  best  a  rough  estimate.  \ 


Table  5-1 


c^.  and  for  Various  Values  of  X(ft. ) 


X 

100 

ESI 

200 

250 

300 

350 

ij.00 

(103)cm 

3-95 

1.76 

0.987 

0.632 

0.U39 

0.322 

0. 21^7 

(with  propellant) 

25 

11 

6.2 

li.O 

2.8 

2.0 

(without  propellant) 

2k 

10.5 

5.9 

3.8 

2.6 

1.9 

i.5 

The  quantity  (c^  +  v^/v^  ) /c^.  shown  in  Table  5-  2  below  is  fundamental  in 
the  estimates  of  both  and  <3^.  It  should  be  noted,  for  the  range  of  Vq  used 
the  ratio  is  the  major  contributer  to  the  quantity.  This  ratio,  in  turn 

as  was  indicated  earlier,  is  determined  by  k^/^.  Since  has  been  estimated 
to  be  approximately  one  foot,  with  no  variation  of  with  X taken  into  account 
the  numbers  shown  in  Table  5*2  should  be  interpreted  as  significant  only  in 
the  following  respects.  First,  the  table  shows  an  overall  order  of  magnitude. 
Secondly,  for  a  particular  X,  the  table  shows  how  the  quantity  varies  with  vq. 


Table  $.2 


Values  of  (c^  +  V'  /v  ^)/cjj  ^or  Various  Values  of  X(ft.)  and  (ft. /sec. 

With  v  =  1000  ft. /sec. 
o 


-ko 

-20 

0 

20  , 

•  ko 

100 

100 

o.57 

0.575 

o.58 

0.585  j 

0.59 

0.605 

150 

.56 

.57 

.58 

.59 

.60 

.61+  - 

200 

•5U 

.56 

.58 

.60 

.62 

X 

.68 

250 

•  52 

.  .55 

. .  .58  .  . 

-61. 

,6k 

.7k  . 

300 

•k9 

.53 

:  .58 

.63 

.67 

.81 

350 

.1+6 

.52 

.58 

.6k 

.70 

.89  - 

koo 

CM 

* 

.50 

.58 

.66 

-7k 

.98  ; 
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If4  were  taken  to  be  0.8  ft.  instead  of  1.0,  the  entries  in  Table  5-2 
would  be  increased  by  approximately  -25  percent  throughout  most  of  the  table 
but  with  smaller  percentage  increases  when  both  land  are  vlarge . 

It  should  be  noted  that 

„  2 

v  ~  G_  c_v  -  g  sin  e  , 
o  1  D  o  &  o 

2 

where  denotes  the  acceleration  due  to  rocket  thrust,  -  the  drag 

deceleration,  and  -g  sin  the  gravity  deceleration  at  launch.  If  the 

rocket  thrust  were  increased  so  as  to  increase  v  significantly  beyond  the 

2  0  2 

100  ft. /sec.  shown  in  the  table,  then  the  term,  v  /cjjV  ,  would  become 

more  significant  and  would  be  the  dominant  term  at  least  for  the  larger 

values  of  X. 

We  now  turn  to  estimates  of  unit  effects  on  angular  deviation  for  the 
various  disturbing  factors.  For  order  of  magnitude  estimates  of  ©  /§q  and 
V4c  ,  we  use  equations  (5*5*1)  and  (5.5*2)  of  Section  5,  namely 


2 

0q/#o  ^CN  +  VVo  VcMV0^rad*  P6r  rad-/sec.  ); 


(5*5*1) 


2 

36//Ao  ““2i<5n^cN  +  VVo  )/cM(mils  per  mil). 


(5*5*2) 


Effect  of  Initial  Cross-Spin 


With  v  =  1000  ft. /sec.,  estimates  of  Q  /i>  in  milliradians  of  deviation 
o  .  q'  o 

per  rad. /sec.  of  can  be  read  directly  from  Table  5*2.  For  a  sustainer 
type  rocket  with  ]v|~  20  ft. /sec. ^  or  less  during  burning,  one  finds  that 
the  effect  of  initial  cross-spin  is 

0^/$o~  0.6  mil  per  rad. /sec., 

with  the  possibility  that  this  estimate  might  increase  to  0.75  if ^  should 
turn  out  to  be  of  the  order  of  0.8  ft.  rather  than  1.0. 

If  v  is  increased,  the  unit  effect  is  increased. 

The  parameter  depends  upon  the  rocket-launcher  combination.  Experi¬ 
mental  data  on  for  boosted  rockets  are  very  scarce.  Extrapolating  from 
data  for  unboosted  rockets,  one  might  expect  |§q  | '  s  up  to  0.8  rad. /sec., 
perhaps,  and  it  would  be  hoped  that  the  magnitude  could  be  kept  under  1.0 
rad. /sec.  Thus  the  maximum  angular  deviation  that  one  might  expect  due  to 

$  should  not  exceed  0.75  mil  for  this  rocket, 
o 
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Effect  of  Initial  Yaw 


With  2q  =  0.03lj.  just  outgide  the  launcher, 

r 

0.003it  if  n  =  0.1  rad. /ft. 

0.0068  if  n  =  0.2  rad. /ft. 

Thus,  for  the  sustainer  type  rocket,  the  effect  of  initial  yaw  is 


2qn  = 


0  /A 
67  o 


-  0.002i  mil/mil  if  n  =  0.1, 
j-  O.OOl^i  mil/mil  if  n  =  0.2. 


Even  with  a  large  initial  yaw  of  6°(  ~ 100  mils),  the  deviation  0g  would  be 
quite  small. 

In  case  there  were  no  spin,  equation  (A. 36)  of  Appendix  A,  with  c^ 
expected  to  be~  5(10  or  less,  shows  that  |@g/A  |  <10 (mil/mil)  if 

1  Vvo2 1  s  10"1*- 


Effect  of  Fin  Misalignment 


Referring  to  equation  (5-5.il)  of  Section  5,  neglecting  ©g/A0  and 
expressing  0^  as  above,  we  have 

V^c  -  cM)][(cN  +  VVo2)/°M]  ^-5-U 

2 

If  the  factor  nc^/(n  -  c^)  is  calculated  (see  Table  5-3  below),  the  second 
factor  is  available  in  Table  5-2. 


Table  5-3 


Values  of  nc,,/(n  -  c.,)  for  Various  Values  of  \  and  n 
M  M 


\ 

n  =  0.1  ■ 

n  =  0.2 

100 

0.65 

0.22 

150 

0.22 

0.093 

200 

0.11 

0.051 

2^0 

0.067 

0.032 

300 

0.0li6 

0.022 

350 

0.033 

0.016 

1*00 

0.026 

0.013 
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If  one  chooses  0.6  as  the  order  of  magnitude  of  (c^  +  v  /v  )/c  M,  then 

the  corresponding  estimates  of  ©u/p-  would  vary  from  0.39i  to  0.0078i.  The 

c 

basic  reason  for  the  wide  variation  is  that  the  torque  resulting  from  fin- 
misalignment  is  directly  proportional  to  c^(see  equation  (5.1.5))  while  c^ 
varies  inversely  as  the  square  of  X.  This  is  one  effect  for  which  the 
deviation  increases  with  increase  in  stability  (i.e.,  as  c^  increases  and 
X  decreases).  Consider  X=  250  ft.  as  a  possible  wave-length  of  yaw,  so 
that  by  use  of  Tables  5.2,  5-3,  and  formula  (5-5-U)  we  have 


f  0.0ti(mil/mil)  if  n  =  0.1, 
0.02i(mil/mil)  if  n  =  0.2. 


It  is  difficult  to  indicate  what  magnitude  of  |p.  |  one  might  expect. 

c 

On  the  other  hand,  one  might  hope  that  with  feasible  manufacturing  tolerances 

an  upper  bound  of  one  degree  (or  17  mils)  on  |  p  |  might  be  maintained. 

c 

It  should  be  noted  that  with  no  spin  the  deviation  due  to  fin -misalignment 
increases  without  bound  as  trajectory  distance  increases.  See  equation  (A.61i) 
of  Appendix  A. 


Effect  of  Constant  Cross-Wind 


From  equation  (5.53),  the  effect  of  constant  cross-wind  is 

®  =  (w  /v  )  (-1  +  ®~/A  +  v  /v), 

T.T  r>'  n'  O  O  o  ’ 


w 


c'  o 


where  wc  denotes  the  cross-wind  velocity  vector  in  the  plane  normal  to  the 

trajectory.  Thus  the  effect  correspondsto  a  combination  of  the  equivalent 

of  an  "initial  yaw”  deviation  with  A  =  w  /v  and  an  additional  deviation 

o  c  o 


w  /v. 
c 


We  can  write 


®w/wc  =  (l/vo)I(vo/v  -  1)  +  fig]. 

Note  that  under  the  assumption  of  slew  spin  with  n  ~0.1  or  greater,  ®g/AQ  is 

predominantly  pure  imaginary,  whereas  the  remainder  of  the  expression  is  real. 

» 

Futhermore,  considering  •  . 

©g/Ao  0.  bo2;i  (rail/mil)  for  nQ  =  0.2, 

■  one  notes  that  if  there  is  any  appreciable  change:  in  velocity  then  the  magnitude 


of  ® w/wc  is  determined  by 

®w/wc  »(l/v  -  1/vq) (rad.  per  ft. /sec.) 

In  the  case  of  a  trajectory  in  which  v  decreases  from  1000  ft. /sec.  to 

900  during  the  first  U4..I4.  sec.,  the  angular  deviation  due  to  w  will  be 

c 

downwind  and  of  the  order  of  0,1  mil  per  ft. /sec.  of  |wj  .  If  the  velocity 

c 

then  increases  from  900  up  to  99 0  during  the  next  lh.b  sec.,  the  angular 
deviation  will  be  approximately  the  same  in  magnitude  but  upwind,  and  the 
® wArc  over  the  total  burning  period  of  the  rocket  will  have  a  real  component 
which  essentially  reduces  to  zero  and  a  small  imaginary  component  -  O.OOlpL 
mil  per  ft. /sec. 

After  the  end  of  burning,  the  deviation  due  to  cross-wind  follows  the 
same  general  pattern.  If  v  decreases  from  990  to  83O  in  another  lli.ii  sec., 
the  deviation  will  again  be  downwind  and  0.2  mil  per  ft. /sec.  Then  if  v 
increases  from  830  up  to  99 0  again  in  30  sec.,  there  will  be  a  corresponding 
upwind  deviation.  The  net  angular  deviation  of  the  trajectory  due  to  a 
constant  component  of  wind  velocity  normal  to  the  trajectory  should  be  small, 
in  this  situation.  The  actual  displacements  (or  deflections)  of  the  rocket 
over  the  various  subintervals  would  probably  not  cancel  out  to  the  same 
extent  as  the  angular  deviations,  since  the  corresponding  displacement  depends 
upon  the  algebraic  difference  between  the  actual  time  interval  and  the  time 
interval  which  would  be  required  if  the  velocity  at  the  start  of  the  inter¬ 
val  had  been  maintained.  Thus,  over  the  final  interval  (30  sec.)  of  velocity 
increase,  the  upwind  displacement  (or  deflection)  of  the  rocket  would  exceed 
in  magnitude  the  downwind  displacement  during  the  preceding  interval  (II4..I4. 
sec . )  of  velocity  decrease . 

It  is  to  be  noted  that  the  over-all  behavior  of  © 'w/wc  'will  depend  upon 
the  angle  of  elevation  and  the  resulting  velocity  variation  over  the  tra¬ 
jectory.  However,  the  sustainer-type  thrust,  by.  means  of  which  the  velocity 
variation  is  held  within  a  relatively  narrow  range,  is  very  effective  in 
reducing  the  wind  sensitivity. 


Effect  of  Linear  Thrust  Misalignment 
From  equation  (5.5.5),  we  have,  on  neglecting  Sl/A  , 


‘00  )®  /$  =--  (i 

o  cr  o 


CO  V 

0.  o 


+  VTo2)/cM] 
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•under  the  assumption  of  slow  spin.  For  boosted  rockets,  this  effect  tends1' 

to  be  negligible  provided  G^(and  thus  v  )  is  relatively  small  and  provided 

reasonable  tolerances  are  maintained  in  the  rocket  construction  so  as  not  to 

allow  |L  |  to  be  unduly  large.  For  a  k-5  in.  rocket  one  would  expect  that 
c 

|L  I  would  be  less  than  0.01  ft. 

1  c  1 

If,  from  Table  5.2,  we  choose  0.6  as  the  order  of  magnitude  of 
2 

(c„T  +  v  /v  )  /c*,,  then  we  note  that 
N  o'  o  '  W 


O.O17G(l0"3)(ft."2)  if  n  =  0.1, 
0.0085G(l0~3)(ft."2)  if  n  =  0.2, 


and  thus 


®T  A 

y  c 


0.01G(mil  per  foot)  if  n  =  0.1, 
0.005>G(mil  per  foot)  if  n  =  0.2. 


2 

A  U5  lb.  thrust  would  correspond  to  approximately  32  ft. /sec.  of  rocket 
acceleration  Gl>  and  for  this  G, , 


©  j/L__  ~  O.32  mil  per  foot  of 


©  /L  ~  O.OO3  mil  per  10  2  ft.  of  L  . 

L  c  c 

This  would  be  negligible,  and  G^  could  be  increased  considerably  before  thrust 
misalignment  would  become  a  serious  source  of  dispersion. 


CHAPTER  6 


EXPERIMENTAL  RESULTS 

« 

Various  experimental  methods  have  been  used  in  connection  with  J^heJjfork 
of  this  project  group  to  obtain  data  from  actual  firings  of  rocketa^  These  „ 
data  were  in  turn  used  to  compare  theoretical  with  experimental. .^g^lts  and 
also  as  a  guide  to  further  developments  in  the  theory  of  J^qgjcet  motion. 

In  this  chapter  the  outstanding  experimental  devices 'are  described 
and  theoretical  results  arising  from  some  of  the  resulting  data  are  given. 

Some  of  the  devices  used  were  the  following:  ^ 

(1)  The  segmented  rail  launcher  for  study  of  action  of  the  rocket  in¬ 
side  the  launcher.  .  „  ,  v 

% 

(2)  An  optical  lever  device  for  determination  of  the  orientation  of 

the  round  during  the  tip-off  period  and  for  several  feet  of  travel 

* 

beyond  the  launcher. 

(3)  A  Fastax  camera  for  determining  linear  and  angular  velocity  during 
the  first  few  feet  of  travel  after  leaving  the  launcher. 

(U)  One  or  more  pairs  of  cameras  positioned  in  such  a  manner  as  to 
give  trajectory  coordinates  during  burning. 

(5>)  An  impact  field  survey  for  obtaining  impact  coordinates. 

The  facilities  of  the  Transonic  Free  Flight  Range  at  Ballistic  Research 
Laboratories  were  also  used  for  experimental  firings,  and  the  layout  and 
instrumentation  there  is  fully  described  in  [Ro]. 

6.1.  Segmented  Rail  Launchers 

The  11-rail,  9.$  ft.  launcher  (ID  1*.600")  is  essentially .a  smoothbore 
launcher.  It  consists  of  11  contoured  rails  insulated  from  each  other  and 
from  the  outside  launcher  tube,  and  each  rail  consists  of  8  sections  of 
lengths,  respectively,  3n>  3n>  18”,  18",  18",  18",  18"  and  15>",  given  from 
breech  to  muzzle.  The  rail  sections  are  connected  to  a  bank  of  neon  lamps 
which  are  photographed  by  three  3^MM  stripped  film  records.  A  wire  through 
the  front  of  the  launcher  to  the  rocket  is  used  to  couplet e  the  electrical 
circuitry.  As  the  rocket  (banded  to  a  specified  O.D.  at  two  appropriate 
positions  about  20  inches  apart)  moves  along  the  launcher,  contacts  of 
the  rear  and  front  bands  with  the  various  rail  sections  are  thu,s  recorded 

100  .  .  .  'V 
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as  a  function  of  time.  Mathematically  the  data  that  result  furnish  cylindri¬ 
cal  coordinates  of  points  of  contact  in  terms  of  time  for  the  rear  and  front 
bands  of  the  rocket.  Sample  records  from  this  “neon"  launcher  are  shown 
in  Figures  6.  5>,  6. 6  of  this  report. 

A  similar  smoothbore  launcher  3  ft.  long  was  used  earlier  in  the 
experimental  work  described  in  [CHM],  Later  a  ten-foot,  four-rail  launcher 
of  somewhat  similar  construction  with  the  rails  positioned  at  1:30,  iiOO, 

7:30,  and  10:30  o'clock  was  used  for  study  of  motion  of  a  spin -stabilized 
rocket  inside  a  rail  launcher.  Experimental  results  of  a  firing  program 
using  this  launcher  are  given  in  [Ca-2]. 

6.2.  The  Optical  Lever  Device 

This  instrumentation  involved  modification  of  standard  rounds  by  replacing 
the  rocket  fuses  by  fuses  equipped  with  two-inch  front  surfaced  mirrors. 

A  high-speed  camera  (lj.0  in.  focal  length)  focused  on  the  muzzle  of  the 
launcher  was  positioned  80  ft.  down  range  from  the  launcher  muzzle  and  a 
few  feet  to  one  side  of  the  line  of  fire,  and  a  coded  screen  was  positioned 
at  the  same  distance  but  on  the  other  side  of  the  line  of  fire  as  shown  in 
Figure  6.1. __ The  camera  photographs  a  spot  on  the  coded  screen  as  reflected 
by  the  mirror  on  the  nose  of  the  rocket,  and  the  angular  orientation  of  the 
rocket  as  a  function  of  time  for  the  interval  during  and  immediately  follow¬ 
ing  tip-off  can  be  determined  from  the  camera  data.  Reduction  of  these 
optical  lever  data  leads  to  graphs  of  the  components  3?_  of  orientation 

x  U 

4?  as  shown  in  Figure  6.2.  In  turn,  these  graphs  are  used  to  determine 

•  v 

experimental  values  of  initial  cross-spin  #q,  and  a  summary  of  such  values 
arising  from  four  different  firing  programs  reported  on  as  indicated  is 
given  in  Table  6.1. 

At  the  time  of  this  writing,  no  firm  data  yielding  direct  experi¬ 
mentally  measured  values  of  for  boosted  rockets  are  available.  There 
is  some  indirect  evidence  (see  [W-l])  which  indicates  that  for  such  rockets 
3?o's  of  rather  large  magnitude,  possibly  somewhat  systematic,  may  develop 
for  a  particular  rocket-launcher  combination. 

6.3.  Schematic  Layout  ..0£.  Cameras 

Determination  of  linear  and  angular  velocity  at  launch,  together  with 
the  average  acceleration  during  the  first  15>  feet  of  travel  outside  the 


FIGURE  6.2 
6  O'CLOCK  TILT  (72') 
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launcher,  is  determined  from  data  from  a  Fas tax  camera  placed  about  1+0  ft. 
to  the  left  of  the  line  of  fire.  Pairs  of  cameras,  one  on  the  line  of  fire 
and  one  at  a  considerable  distance  along  to  one  side,,  are'  used  to  determine 
position  of  the  rocket  as  a  function  of  time  during  the  burning  period  and 
for  several  hundred  feet  beyond.  The  relative  positioning  of  these  cameras 
is  shown  in  Figure  6.3  below;  some  of  them  may  also  be  seen  in  Figure  6.1. 
Sample  figures  from  the  resulting  data  are  to  be  found  in  1C],  I CD],  and. 
IW-2]. 

Typical  results  for  end-of-burning  parameters  arising  from  reduction 
of  camera  data  and  impact  spotting  are  given  in  Table  6.2  (see  [¥- 2]  for 
further  details  of  results  of  this  firing).  A  summary  of  information 
obtained  from  side  Fastax  data  for  four  different  firing  programs  is  given 
in  Table  6.3. 

Key  To  Notation  In  Tables  6.1,  6.2,  6.3 


Y0 


component  of  initial  cross-spin  in  the  vertical  direction. 


ZO 


*f  = 
*b  = 


vb  = 

co  = 

o 

G  = 

a  = 

n  = 


component  of  in  horizontal  direction. 

time  from  ignition  to  end  of  burning. 

t^  -  time  from  ignition  to  1. 8  in.  of  motion. 

distance  along  boreline  from  launcher  muzzle  to  position  of 
rocket  c.g.  at  bumou:. 
linear  velocity  at  launch. 

linear  velocity  at  burnout. 

spin  angular  velocity  at  launch. 

acceleration  of  rocket, 
spin  angular  acceleration. 

ratio  of  spin  angular  velocity  to  linear  Velocity, 
vertical  component  of  angular  deviation  0  at  burnout. 

horizontal  component  of  0  at  burnout.' 

vertical  distance  from  boreline  at  burnout*  -  ;  , 

horizontal  distance  from  boreline  at  burnout. 


Zh 

|  rJ  =  (Y^  +  distance  from  bore-line  at  burnout.  . 


r.  = 


arc  tan(Z, /£.)  =  direction  in  which  [  R,  |  is  measured  relative  to 
.  '  .  D  vertical.  D  ^  :  ;  i 


FIGURE  6.3 
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6.h.  Rocket  Motion  Inside  A  Smoothbore  Launcher 

In  this  section  results  of  data  obtained  from  neon  records  and  optical  lever 
records  lead  to  a  mathematical  model  for  determining  rocket  motion  inside 
a  smoothbore  launcher. 

(l)  Some  Experimental  Data 

Attention  is  to  be  given  first  to  certain  experimental  data  concern¬ 
ing  the  motion  of  some  M3 3  rockets  which  were  deliberately  misaligned  by 
a  1°12'  tilt  of  the  nozzle-plate.  The  tilt  was  initially  oriented  in  the 
launcher  at  one  of  four  positions  spaced  90°  apart,  to  be  referred  to  as 
the  "3  o'clock",  "6  o'clock",  "9  o'clock",  and  "12  o'clock"  positions. 

For  these  rounds  data  were  obtained  on  the  motion  inside  the  9.5  ft. 
smoothbore  launcher  (neon  records),  on  the  motion  during  tip-off  and 
during  the  first  few  feet  of  travel  outside  the  launcher  (optical  lever 
records),  and  on  the  motion  at  and  near  the  end  of  burning.  For  a 
description  of  the  test  program  and  the  instrumentation  used  at  the 
Ordnance  Missile  Laboratories  of  Redstone  Arsenal,  see  the  report  of 
Horn  and  Cone  [HC],  The  reduced  data  from  the  program  are  presented  in 

IW-2]. 

The  theoretical  amount  of  linear  thrust  misalignment  produced  in 
the  M33  rocket  by  a  l°12r  nozzle-tilt  is  expressed  by  an  |  Lj  (linear 
thrust  misalignment  parameter)  of  at  least  0.026  ft;  see  the  report  of 
Cell  on  nozzle -plate  tolerances  [  C— It-  ] . 

This  amount  of  misalignment  was  sufficiently  large  so  as  to  be  a 
dominant  source  of  maltorque  even  though  it  might  be  combined  with  addi¬ 
tional  misalignment  or  unbalance  inherent  in  a  basic  standard  round.  One 
purpose  of  the  test  was  to  trace  the  effect  of  this  misalignment  upon  the 
motion  of  the  rocket  during  its  burning  period  in  an  attempt  to  determine 
whether  its  effect  upon  the  direction  of  motion  at  the  end  of  burning  was 
in  general  agreement  with  the  theoretically  predicted  effect.  Our  interest 
in  this  section  concerns  the  misalignment  effect  upon  the  rocket  motion 
inside  the  launcher. 


no 


Neon  Records 

In  Figure  6.5  are  shown,  side  by  side,  the  neon  records  of  the  rear 
band  contacts  of  three  rounds.  All  three  had  the  nozzle-tilt  initially 
oriented  at  the  3  o'clock  position.  The  sketch  below  shows  the  positions 
of  the  eleven  rails  in  the  segmented -rail  simulated-smoothbore  launcher 
as  viewed  from  the  rear  of  the  launcher. 


Fig.  6. hi  The  Positions  of  the  Eleven  Rails  of  the  Launcher 

The  clockwise  processional  patterns  of  the  rear  contacts  in  the  three 
records  are  very  similar.  All  show  the  final  contacts  (at  position  8  on  the 
records)  to  have  been  approximately  on  rail  k»  and,  if  one  traces  tho  con¬ 
tacts  back  through  the  launcher,  there  is  a  reasonably  def initio  nstarting 
position"  of  the  clockwise  pattern  back  in  the  second  section  of  the  launcher 
(after  about  inches  of  forward  motion).  This  starting  position  is  in  the 
neighborhood  of  rail  3,  and  this  corresponds  very  closely  to  the  initial 
3  o'clock  orientation  of  the  nozzle -tilt.  In  fact,  the  over-aid.  rate,  of;, change 
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of  orientation  of  contact  with  respect  to  distance  moved  down  the  launcher 
(i.e.,  a  mean  precessional  rate  in  radians  per  foot  of  travel)  is  very- 
close  to  the  spin  rate  of  the  rocket.  The  rockets  completed  one  spin- 
rotation  in  5  feet  of  travel. 

Similar  reproducibility-  of  motion  inside  the  launcher  was  found 
within  each  of  the  other  three  groups.  From  group  to  group,  the  main 
observable  differences  were  due  to  the  difference  in  the  initial  orienta¬ 
tion  of  the  nozzle-tilt.  This  is  shown  in  Figure  6.6,  in  which  is  presented 
a  representative  neon  record  from  each  of  the  four  groups.  Note  that  the 
"starting  positions"  of  the  precessional  patterns  correspond  closely  to 
the  initial  orientation  of  the  tilt.  Also,  note  that  the  positions  of 
final  contact  are  phased  approximately  at  90°  intervals,  and  that  the 
precessional  rates  are  quite  reproducible. 

(2)  A  Mathematical  Model  And  Theoretical  Considerations 

Figures  6.6  and  6.6  indicate  a  significant  dependence  of  the  preces¬ 
sional  pattern  of  contact  between  rear-band  and  launcher  upon  the  orienta¬ 
tion  of  the  nozzle -tilt  as  the  rocket  moved  down  the  launcher.  In  these 
rounds  the  nozzle-tilt  corresponded  to  an  intentionally  large  thrust  mis¬ 
alignment  which  can  also  be  interpreted  in  terms  of  a  cross-force  of 
approximately  70  lb.  acting  at  the  rear  band.  In  an  attempt  to  set  up 
a  simple  mathematical  model  to  account  for  the  observed  precessional  motion 
of  the  rear  of  the  rocket  -under  the  action  of  such  a  large  cross-force,  we 
have  produced  a  reasonably  simple  mathematical  representation  which  has 
broader  implications  than  those  which  were  being  sought.  Not  only  did  the 
resulting  mathematical  theory  call  attention  to  some  aspects  of  the  above- 
mentioned  records  whose  significance  had  not  been  previously  appreciated, 
but,  in  addition,  simple  modifications  became  evident  which  provided  a  basis 
for  understanding  other  patterns  of  contact  of  both  front  and  rear  bands 
which  had  been  observed  in  many  neon  records. 

There  have  been  previous  theoretical  studies  of  the  motion  inside  a 
smoothbore  launcher,  and  general  equations  of  motion  have  been  set  up. 

See  Hall  [Ha]  and  Cell,  Herz,  Menius  ICHM].  None  of  these  has  yielded 
a  satisfactory  understanding  of  the  motion  as  depicted  by  the  neon  records. 
One  such  representation  was  presented  by  Herz  in  Chapter  VI  and  Appendix  I 


of  the  reference  [GHM],  He  pointed  out  that  because  of  certain  physical 
characteristics  of  such  rockets  as  the  Ml 7  and  M33  it  is  possible  to 
approximate  the  physical  situation  in  terms  of  two  uncoupled  particles, 
one  at  the  center  of  the  rear  band  and  one  at  the  center  of  the  front 
band.  In  such  a  representation,  one  would  have 


2  2  u  2 

m  =  mk  /(k  +  JL  )  =  the  'reduced  mass  at  the  rear  band' 

r  r 

m^  =  mk^/(k^  +  =  the  'reduced  mass  at  the  front  hand1. 


where 


=  B  =  transverse  moment  of  inertia  about  the  c.g.. 


=  distance  from  c.g.  to  center  of  rear  band, 

=  distance  from  c.g.  to  center  of  front  band. 


a)  We  wish  to  use  such  an  approximate  representation,  but,  as  an  aid 
to  visualization,  we  shall  replace  the  two  uncoupled  particles  by  un¬ 
coupled  solid  cylindrical  disks  of  appropriate  mass  and  with  diameters 
equal  to  the  respective  band  diameters. 

b)  We  shall  assume  that  frictional  effects  arising  from  contact  bet¬ 
ween  the  disks  and  the  launcher  can  be  ignored.  Thus,  we  picture  the 
disks  as  sliding  freely  on  the  inner  surface  of  the  launcher. 

c)  We  shall  be  primarily  concerned  with  the  nature  of  'continuous  con¬ 
tact'  patterns  between  disk  and  launcher,  although  it  is  possible  also 

to  consider  conditions  which  may  cause  a  discontinunity  in  the  contact. 

d)  Many  neon  records  have  indicated  a  sudden  lift  of  the  rear  of  the 
rocket  (and  to  a  lesser  extent  of  the  front  also)  immediately  after 
ignition.  This  has  been  analyzed  in  [CHM]  in  terms  of  a  Bernoulli 
effect.  In  numerous  cases  this  lift  has  resulted  in  the  rear  of  the 
rocket  striking  the  top  of  the  launcher.  In  some  cases,  this  leads  to 

a  balloting  type  motion  of  the  rear  as  the  rocket  moves  down  the  launcher 


More  often,  as  the  rear  falls,  it  'slides'  down  a  portion  of 
wall  and  there  is  set  up  a  'continuous  contact'  motion  with 
cross-- velocity  of  significant  magnitude.  Such  a  record  is  s 
Figure  6.9.  We  shall  be  interested  in  the  above  lift  effect 
source  of  significantly  large  initial  cross-velocities. 

The  Basic  Model 
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0:  center  of  outer  circle  (i.e.,  on  launcher  axis) 
G:  center  of  sliding  disk 
P:  point  of  contact 

0C  =  one  half  of  clearance  =  g,  (e.g. ,  0.01  in.) 
CP  =  radius  of  disk,  (e.g.,  2.29  in.) 

OC  =  complex  variable  p(t)£=  e^  ^ 

?  =  a  force  rotating  with  the  rocket  =  Fe^^ 


Considering  the  forces  acting  on  the  disk,  we  obtain  in  complex 

form 

nip*  =  -  mg  +  Fe^  -  Ne^  .  (6,1) 

With  p  =  £ e^  ( £  assumed  constant  for  continuous- contact),  we  obtain 
as  component  equations  along  OP  and  normal  to  OP, 

K  g  sin  X  +  (F/m)  sin  (p  -X) , 

£ 

£>  l.  mg  cos  X  -(F/m)  cos  (p-X)  +  N/m.  (6*2) 

Three  cases  will  be  considered  briefly.  There  is  first  the  case  where 
F  is  negligible  in  comparison  with  the  gravity  effect.  Secondly,  we  consider 
F»mg,  constant  in  magnitude,  but  spinning  with  the  rocket.  Lastly,  an  F 
of  variable  magnitude  can  arise  from  unbalance  (dynamic  and/or  static)  in  the 
rocket,  with  the  magnitude  of  F  increasing  as  the  square  of  the  spin-rate. 

Case  1:  F  =  0. 

In  the  sketch  of  our  model,  it  is  convenient  in  this  case  to  orient  P 
relative  to  the  bottom  of  the  launcher,  usings  =  180°  -X.  Equations  (6.2) 
become 

£  V+  g  sin  \l*  =  0 

K  t+  g  cos  t  =  N/m  (6.3) 

The  first  of  these  is  the  familiar  pendulum  equation  which  may  be  integrated 
to  yield  the  energy  equation, 

l/2(Sf)  +  2g  g  sin2  |  =  1/2 (?tQ)2  +  2g  ?  sin2  =  Eq,  (6.10 

in  which  Eq  has  the  dimensions  of  energy  per  unit  mass.  The  second  equa¬ 
tion  may  be  written  in  the  form 

N  =  — jr  +  |g  K  -  3g  K  sin2  || .  (6.5) 

In  order  to  insure  continuous  contact,  N  must  be  non-negative.  Two  such 
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possible  cases  are  of  interest. 

a)  If  Eq  =  g£,  one  has  the  oscillatory  pendulum  type  motion  with 
amplitude  =  90°.  The  period  of  such  oscillations  is  dependent  upon  g, 
and  the  amplitude.  For  the  neon  records  shown  in  Figures  6.8,  6.9 , 

S,  =  l/2(0. 02)  in.  =  0.0008  ft.  and  g/£  =  l;0.000(sec.-^).  In  order  for  such 

an  oscillation  to  show  up  on  a  neon  record  there  must  be  an  amplitude  of 
at  least  17°.  Under  these  circumstances  the  period  would  range  from  32 
milliseconds  up  to  37  ms.  for  a  90°  amplitude. 

Such  oscillations  appear  primarily  in  the  neon  records  of  front 
band  contacts.  In  Figure  6.8,  such  an  oscillation  is  seen  in  the  nose  con¬ 
tact  record  with  a  period  of  approximately  35  ms.  There  is  also  an  indica¬ 
tion  of  such  an  oscillation  in  the  front  contact  record  shown  in  Figure  6.9. 

b)  If  Eq  =  5g€/2,  one  has  a  'circulatory'  type  motion  in  which  the 
disk  (in  the  model)  would  slide  completely  around  the  launcher  in  a  periodic 
fashion.  With  §  =  0.0008  ft.,  (i.e.,  0.02  inch  clearance  in  the  launcher), 
one  finds  that  if  Eq  is  interpreted  in  teres  of  the  preoessional  rate  (clock- 
wise  or  counterclockwise)  at  the  bottom  of  the  launcher,  the  minimum  such 
precessional  rate  required  for  a  'continuous  contact1  circulatory  motion  is 
approximately  U5o  rad. /sec.  with  a  corresponding  maximum  period  of  approxi¬ 
mately  20  ms. 

Such  motion  has  been  observed  many  times  in  the  neon  records  of 
rear  band  contact.  In  most  cases,  it  appears  as  a  counterclockwise  pattern 
(i.e.,  as  viewed  from  the  rear  of  the  launcher).  Figure  6.9  shows  such  a 
counterclock  motion  of  the  rear  of  the  rocket  with  fairly  uniform  period 
of  about  10  ms.  This  type  of  neon  record  usually  indicates  that  the  rear 
of  the  rocket  has  experienced  the  Bernoulli  lift  effect  and  has  'bounced' 
off  the  top  of  the  launcher  (once  or  twice)  and  has  acquired  a  significant 
cross-velocity  as  it  goes  into  the  circulatory  motion. 

c)  If  g £<  Eq<  (5/2) gg  ,  the  disk  has  sufficient  energy  to  slide  above 
the  90°  level  (on  either  side)  but,  at  some  point  before  it  reaches  the  top  .. 
of  the  launcher,  the  reaction  force  N  will  reduce  to  zero,  and  contact  will  . 
be  broken  and  a  'continuous  contact'  type  of  motion’  will  not  prevail.  - 


,V3¥g3SfSa®^385? 
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Case  2;  F  =  a  constant»mg  and  -with  F  rotating  with  the  rocket. 

As  a  first  approximation  to  this  case,  we  use  equations(6, 2)  with  the 
gravity  terms  omitted.  Since  F  rotates  with  the  rocket,  the  angle  p(t)  is 
such  that  ft  =  ?)  and  fc  =  Vf,  where  ?)  and  rf  are  the  spin  rate  and  angular 
acceleration  respectively.  These  are  assumed  to  be  known,  and  in  the  M33 

p 

case  rf  is  taken  to  be  a  constant  (approx.  3000  rad. /sec.  ).  Letting 


0  =  X  -  n,  (6.6) 

the  first  of  equations  6-2)  is  expressible  approximately  as 

0  +  ■—  sin  0  =  -  M-  (6.7) 

mg 

With  F  ~  70,  m~0.28,  and  g  =  0.0008,  corresponding  to  the  cases  depicted 
in  Figures  6.5  and  6.6,  we  find  that 

F/m  £**  300,  000, 

and  thus  we  choose  to  ignore  the  tern  -  ±iequation(6.7),  and  consider 

0  +  sin  0=0.  (6,8) 

mg 

This  is  the  pendulum  equation  again.  Because  the  size  of  F/mg  (as 
indicated  above)  would  tend  to  prohibit  a  circulatory 'type  solution  •under, 
the  conditions  of  our  rocket  problem,  we  consider  only  the  oscillatory  type 
of  solution  with  amplitude  <  90°. 

Note  that  an  oscillatory  solution  for  0  would  have  the  physical  signi¬ 
ficance  of  an  oscillation  of  the  'contact -angle1  about  the  'prescribed' 
spin  angle  j i;  that  is,  in  a  neon  record,  an  oscillation  of  the  contact  about 
the  rotating  orientation  of  the  cross-force.  'For  small  0  oscillations,  the 
'  period  would  be  given  by 

-  . T  “  2ir/V  F/m  g . (  ~11  ms.  for  values  cited-  above).  . 

For  oscillations  of  larger  amplitudes  the  periods  would  be  slightly 

larger.  Y,.v'  '::Vw.Y  ./;Y 
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One  can  see  some  evidence  in  Figures  6. $  and  6.6  of  oscillations  of 
the  rear  contact  about  some  rotating  'equilibrium  position'.  This  type 
of  motion  shows  up  very  clearly  in  the  record  of  Round  5U 1  as  shown  in 
Figure  6.8.  Round  SUl  was  one  of  the  rounds  with  the  nozzle -tilt  initially 
at  9  o'clock.  The  oscillation  shown  in  the  lower  portion  of  that  figure 
was  obtained  by  subtracting  from  the  contact  record  the  spin -angle  p(t) . 

The  resulting  approximate  curve  clearly  shows  an  oscillatory  pattern  about 
the  initial  9  o'clock  position  and  its  'period'  is  approximately  equal  to 
the  11  ms. 


Case  3:  Effects  of  Unbalance 

In  order  to  represent  the  unbalance  effects  in  equations (6. 2),  it  is 
necessary  to  introduce  two  unbalance  parameters  relating  to  the  rear  and 


front  bands  respectively: 

U 

r 


r 

c 


> 


where  p  =  the  measure  (in  radians)  of  dynamic  unbalance, 
c 

rc  =  the  measure  (in  feet)  of  static  unbalance. 


In  the  'uncoupling'  of  the  equations  of  motion  mentioned  above,  it  is  found 
that  the  corresponding  forces  on  the  'uncoupled'  disks  at  the  rear  and 


front  are  given  by 


=  -m  U  co 


Ff  = 

2 

if  ®(the  angular  acceleration)  is  small  in  comparison  with  co  (the  square 
of  the  spin  rate). 

If  dynamic  unbalance  is  dominant,  then  these  two  forces  are 

essentially  180°  out  of  phase  in  their  orientations.  In  the  firing  pro-  - . 

gram  referenced  above,  some  rounds  were  deliberately  unbalanced  dyna¬ 
mically  sc  as  to  produce  a  Jp  |  ~;2.2  mils  without  magnifying  the  static 
unbalance.  With  a  clearance  of  0.02  in.,  with  ^  =  (5A)  ft*  and 
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/'  =  (5/12 )  ft.,  one  obtains  (with  r  =0)  the  following  orders  of 

X  C 

magnitudes : 

At  rear  band  At  front  band 

U  ~  0.0028  ft.  U_  ~ 0.0009 (ft.) 

r  f 


Urn fe  ~  21,000  X 


~ 7,000  X(sec.-2), 


where  X  denotes  the  distance  (ft.)  through  which  the  band  has  moved  down 
the  launcher. 

2 

Thus,!  at  the  r  ear  band,  U  ca  /g  increases  from  zero  to  approximately 

_2 

200, 000 ( sec. ~  )  when  the  band  is  at  the  muzzle  of  the  9-5  ft.  launcher. 

After  two  feet  of  travel,  U^co  ~  i|2, 000  and  is  comparable  in  magnitude 
with  the  corresponding  gravity  parameter  (namely,  g/t  =  lj.0,000).  After 

p 

U  or  5  feet  of  travel,  U  co  definitely  dominates  the  gravity  effect  at 
the  rear  band. 

2  —2 
At  the  front  band,  U^xo  /£  increases  only  to  5U,000(sec.“  )  as  the 

front  band  reaches  the  muzzle  after  travelling  7-75  ft.  Only  after  some 

5  feet  of  travel  does  U^n  /  £  become  comparable  with  g/^  .  However,  the 

following  is  worth  noting.  If  the  contact  at  the  front  band  remains  near 

the  bottom  of  "the  launcher  under  the  influence  of  gravity  (with  X  near 

180  ),  then  in  the  first  of  equations  $.2)  as  F/mg  (i.e.  U^o  /£)  becomes 

comparable  in  magnitude  with  gfe  and  its  orientation  (spinning  with  the 

rockets)  causes  |p  -X|to  tend  toward  90°,  the  unbalance  term  can  become, 

at  least  temporarily,  dominant.  This  can  cause  the  contact  to  move  toward 

the  orientation  of  the  unbalance  force,  and  with  this  force  increasing 

in  magnitude,  it  can  affect  the  pattern  of  the  front  contact  during  the  last 

two  or  three  feet  of  travel.  This  effect  has  been  noted  often  in  the  neon 

records  of  the  front  band  contacts.  There  is  some  indication  of  such  an 

effect  in  the  front  band  contact  record  shown  in  Figure  6.8,  with  a 

discontinunity  in  contact  following  a  simple  gravity  oscillation.  - 


Figure  6.10  shows  the  neon  records  of  one  of  the rounds  which  were 


dynamically  unbalanced  with  |J3  |  ~ 2 . 2.  Note  that  at  the  rear  the  contact 
remained  near  the  bottom  of  the  launcher  (rails  No.ll  and  No.l  are- at  the- 
bottom  of  the  launcher-)  during,  the  first  2v5  feet  of  travel  to  beyond  posi¬ 
tion  3.  By  that  time  the  unbalance  effect  had  become  significant,  and 
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starting  at  position  U,  there  is  a  continuous  contact  dominated  by  the  un¬ 
balance  force.  Note  the  oscillation  of  the  contacts  about  the  indicated 
rotation  of  the  unbalance  force  which  was  initially- directed  toward  the 
bottom  of  the  launcher,  and  which  had  rotated  again  ‘£0  that  direction  just 
beyond  position  3  in  the  record  (after  feet  of  travel) .  At  the  front 
band,  the  contacts  show  a  simple  gravity-type  oscillation  about  the  bottom 
of  the  launcher  up  through  position  C.  From  position  D  through  F,  the  con¬ 
tact  is  dominated  by  the  unbalance  force  at  the  front  band.  Note  that  the 
indicated  rotating  orientation  of  the  unbalance  force  at  the  front  band  is 
180°  out  of  phase  with  that  at  the  rear  band. 

Mathematically,  the  equations  of  motion  arising  from  equations  (6,2)  dn 
this  case  are  not  as  simple  as  in  Cases  1  and  2  since  the  magnitude  of  F  is 
not  constant  for  an  accelerating  rocket.  However,  qualitatively  the  ideas 
of  Case  2  apply  after  the  unbalance  force  becomes  the  dominating  effect. 
Instead  of  oscillations  periodic  in  time,  one  can  easily  show  that  if  the 
rocket  has  constant  linear  and  angular  accelerations,  there  is  an  approximate 
periodicity  in  the  distance  variable  X.  The  first  of  equations  (6.2)  is  trans 
formed  into  an  equation, 

0"  +  (l/2X)0*  +  (Un2/t)  sin  0  =  0, 

comparable  with  equation (6 . 8 )but  with  differentiation  with  respect  to  X. 

In  the  coefficient  of  sin  8 ,  n  denotes  the  constant  ratio  of  angular  to 
linear  acceleration.  Thus,  for  small  0,  approximate  solutions  are  ex¬ 
pressible  in  terms  of  Bessel  functions  of  orders  -  (l/lj.).  The  corresponding 
oscillations  would  not  be  strictly  periodic  and  would  be  subject  to  some 
damping.  But,  for  the  neon  records  obtained  in  the  test  program : referenced 
above,  the  data  could  be  expected  to  exhibit  at  most  two  successive  oscilla¬ 
tions  since  the  expected  period  can  be  estimated  to  be  nearly  3  feet.  In 
Figure  6.10,  the  distance  along  the  launcher  between  positions  3-  and  7  is  3 
feet. 

With  regard  to  the  three  cases  which  have  been  briefly  presented,  it 
should  be  kept  in  mind  that  the  fundamental  equations(6.2)  are  non-linear, 
and  thus  one  cannot  think  in  terms  of  a  linear '  superposition  of  ^effects. 

The  above  treatment  sheds  some  light  on  what  may  be  expected  if -some  one 
disturbing  factor  is  dominant.  =  1  ' 


APPENDIX  A 


ASYMPTOTIC  ESTIMATES  OF  THE 
DEVIATION  OF  BOOSTED  FIN -STABILIZED 
ROCKETS  DURING  THE  BURNING  PERIOD 


In  this  appendix  the  mathematical  procedure  for  arriving  at  the  results 
given  in  Section  5-5  will  be  presented.  The  techniques  involved  here  are 
somewhat  unusual,  and  the  implications  relative  to  applications  in  other 
areas  should  be  significant. 

The  take-off  point  for  the  discussion  here  is  the  fundamental  system 

of  differential  equations  given  in  Chapter  5  as  equations  (5-1.7)  and  (5.1.8) 

in  which  $  is  measured  relative  to  the  ideal  trajectory.  Again  M  and  F 

c  c 

are  given  by  (5.1.5)  and  (5.1.6)  and  other  notation  is  the  same  as  in  Chapter 
5,  except  for  the  fact  that  in  the  treatment  the  ratio  n  of  spin  rate  to 
linear  velocity  is  allowed  to  vary.  This  extends  the  results  obtained  here 
to  cases  covered  by  the  more  general  theory  given  in  [CH],  and  the  results 
may  readily  be  specialized  to  the  case  of  constant  n  as  was  done  in  writing 
the  formulas  of  Section  5«5- 

By  changes  of  independent  variable  in  equation  (5.1.7)  and  (5-1.8),  the 
following  basic  equation  result: 

•  "  • 

-  ^iqn-c.-J#  +  c.,(vA  +  w  )=  -iM  /Bv 

^  H  M  c  c 

(A.l) 

2  2 

-  [(g  cos  e)/v][2iqn-cH  +  v/v  -(g  sin  e)/v  ], 


-  (vA)'  -  c„(vA  +  w  )  =  F  /mv, 

N  C  C 

wherein  the  dots  indicate  differentiation  with  respect  to  time  t,  and  the 
primes  indicate  differentiation  with  respect  to  trajectory  arc  length  s, 
that  is. 


i.  =  JL  («) 

ds  W 
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Equations  (A.l)  yield,  on  elimination  of  #,  a  single  differential 
equation  in  (vA).  The  corresponding  homogeneous  (or  reduced)  equation 
obtained  from  this  new  equation  in  (vA)  is  the  following: 

(vA) "  -  (2iqn  -  cH  -  c^HvA)'  +  (cM~  2iqn  cn)(tA)  =  0.  (A.  2) 

In  the  case  of  no  spin  (i.e.,  n  =  0),  if  the  damping  of  (cH  +  c^.)  is 
ignored,  (vA)  is  a  periodic  function  of  distance  with  a  “period”  (i.e., 
the  wave-length  of  yaw)  equal  to  2tt/  If  X  denotes  this  wave-length 

of  yaw, 

X  =  2tr/  Vc^" .  (A. 3) 

This  basic  wave-length  is  still  significant  even  when  slow  spin  and  damping 
are  taken  int o  account. 

Likewise,  if  only  the  gravity  forcing  term  is  retained  in  the  equation 
resulting  from  eliminating  4  in  equations  (A.l),  and  if  A  is  taken  as  the 
dependent  variable,  one  obtains  the  differential  equation 

A”  -  (2iqn  -  CH  “  CN  -  1  +  lew  -  2iqn(cN  +  -^)]  A  = 

v  v 

g  cos  6/„.  .  G  2g  sin  e\  /.  i  \ 

-  £>~2 - (2iqn  +  ~  "  CH  ~  cd - 2 - )•  (A-h) 

v  v  v 

This  equation  yields,  as  a  good  approximation  for  the  yaw  of  repose, 

—2  —2  —2 
-g(cos  e)v  [2iqn  +  Gv  ”  ch  ”  CD  ”  2g(sin  e^v  ] 

A  —  n  __o  • 

cM  -  2iqn[Gv”  +  Cjj  -  ^  -  g(sin  e)v“  ] 

2 

This  differs  from  that  given  in  [MKR]  only  in  the  presence  of  G/v  ,  expressing 

2 

the  effect  of  the  rocket  thrust  in  v/v  . 
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The  gravity  effects  will  not  be  further  considered  in  this  appendix. 

It  can  be  assumed  that  in  equations  (A.l),  after  v  and  co  are  replaced  by 

v(s)  and  co(s)  as  obtained  for  the  approximate  trajectory,  the  forcing 

terms,  M  /Bv  and  F  /mv,  expressing  the  misalignments,  will  be  independent 
c  c 

of  further  significant  gravity  effect.  Thus,  in  solving  the  linear  system 
(A.l),  the  superposition  principle  enables  one  to  consider  separately  the 
effects  of  initial  launch  conditions,  of  the  misalignments,  and  of  gravity. 
For  a  given  type  of  rocket,  to  the  extent  that  v  ,  coq,  and  e  can-  be 
reproduced  from  round  to  round,  the  gravity  effect  upon  the  deviation 
will  be  reproducible  and  is  not  a  source  of  dispersion.  The  parameters 
associated  with  the  other  effects,  (initial  yaw  Aq,  initial  cross-spin 
initial  deviation  @q,  and  the  misalignment  parameters,  L^,  p^,  and  cO 
are  not,  in  general,  reproducible  from  round  to  round,  either  in  magnitude 
or  in  orientation,  and  need  to  be  considered  as  possible  sources  of 
dispersion. 

We  shall  make  the  following  basic  assumption  regarding  stability. 
Discussions  of  relevant  statility  conditions  may  be  found  in  [MKR], 

[DFB],  and  iMu]. 

•  - 

Stability  Assumption:  That  the  solution,  and  (v A),  of  the  homo¬ 
geneous  system  resulting  from  (A.l),  namely 

(  * 

-  (2iqn  -  c^)#  +  c^C vA)  =  0, 

(A.  6) 

L  $  -  (vA) 1  -  cn(vA)  =  0, 

•  • 

as  produced  by  initial  conditions,  and  Vq Aq,  are  such  that  and  (vA) 
remain  bounded,  with  A  and  $/v  damping  out  as  trajectory  distance,  s, 
increases,  the  damping  being  produced  by  an  increasing  v,  or  aerodynamic 
damping,  or  both. 

In  the  considerations  to  follow,  the  solution  of  the  homogeneous  system 

(A. 6)  will  be  expressed  in  terms  of  two  basic  solutions  determined  by  the 

•  • 

two  sets  of  initial  conditions:  ($  =  1,  A  =  Q)  and  (§  =  Q,  A  =  l). 

o  o  ,  o  o 
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Thus  a  general  solution  of  equations  (A. 6)  is  to  be  expressed  as 


(  • 


§  =  $  §  (s  ,s)  +  A  3?_(s  ,s), 

O  cj  N  o  0 v  o 7  7 


A  =  ^Ajs^s)  +  +  Ag(s0,s)], 


o  q  o 


in  which  the  subscripts,  q  and  6,  are  used  to  denote  respectively  the  first 

and  second  of  the  above  sets  of  initial  conditions.  Note  that  the  functions 

•  • 

3>g,A^,  and  Ag  have  initial  values  respectively  of  1,  0,  0,  and  0.  In 
turn,  one  can  produce  functions,  and  #g,  with 

r>t 


3?  =  4?  + 
o 


[  4?  dt  =  #  +  —  ds  =  §  +  4?  #  (s  ,s)+  A  ^c(s  ,s), 

Jt„  o  i  V  o  c  q  o’  O  6  o’  (A_8) 


in  which  the  functions  4?  and  4?,  are  initially  zero.  Likewise, 

q  o 

• 

0=  #  -  A=  0  +  #  0  (s  , s)  +  A  0K(s  ,s),  (A. 9) 

o  o  q  o  o  o  o  * 

in  which  0  =  $  -  A,0  =  4>  -A,  and  0C  =  <£..  -  Ac.  Note  that  bothJS) 

o  o  o’  q  q  q*  o  o  6  -  q 

and  are  initially  zero,  with  the  initial  yaw  A^  having  been  made  a 
part  of  the  initial  parameter  0q. 

We  shall  next  obtain,  under  conditions  applicable  to  rockets  launched 
at  high  velocity,  asymptotic  estimates  for  0^  and  0^.  Then  later  some 
approximations  will  be  derived,  expressing  the  effects  of  the  forcing 
terms  (cross-wind  and  misalignments)  in  terms  of  "equivalent  initial 
conditions"  and  thus  the  corresponding  effects  upon  angular  deviation 
®  will  be  expressible  in  terms  of  0^  and  0^. 

Asymptotic  Estimates  for  _0^  and 

WeTIow  inpose  restrictive  conditions’  on  the  launch  velocity  and  the 
acceleration  due  to  rocket  thrust  outside  the  launcher.  It  will  be  assumed 

that,  during  the  burning  period  outside  the  launcher, 

/ 

v  =  3>00  ft. /sec., 

(A.  10) 

f/r2  ]  =  10“3  ft.-1  . 

t  -■ 

It  will  also  be  necessary  to  know  something  about  the  spin  in  order  to 

I 
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estimate  the  magnitudes  of  terms  in  which  n  (i.e.,  co/v)  appears.  ■  For 

the  cases  where  spin  is  involved,  we  shall  assume  that  the  spin  rate  is 

such  that  lOc^  =  n  <  0.23 (rad. /ft. ) .  The  upper  bound  on  n  simply  means 

that  the  spin-rate  is  not  permitted  to  attain  the  high  values  characteristic 

of  ospin-stabilized  projectiles  (for  which  normally  1<  n<  2).  A  word  of 

explanation  is  in  order  concerning  the  lower  bound.  Recalling  that  the 

wave-length  of  yaw  is\=  2Tr/yrc^  ,  we  note  that,  for  n  constant,  distance 

transversed  by  the  rocket  during  one  spin-revolution  of  the  rocket  is 

expressible  as  (2-rr/n)  ft.,  and  thus  the  number  of  spin -revolutions  per 

wave-length  of  yaw  is  given  by  n/  yG^.  To  avoid  the  possibility  of 

resonance  effects,  introduced  by  forcing  terms  which  rotate  with  the 

rocket,  it  is  desirable  to  avoid  having  n/  y£^  close  to  unity.  It  is 

expected  that,  if  spin  is  employed,  then  n  would  be  maintained  at  values 

2 

exceeding  yE^.  The  condition  n  =  10c^  means  that  the  rocket  will  complete 
at  least  3  spin-revolutions  per  wave-length  of  yaw.  With  cM~  0.001,  we 
shall  consider  n  to  be ~  0.1  or  greater. 

It  is  also  expected  that  the  possible  variation  in  n  will  tend  to  be 
either  monotanic  increasing  or  decreasing  during  the  burning  period.  With 
the  assumption  of  a  burning  distance  an  the  order  of  1000  ft. ,  we  shall 
thereby  expect  that  |  n'  |  will  not  exceed  5(10-^)  and  that  |  n"  |<10~^. 

The  three  most  specific  and  simplest  patterns  are  respectively: 

(1)  no  spin,  with  a>  =  n  =  0; 

(2)  pre-spun,  with  only  aerodynamic -spin  deceleration  outside  the 

launcher,  with  oo  =  co  -Lc.cods: 

o  A 

(3)  spin  proportional  to  velocity,  i.e.,  n  constant. 

For  all  of  these  cases,  we  can  write 

n'  =  (db/v2)  -  n(v/v2)>  (A.ll) 

with  this  quantity  reducing  to  0  in  cases  (l)  and  (3)  and  to  (-n)  (c  +  v/v2) 

\  & 
in  case  (2). 

From  the  second  of  equations  (A. 6), 

v$'  -  vA  '  =  v©  '  =  (Cjj  +  v/v^HvA), 

and  thus  ‘  ‘  (A. 12) 

©  1  =  (cjj  +  Vy2)  a. 
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Under  our-  stability  assumption,  A -HD  and  (cR  +  v/v  )  remains  bounded  as 
S-+  00 .  Hence  (0  -  @  )  will  approach  a  limiting  value  which  is  also  the 
limiting  value  of  (§  -  +  Aq).  It  is  convenient  to  seek  an  estimate  of 

the  limiting  value  by  considering  the  relation 

*  -  =  f  |  ds.  (A.13) 

so 

The  homogeneous  system  (A. 6)  can  be  rewritten  in  the  form 
4>"  -  (2iqn  -  cR  -  cR)$'  +  (cM  -  2iqn  cw  -  2iqn’)3>  =  0, 


#'  =  (2iqn  -  cR)#  -  Cjj(vA). 


M 


N 


Let  L  denote  the  linear  operator 

L  =  IU2  -  (2iqn  -  cR  -  cn)D  +  (c^  -  2iqn  cN  -  2iqn')], 
and  M  its  adjoint  operator 


(A.lir) 


(A.  15) 


M  »  ID  +  (2iqn  -  cR  -  cN)D  +  (cM  -  2iqn  cN)]. 

A  standard  Green's  formula  may  be  expressed  as  follows: 

r®  •  •  •  •  • 

J  lL(#)f(s)  -  M(f)3>]ds  =  If  #'  -  f'§?  -  (2iqn  -  c^-Cjj)^] 
so 


(A. 16) 


=  HcNf  -  f')^  -  cM(vA)f]g 

o5 


(A.  17) 


in  which  f  denotes  an  arbitrary  function  such  that  the  integral  in  (A. 17) 
is  meaningful. 

To  evaluate  the  integral  of  (A.I3),  we  wish  to  have  m(f)  =  -l/v. 

If  one  designates  that 


*00-^7 

CM  CN  ”  v 


(° u  +  cT  +  3G/v  ) 


°M 


ft)  3. 


(A.  18) 


in  which  cR  =  cR  -  cR,  then  the  more  significant  terms  of  f  and  f"  are 
given  by 

*  ^ + ^  *  3  ^ 

M  v  M  v  v  m  v 
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1  {SJ  c Mvv  2;  * 

M  v 


(A. 20) 


Considering  M(f),  one  finds  that,  within  the  conditions  which  have  been 
imposed, 

M(f)  =  -  +  Q(10"2)]« -1/v.  (A. 21) 


•  ~j 

Since  L($)  s  0  and  M(f)  =  -  one  obtains,  from  (A. 17)  and  (A.13), 


$  -  4> 


-  r±  as. 

J  S„  V 


l(cNf  -  f')§  -  cM(vA)f] 

o 


(A.  22) 


Using  (A. 18)  and  (A. 19),  one  finds  that  the  more  significant  terms  of 


(c^f  -  f ' )  are  as  follows: 


CNf-f’  (cN+~)  + 


v  v  ^  2iqn 


(A.23) 


and  thus  that 


o 


4^*4>-(nr3SA)^  I 


(A.2U) 


Under  the  stability  assumption,  #/v  and  A  approach  zero  as  s  increases. 
Thus  from  (A. 2b)  we  find  that 

$  [  v  2iqn'  2iqn  v  _  v  0  1 

-  *o>”  ^  {  (0,  +  — S>CL  -  —Si*  — -2  V  *  *«-?>*] 


!2iqn  ^  (c„  +  cT  +  3G/v  ^)  ^  j 

1  +  ^  - ^r>  •  tt-ao 

M  V  M  V 

o  o  J 


132 


and 

1 3 m  ( ©  -  @  )  =  lim(#  -  #  +  A  ) 

O  0  0 


# 


C„V 

M  o 


V  2iqrL'  1  Ir  f 

v  M  v  v 

o  oo 


-A 


2iqn  ^  (c^  +  c,  +  3G/v  )  v 

2(%  *  -2*)  -  -S— i - Gs_(  o  ) 


c„  '"N  2 

M  v 

o 


'M 


■T' 

V 

o 


(A.  26) 


We  thus  obtain,  as  the  asymptotic  estimates  for  andQg,  the 
following: 

,  v  2iqn  *  2iqn  v  „  v  l 

l  (V  -T>&  -  -S7->  . 

v  MM  v  v  J 

o  o  o 

(A. 27) 
(A.  28) 


q  oJ 


owv 
M  o 


2iqn  t 

VV°  ’  --^-2  (o»  * 

M  v 

o 


(ch  +  CT  +  3G/v^)  ir 


"M 


(_2_) 


in  which  vq  and  n^(i.e.,  no/vQ)  are  to  be  interpreted  as  t-»-ro  +  or 
physically  simply  as  the  rates  of  change  just  outside  the  launcher, 
assuming  that  the  rocket  thrust  (as  opposed  to  that  of  the  boost  charge 
in  the  gun  tube)  has  already  been  initiated  at  launch.  However,  if  the 
rocket  propellent  is  not  ignited  until  after  the  rocket  has  traversed  a 
significant  distance  along  its  trajectory,  then  v  ,  n^,  v^,  and  hQ  are 
to  be  evaluated  at  that  time. 

For  the  case  of  no  spin,  one  notes  that 


Vv”5  '  (°N  +  VT0  )/oMTo 


®8(so>“)  ‘ 


(cH  +  cT  +  3G/v^  )  -V 


H 


2-(-4) 


(A.  29) 


o  / 


For  the  case  of  the  pre-spun  rocket  with  no  angular  acceleration  after 

2 

launch  except  aerodynamic  deceleration,  (with  ca  of  lower  order  than  vq/vo  ) 


®  (s  »)  = 


’q  '  o'  '  c „v 


v  2iqn  _  v  v  „ 

(c  +  _2_)  +  — °  r.  2.  o  ,i./  o^2 

'CW  O' 


Mo 


"N  2J 
v 

0 


"M 


rV+3cN“T  *M-2zr)  >  fc* 


Y, 


with  ®-(s  ,“).as  in  (A. 28). 
o  o 
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For  the  case  of  n  constant,  we  obtain  the  estimates  for  ©  and@_ 

- i  q  8 

directly  from  (A. 27)  and  (A. 28),  setting  nQ  =  0  in  (A. 27). 

In  the  cases  where  the  rocket  is  spinning,  one  notes that  the 
estimate  is  predominantly  real  and  positive.  As  viewed  from  the  rear 
of  the  launcher, a  (initial  cross-spin)  which  is  directed  upward  produces 
a  limiting  angular  deviation©^  which  is  upward  and  slightly  to  the  right. 
Similarly  the  0  g,  due  to  an  initial  upward  yaw,  is  .predominantly  imaginary  and 
is  directed  to  the  left  and  slightly  upward. 

In  the  case  of  no  spin,  both  estimates  are  real  and  positive. 

If,  in  the  cases  where  spin  is  present,  one  is  Interested  primarily 
in  the  magnitudes  of  ©^  andQg  then  the  estimates  can  be  given  in  the 
simpler  forms, 


®  a (c„  +  v  /v  2)/c„v  , 
q  N  o'  o  "  M  o> 


(A.31) 


@o  a-  2iqn  (c„  +  v  /v  2)/cM  =  -  2iqoo  0  .  (A. 32) 

o  ^  o  N  o'  o  'M  oq 

_  2 

In  (A.  31)  if  one  replace  v  by  (G  -  C-v  -  g  sin  e  )  and  employs  the  "lift" 


parameter,  cp  =  “  cpi  then 


g  sm  e  p 

q  «-L_  (c - o  +  G  ). 

a  c.,v  L  2  2 

^  Mo  v  v 

o  o 


(A.  33) 


The  one  significant  respect  in  which  (A. 33)  and  the  corresponding  estimate 

for  ©6  differ  from  the  analogous  estimates  for  ordinary  projectiles  (as 

2 

derived,  for  instance  in  [MKR])  is  the  presence  of  the  term.  G/v  arising 

p  o 

from  the  rocket  thrust.  When  G/vq  is  of  the  same  general  order  of  magni¬ 
tude  as  cL(cL  =  cN  -  Cp)  and  cM,  as  we  have  assumed  in  this  development,  then 
the  angular  deviations  for  the  boosted  rockets,  due  to  unit  amounts  of  initial 
cross-spin  and  initial  yaw,  become  comparable  in  magnitude  with  the  corres¬ 
ponding  deviations  for  shells. 


The  Effects  of  w  ,  M  ,  and  F 
- c5  — c’  - c 

Parameters 


Dressed  in  Terms  of  Equivalent  Launch 


1.  Cross-wind,  w  :  (with  cross-wind  velocity  w  assumed  constant) 

_  Q  __  .  C . 

In  equations  (A.®  we  set  M  =  F  =  0  and  remove  the  gravity  term. 

c  c 


13U 


Since  wc  is  constant,  we  may  write 


and  thus 


(vA) '  =  (vA  +  w  ) 1  , 
c 


-  (2iqn  -  c^)#  +  c^CvA  +  wc)  = 


#  -  (vA  +  w  ) '  -  c..(vA  +  w  )  =  0 
c  N  c 


(A.3U) 


This  can  be  considered  as  the  homogeneous  system  (A. 6)  in  #  and  (vA  +  wc). 
Taking  4>o  =  $  =  ~  we  consiher  the  solution  of  (A.3i+)  in  the  frame¬ 

work  of  (A. 7),  (A. 8),  and  (A. 9),  with  A  replaced  by  (  A  +  w  /v)  and  with 

c 

the  corresponding  initial  value,  w  /v  .  Thus,  using  a  subscript  w  to  denote 

c  o 

the  solutions  due  to  w  in  the  absence  of  other  disturbing  factors,  one 
•  c 

obtains,  with  $  =  0  and  (A  +  w  /v  )  =  w  h  . 

*  O  c  O  O  C'  O* 


S  =  —  $K(s  , s) , 
w  v  5  o*  * 


(Aw  +  wc/v)  =  ^ll  +  A6(so,s)]J, 
o 


and  thus 


$  =—  9.  , 

W  V  o  3 

0 


w  w 

A  =  -£.  ti  +  A  ]  -  -2.  ,  . 

w  v  v  6J  v  3 
o 


(A.  35) 


w  w 

©  «  $  .A  =  -  —  [l  -  0j  +  — 

www  v  5  v 

o 


Note  that  the  angular  deviation  can  be  written  in  the  form. 


w  v 

®  =-^n--£  -aj. 


(A.  36) 


Except  for  the  assumption  that  w  remains  constant,  this  result  is  valid 

c 
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for  rockets  and  projectiles  in  general,  fin-  or  spin-stabilized  sp  long 

as  w  A  remains  sufficiently  small  to  be  considered  in  the  linear,  "small 
c  o 

yaw"  equations.  The  behavior  of  ®g,  of  course,  will  depend  upon  the 
characteristics  of  the  particular  projectile  in  question. 

Restrictive  Assumptions  to  be  made  in  the  considerations  to  follow: 


|  vA2|  ^  0.005; 


(A.  37) 


'  10cM  =  n2  =  0.25  \ 

j  n1 j  =  0.0005  ) 


n"  ~  10^ 

2.  Fin  Misalignment,  . 

From  equations  (5.1.5)  and  (A.l),  we  consider 


37 


-  (2iqn  -  cH)$  +  cm(vA)  =  c^v  e  , 
$  -  (vA)1  -  Cjj(t4>  =0. 


(A. 38) 


(A. 39) 


We  consider  first  the  case  where  spin  is  present.  For  convenience,  let 


H  =  c#c  T  6 


ip 


(A.liO) 


and  note  that 

H*  =  (i  n  +  v/v2)H.  (A.lpL) 

If  one  substitutes  in  the  left-hand  members  of  equations  (A. 39) 

#  = - —  (1  -  i  n'/n2],  (A.l*2) 

(l-2q)n  -cM 


vA  = 


-H 


(l-2q)n  -cM 


3n» 


n 


-)],  (A.U3) 


and 
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he  finds  that  the  left-hand  member  of  the  first  equation  differs  from  H 
by  a  quantity  of  the  order  of 


-  i  ti  (CH  *  ±)  ♦-?(=„- ^i)]  H, 


3  N  n 


which  is,  in  magnitude,  less  than  10°/o  of  |  H  |  .  Similarly,  the  left  member 

icN 

of  the  second  equation  reduces  essentially  to  — j(c^  -  3nr  /40H,  This  could 

•  H 

be  readily  absorbed  as  a  modification  of  the  §  given  in  (A.l|2),  but  it  is 
negligible  in  comparison  with  the  basic  #  of  (A- 2+2 ) . 

Thus,  to  a  good  approximation,  we  can  replace  the  n on-homogeneous  system 
(A. 39)  by  the  homogeneous  system, 


in  which 


(4  +  P)'  -  (2iqn-c„)($  +  P)  +  cM(vA  +  Q)  =  0, 


(S  +  P)  -  (vA  +  Q)»  -  cn(vA  +  Q)  =  0 


inH  .  ,  ,  21 

- 5 -  11  -  in'/n  ], 

(l-2q)n  -  cM 


(AM 


(AM 


(l-2q)n  -  c. 


(A.U6) 


Taking  =  Aq  =  0,  we  express  the  solution  of  (A. 14*)  in  the  framework 

of  (A. 7),  (A. 8),  and  (A. 9),  considering  §  replaced  by  ($  +  P)  with  its 
initial  value,  Pq,  and  A  replaced  by  (A  +  Q/v)  with  initial  value  Q0/v0» 
Thus,  using  the  subscript  p.  to  denote  the  fin-misalignment  effects,  we 
obtain 

i  ♦  p  -  p  i  +  (q  /v  )4_. 

ft  o  q  o'  o  6 

To  a  good  approximation,  since  (i  n  v  e1T1)  =  d(e3ri)/dt. 


f  P  dt  **  —  ] 

J+  nr 


-1H  .t  _  -IP  [  o 

nv  Jt  cd  co 
o  o 


(A.l*7) 
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Thus 

$  =  iP/co  -iP/co  +  P  g  +  (Q  /v  )3C,  (A.lj.8) 

p  '  o'  o  o  q  o'  0  6' 

and  similarly 

A  =  -Q/v  +  P  A  +  (Q  /v  )(1  +  Ac).  (A.U9) 

J-*  □  Cj  o  o  u 

Since  0  =  4>  -  A  ,  one  obtains 

P  P  PJ 

0  =  P0  +  (Q  /v  )  @  +  [Q/v  -  Q  /v  1  +  i[P/o>  -  P  co  ].  (A. 50) 

p  0  q  0  00  00  00 

Equation  (A.5>0)  is  the  simplest  in  the  case  where  co  =  n?  with  n  constant. 
Taking  n'  =0  and  (1  +  ic^/n)  =  1,  one  notes  that 

iP/co  =  -  Q/v.  (A.  5l) 

Thus,  where  the  spin  rate  is  approximately  proportional  to  the  velocity. 


0 

P 


1  n  VM^c 

2 

n  CM 


0  + 
q 


c  p 
me 


2 

n  -  c 


M 


(A. 52) 


In  (A. 52),  (l-2q)  has  been  replaced  by  1  for  convenience. 

It  should  be  noted  that  if  one  is  primarily  interested  in  an  estimate 
of  the  magnitude  of  the  angular  deviation  0  and  in  those  parameters  which 
significantly  affect®  ,  then  equation  (A. 52)  is  relevant  to  all  cases 

r 

where  spin  is  present,  provided  the  spin-pattern  satisfies  the  general 
restrictive  conditions  as  given  in  (A. 38).  This  can  be  recognized  by 
noting  that,  in  (A.lj.5)  and  (A.I46),  the  imaginary  terms  within  the  brackets 
do  not  affect  the  magnitudes  of  P  and  Q  but  rather  correspond  only  to 
slight  rotations  in  the  complex  plane.  Thus,  in  regard  to  magnitudes, 
relation  (A.5l)  is  a  good  approximation  in  general  and  likewise  (A. 52) 
with  n  replaced  by  n^  in  case  n  is  variable. 

Where  the  magnitude  of  ®  is  the  primary  concern,  we  can  also  use  the 
approximation  given  by  equation  (A. 32)  and  replace  (A. 52)  by  the  simpler 
result 


in  v  c„(l»2q) 

•A--**4 - •- 


an  v  c„ 
0  o  M 


e 


n  -  On, 
0  M 


'  2  q 

no  -•CM 


(A. 53) 
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Since  we  have  assumed  that  n  =  10c...  we  can  note  further  that,  since  the 

o  M 

basic  estimate  for®  ,  as  given  in  (A. 31)  or  (A. 33)*  is  independent  of:  n, 
®^/Pc  is  essentially  inversely  proportional  to  nQ.  As  a  first  order 


approximation,  we  can  write 


g  sin  e 


8  A  nr  (ol  +  -5 - r-2  > 


a.a) 


Thus  ®  /p  is  predominantly  imaginary  and  inversely  proportional  to  n  . 

V  c  0 

This  general  behavior  of  0  /p  also  shows  up  in  the  results  for  unboosted 

c 

firmer  rockets  with  slow  spin,  as  may  be  seen  in  Figure  20  of  [P]  which 
relates  to  unboosted  rockets. 

Although  we  are  not  primarily  concerned  with  the  no  spin  case,  it  is 
of  interest  to  note  that  our  approximation  methods  are  readily  applicable 
to  that  case  also.  With  n  =  0  and  the  exponential  e^  deleted,  equations 
(A. 39)  become 


+  cR$  +  cH(vA)  =  CflM^v 


(A.  33) 


^  -  (vA)«  -  Cjj(vA)  =0  j  .  (A.5 

In  this  case,  we  can  simply  absorb  the  forcing  term  as  a  (vA)  and  write 

•  • 

+  cR$  +  Cjjlv(  A  -  Pc)]  =0, 

*  -  M  A  -  Pc)]'  -  cN[v(  A  -  Pc)]  =  P-cv(cN  +  -&A2)  . 

Absorbing  the  right  member  of  the  second  of  these  equations  as  a  #,  one 
obtains 

[  t$  -  P“CV(CN  +  *A2)]*  +  CHI^  "  ^ct(cn  +  V"v2)]  +  cM[v(  A  -  Pc)]  -  0, 


-  Pcv(cN  +  t/v2)]  -  [v(  A  -  pc)]  -  cN[v(  A  -  pc)]  =  0, 


(A. 36) 


where  the  right  member  of  the  first  of  these  equations  should  be  essentially 


-P-cv(cN  +  Vv2)(ch  +  v/v2). 
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However-,  the  ratio  of  this  to  the  original  right  member  in  (A. 55)  is 

— (cM  +  t/V2)(cu  +  v/V2) 


and  this  is  negligible  in  comparison  with  1. 

2 

Treating  (c^  +  v/v  )  as  a  constant,  one  obtains  as  basic  estimates  in 
the  case  of  no  spin, 

*  2  2  ' 

%  "  >V(oN  *  )  -  Pc'o(ol  +  )#q'  “cV 

%  ’  (°N  *  Vv2)(s-So)  -  novo(o1}  *  Vv2)®.-  “0*6r 


V  *  %  -  Vo(cN  *  *A  )  4q  -  H0(l  *  »6); 


and  thus 

9 A  *  (cN  *  Vv2)(3-S0)  -  vo(dn  +  Wl  ®q  -  ®6-  {i>57). 

The  significant  term  in  (A. 57)  is,  of  course,  the  first,  which  increases 
with  trajectory  distance  (s-sq). 

3.  Linear  Thrust  Misalignment,  L_,  and  Angular  Thrust  Misalignment,  a  . 

f  c  ^  c 

We  shall  deal  with  these  effects  briefly  to  obtain  the  approximations 
primarily  for  the  case  where  co  «  nv  with  n  constant.  However,  the  same 
considerations  as  above  would  show  that  the  results  are  relevant  in  general 
(within  the  scope  of  the  stated  restrictive  assumptions)  to  other  spin 
patterns,  if  one  is  primarily  concerned  with  the  magnitude  of  the  resulting 
deviations. 

In  the  case  of  linear  thrust  misalignment,  one  has  the  basic  equations 


-  (2iqn  -  cH)#  +  cm(vA)  =  -  e1T1  , 

k  v 


$  -  (vA)1  -  Cjj(vA)  =  0. 


(A. 58) 


Noting  that 


£0  1W 


we  absorb  the  right  member  of  the  first  equation  of  (A. 58)  as  a  #,  replacing 

#  by  [3>  -  i(GLc/K  co)e11'] .  To  balance  the  second  equation,  we  replace  (vA  ) 

by  v[  A  -  (GL  /k2co2)eir|] . 
c 

One  obtains  an  approximate  homogeneous  system  and  thus  the  approximate 
results'  -GL  GL  GL  1  ip 

$T  =  -O-2-  *  -  -o  °o  **  ~  -o “(  O  ~  T)  J 

kco  ^  k  co  k  co  co 


At  =  — 5——  A  -  -5-m 
L  .2  q  .  2 

k  co  ^  k  co 

0  o 


GLeiT1 

h  (1  +  V  +  -hr  * 


and  thus 


®l/Lc  ’  ~2  2^  “o®q  +  ®6]' 

k  co  M 


(A.59) 


Using  (A. 32)  to  express  Sg  in  terms  of  0  ,  one  obtains 


slAc  -  f- 

k  co 


(1  -  2q)  ®  . 

q 


(A. 60) 


Similarly,  for  angular  thrust  misalignment,  one  can  obtain  as  good 
first  approximations,  _<iGa 


4-  » - —  $•, 

a  v  co  6’ 
o  0 


and  thus 


-iGa  iGa 

A  - - — (1  +  Ar)  + - —  e1T1  : 

a  v  co  '  n  t  co  ’ 

o  o 


0 Ja  =  0K  +  -^-(l  -  eiT1). 


a'  c  v  co  5  v  co  v 
00  00 


(A.  61) 


For  boosted  rockets  0g  will  be  small  in  comparison  with  1  and  the  behavior 


of  ®a/ac  will  be  determined  primarily  by 


.  _  Y  CO 

lG  (  1  _  -.0  0  JWh 

V  co  ^  vto 

O  0 


If  (vzo)  increases  during  burning,  the  effect  of  the  rotating  term  will 
damp  out.  Otherwise  the  magnitude  can  vary  approximately  from  0  to  SG/vko^* 


lUi 


Summary 

Under  the  assumption  of  a  spin  pattern  within  the  scope  of  restrictive 
conditions  (A. 38)  and  with  v- =  500  ft. /sec.  and  v/v  =  0.001  (ft.  ),  there 
have  been  obtained  the  asymptotic  estimates, 


and 


®  « (c„  +  v  /v  ) /cMv  , 

q  N  o  o  77  M  o* 


8  -  21,*“oV 


together  with  the  approximations, 

8A  -  -V  Crf“oSq  < 


s./L  a 

L7  c  ,  2  q* 

k  co  ^ 


0  /a  *  J£L  (1  _  eJ‘r|  ). 

a7  c  v  co  v  vco 
o  o 


v  co 

-Q  o  a1Tl 


(A. 62) 

(A. 63) 

(A. 61;) 

(A. 65) 

(A. 66) 


As  the  combined  angular  deviation  (with  the  gravity  effect  excluded),  one 
can  write 

* 


=  0  +  ($  + 
o  o 


n 


iGL  w 

Q  +  (  A  +  -S.  )  0 
,2a  o  v  o 

k  co  ^  o 


w  v  iGa 

-  — (1  -  — )  + - £(1 

V  V  V  co 

o  0  o 


V  CO 
0  0 

vco 


eiTl 


). 


(A. 67) 


l?!ram  (A. 62)  and  (A. 63),  one  notes  thatioughly  0^ 

l/v  ,  and  ©_  of  the  order  of  (2iqn  ). 

7  o’  6  ^  0 


is  of  the  order  of  magnitude 


APPENDIX  B 


SOME  CONSIDERATIONS  CONCERNING  THE  EFFECT  OF 

• 

LAUNCHER  MOTION  UPON  THE  LAUNCH  PARAMETER,  9 

o 

The  simple  model  to  be  considered  below  was  formulated  in  an  attempt 
to  obtain  some  order-of -magnitude  estimates  of  possible  effects  of  launcher 
motion  upon  #  .  The  situation  under  consideration  involves  a  boosted 
spin-stabilized  rocket  launched  from  a  rifled  tube.  Experimental  data, 
secured  at  Redstone  Arsenal  and  at  Aberdeen  Proving  Ground,  have  exhibited 
high  frequency  oscillatory  motion  of  the  muzzle  of  the  launcher  in  both 
the  vertical  and  horizontal  planes. 

The  model  is  to  be  considered  in  the  following  context.  Presumably 
the  rocket  is  well  constrained  by  the  rifled  launcher  during  most  of  the 
interval  of  emergence  from  the  launcher,  and  the  launcher  motion  is 
transmitted  to  the  rocket  as  a  whole,  with  no  freedom  for  rotational 
motion  on  the  part  of  the  rocket  about  the  end  of  the  launcher.  However, 
suppose  that  the  rocket  and  the  launcher  are  so  constructed  that  there 
can  be  a  short  length  at  the  rear  of  the  rocket  (less  than  an  inch, 
perhaps)  such  that  one  might  expect  continued  contact  (of  some  sort) 
between  launcher  and  rocket  as  the  rear  emerges  but  with  sufficient 
clearance  (or  partial  freedom  from  constraint)  so  that  the  rocket  is 
free  to  rotate  about  a  point  at  the  muzzle  of  the  launcher.  In  such 
a  situation,  high  frequency  oscillations  of  the  muzzle  might  have  a 
significant  effect  upon  #  . 

The  considerations  below  are  limited  to  a  plane  and  the  vertical 
plane  is  specifically  considered.  However,  the  results,  except  for  the 
gravity  effect,  would  be  applicable  to  the  horizontal  plane.  The  dia¬ 
grams  below  (some  of  them  exaggerated)  indicate  the  geometrical 


Hi  2 


framework  to  be  considered. 


1.  Let  y  denote  the  displacement  vector  of  the  muzzle  and  £  the  position 
vector  of  the  c.g.  Positive  displacements  will  be  measured  upward. 

2.  Let  x  denote  the  horizontal  distance  from  the  muzzle  to  the  c.g. 
(positive  to  the  right). 

3.  Let  #  denote  the  orientation  angle  of  the  rocket  axis  relative  to  the 
reference  axis  (counter-clockwise  angles  positive). 


H|l| 


X 


For  small  angles  §, 

x  <£  =  £  -  y.  (B.l) 

To  obtain  equations  of  motion,  we  consider  the  motion  of  the  c.g. 
at  right  angles  to  the  reference  axis  and  the  rotational  motion  of  the 
rocket  about  its  c.g.  The  forces  and  torques  to  be  considered  are  those 
due  to  gravity  and  to  the  reaction  force,  R,  acting  upon  the  rocket  at 
the  muzzle.  This  reaction  force  is  assumed  to  be  at  right  angles  to  the 
reference  axis. 


Reference  Axis 


Thus 


1 h$ 


I  m  X  -  -  mg  +  E 

2  ** 

ink  §  =  -  xR 

2  * 

in  which  ink  denotes  the  transverse  moment  of  inertia  (about  a  transverse 
axis  through  the  c.g.). 

From  eq.  (B.l), 

K  =  x  #  +  y. 

Over  the  short-time  interval  under  consideration,  we  assume  that 


Then 


and 


x  =  v  =  a  constant,  v  ,  the  launch  velocity. 


£  =  x  $  +2 v  9  +  y  , 

o 


x  9  +  2vq#  =  -  g  -  y  +  R/m 


2  '  ’ 
k  $ 


Eliminating  R,  we  obtain 


*  -  xR/m 


(B.3) 


(B.U) 


Q  O'*  *  *  * 

(x  +  k  j  $  +  2vqx  #  =  -gx  -  xy  .  (B.f>) 

In  the  derivation  of  the  above  equations,  it  is  evident  that  no 
attempt  has  been  made  to  take  account  of  the  rocket's  spinning  motion 
(about  its  longitudinal  axis).  If  the  gyroscopic  effects  of  the  spin  were 
included,  there  would  be  a  gyroscopic  tGrque  to  rotate  the  rocket  out  of 
the  plane  under  consideration.  In  notation,  relative  to  a  complex  plane 
perpendicular  to  the  reference  axis,  (as  in  the  3-dimensional  equations 
of  motion  of  spin-stabilized  rockets),  this  gyroscopic  term  would  appear 
in  eq.(B.5)  as  (-2iqco  ^  in  which  co  denotes  the  spin  rate  and  2q  the 
ratio  of  the  moments  of  inertia.  With  2q  ~0.0lj.  and  co  =  nvQ(rad./sec, ) 


llj.6 


with  l<n  <2,  the  coefficient  of  $  in  the  gyroscopic  term  would  have  a 
magnitude  lees  than  0.08vq,  whereas  the  term  in  (B.5>),  2vqx  3>,  with 
x  ~1  ft.,  would  have  a  coefficient  (of  $)  ~  2vq.  The  gyroscopic  effect  is 
much  smaller  and,  over  the  short  time  Interval  involved,  the  gyroscopic 
effect  would  be  insignificant,  particularly  if,  at  the  start  of  the  interval, 
no  appreciable  $  had  as  yet  developed. 

Equation  (B.fJ)  is  readily  integrated.  Since  we  are  particularly 
interested  in  the  effect  over  some  Ax  distance,  as  the  rear  of  the  rocket 
emerges,  we  transform  (B.f>)  into  the  equation 


(x2  +  k2)-g+2xi=-f  x-f  V. 

O  0 


Thus, 


(B.6) 

(B.  7  ) 


Note  that  equation  (B.6)  is  linear  and  hence  that  we  can  consider 
separately  the  effects  of  gravity,  the  initial  $  (at  the  start  of  this  r ' 
interval  under  consideration),  and  the  launcher  motion. 

For  boosted  launch  velocities,  vq  =  £00  ft. /sec.  perhaps,  the  gravity 
effect  is  quite  negligible.  Writing  x^  =  x^  +  Ax,  with  x^~l  ft.  and 
x  ~1  inch  or  less. 


x,  +  x. 


(*22-  *2>V  -  if  (I22-  *12)  ■  -  f 


~  -  0.00£  (or  less  in  magnitude). 

With  k2  ~0.5  ft.2, 

|  |~  O.OOI4.  or  less.  (rad. /sec.) 

4  4 

For  the  effect  of  $^(the  initial  $),  we  have 


31*7 


.  x  2+  k2  .  .  (x  2-  x  2)  . 

&  =  —± - §  =  # _ _ _ _ —  # 

£  2  2  X  X  2  2  1 


x2  +  k 


x2  +  k 


?x  *  +  xo 

=  ^x - 2 - 2  \*  W±th  X  =  — 2 -  * 

x2  k 

<  2x 

Note  that  for  Ax  =  1  inch,  — - - -  Ax  ~  0.1  or  less.  Thus  the  resulting 

x  +  k 

•  C  i 

$2  is  of  the  same  order  of  magnitude  as  being  slightly  less  in  magnitude. 
If  the  constraints,  prior  to  this  final  Ax  of  emergence,  are  such  as  were  pic¬ 
tured  at  the  start  of  this  appendix  namely  such  that  the  would  essentially 
represent  a  rotational  rate  of  the  launcher  tube  itself  (or  of  the  more 
flexible  portion  of  the  tube  near  the  muzzle),  we  are  here  also  presented 
with  the  question  as  to  whether  the  launcher  motion  itself  could  involve 
a  significant 

Turning  our  attention  to  the  effect  of  y,  we  consider 


d  r  /_2  .  ,2,1,  _  x  ‘  ‘ 


o 


(B.3) 


with  =  0.  Since 


J  ^2  dx  =  Jx  d(y)=  xy  -  J y  dx 


=  xy  -  VQy  , 


we  obtain 


(x22+  k2)$2  =  -  £*7  -  VQy]x2 


(B.9) 


If  the  displacement  of  the  muzzle  is  represented  by 
y  =  A  +  Bt  +  C  sin(ht  +  a), 


(B.10) 


li|8 

i.e.,  a  sinusoidal  oscillation  superimposed  upon  a  linear  displacement, 

(A  +  Bt),  with  t  =  0  corresponding  to  x  =  0,  then  we  note  that 

xy  -  vQy  =  fix  +  Chx  cos(ht  +  a)  -  Avq-  B/Qt  -  CvQSin(ht  +  a) 

=  -  Avq  +  CvQ[ht  cos(ht  +  a)  -  sin(ht  +  a)] 

=  Cv  ^  +  0  cos  (0  +  a)  -  sin(0  +  a)]  , 

O  L» 

in  which  0  denotes  ht. 

Thus 

0 

(x22+  k2)*2  =  -  Cvq  [-£  +  f(e)  e2  ,  (B.ll) 

where  -  - 

f  (0)  =  0  cos(0  +  a)  -  sin(0  +  a). 

It  is  clear  that  the  term  (-A/C)  does  not  contribute,  and-that  one  is  ' 
concerned  with  the  change  in  f(0).  The  magnitude  (and  algebraic  sign) 
of  the  change  in  f(0)  depends  upon  the  location  of  the  interval, 

A0  =  @2  _  and  also  upon  the  phase  angle,  a. 

One  can  quickly  verify  that  f(8)  has  successive  maxima  and  minima  at 
9  =  0,  it  -  a,  2 it  -  a,  3tt  -  a,  etc.,  with  the  graphical  representation  of 
f(0)  somewhat  as  follows: 
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It  is  clear  from, the  above  curve  that  where  the  original  sinusoidal 
contribution  to  the  displacement  is  near  zero  (i.e.  sin(e  +  a)  »  0),  the 
corresponding  change  in  f  (0)  (for  A0  not  too  large)  would  not  be  as 
serious  as  Af  (  8)  might  be  where  J"sin(e  +  a) j  ~  1. 

If  one  were  to  consider  Ax  and  the  corresponding  A0  to  be  sufficiently 
small  so  that  a  differential  approximation  might  be  used,  one  can  return 
to  equation  (B.8)  arid  consider 

•  X 

A[(x2+  k2)$]  »  -  ^  y^  Ax.  (B.12) 

o 

Since  A[(x2+  k2)$]  *  (x.^2+  k2)  A§  +  A(x2+  k2),  then  with  ^  =  0, 

(xx2+  k2)  A 4  =  (Xl2+  k2)$2  *  -  ^  y^ Ax., 

o 

and  with 

y  =  A  +  Bt  +  C  sin(ht  +  a), 

2  2 
y  =  -  Ch  sin(ht  +  a)  =  -  Ch  sin(0  +  a), 

one  obtains 

P  2  •  Ch  X-. 

(x-j^  +  k  )$2  *  — — —  sin(  0  2  +  a)  Ax. 
o 

For  some  of  the  relevant  data,  C~  0.002  fr.,  h~2600  rad. /sec., 

2 

x,  ~1  ft.,  k  ~  0.£,  and  v  ~  $00  ft. /sec.  If  one  restricts  Ax  to  l/2  inch, 
h  ^ 

then  A0  =  —  A  x~  0.22  rad.  and  the  above  approximation  is  reasonably  good. 
One  finds  tfiat 

^2  ~(3A)sin(6  j  +  ra(i./sec. 

Note  that,  within  the  scope  of  this  differential  approximation,  the 
resulting  contribution  to  $  is  proportional  to  Ax  and  depends  significantly 
both  upon  the  frequency  and  amplitude  of  the  sinusoidal  oscillation  and 
also  upon  the  particular  portion  of  the  sinusoid  involved. 
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